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1.  INTRODUCTION 


The  response  of  soil  to  high  amplitude ,  short 
duration,  impulse  loads  characteristic  of  the  airblast  and 
ground  motion  from  a  weapon  explosion  is  an  important 
problem  in  protective  construction  design  and  analysis 
(Calhoun  and  Kraft,  1966;  Schindler,  1968).  A  means  of 
understanding  the  soil  response  is  through  laboratory 
investigation.  However,  the  laboratory  environment  must  be 
able  to  reflect  the  type  of  confinement,  magnitude  of 
stress  change,  and  the  time  scale  of  loading  expected  in 
the  problem  (Whitman,  1970). 

A  reasonable  assumption  is  that  the  soil  response  in 
the  superseismic  region  of  an  airblast  will  be 
approximately  one-dimensional  (Crawford,  Higgins,  and 
Bultmann,  1974).  The  time  scale  of  loading  is  essentially 
instantaneous  (Crawford  et  al.,  1974).  Peak  stresses  can 
be  in  excess  of  3  000  MPa  (Brode,  1984).  To  model  the  soil 
response  to  this  type  of  loading,  uniaxial  strain  devices 
have  been  developed  (e.g.,  Schindler,  1968).  Current 
devices  can  apply  stresses  up  to  400  MPa  with  loading 
tirues  on  the  order  of  0.3  milliseconds  (e.g.,  Jackson, 
Ehrgott,  and  Rohani,  1980). 

Over  the  last  thirty-five  years  the  split-Hopkinson 


presura  bar  (SHPB)  technique  has  been  used  as  a  tool  for 
investigating  the  response  of  metals,  rocks,  ceramics, 
foams,  and  other  materials  to  short  duration  compressive 
impulse  loads  (e.g.,  Lindholm,  1964;  Hodge  and  Wasley, 
1969;  Christensen,  Swanson,  and  Brown,  1972).  Some  of  the 
SHPB  devices  in  use  can  apply  stresses  in  excess  of  1000 
MPa  with  loading  times  on  the  order  of  0.04  milliseconds 
(Gaffney  and  Brown,  1984).  Until  recently,  the  SHPB 
technique  has  not  been  readily  applied  to  the  field  of 
soil  mechanics.  Because  of  the  ability  to  apply  high 
stresses  at  a  high  rate  of  loading,  the  adaption  of  the 
SHPB  technique  to  measure  the  dynamic  response  of  soil 
seems  to  be  a  natural  extension.  The  objective  of  this 
research  has  been  to  determine  whether  the  SHPB  technique 
can  be  effectively  used  to  measure  the  dynamic  response  of 
soil. 

The  use  of  soil  as  specimens  in  a  SHPB  experiment  is 
not  a  trivial  matter  because  soils  hava  very  low  wave 
speeds  (  =  300  m/s)  in  comparison  to  the  traditional 
materials  tested  in  the  SHPB  (e.g.,  steel,  5000  m/s). 
Soils,  also  exhibit  nonlinear  hysteretic  behavior  which 
will  cause  the  amplitude  of  a  stress  wave  to  attenuate  as 
it  propagates  through  it  (Hendron  and  Auld,  1968).  In 
addition,  the  relatively  low  unconfined  compressive 
strength  of  the  soil  (e.g.,  <  0.1  MPa)  creates 
difficulties  in  controlling  boundary  conditions. 
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This  dissertation  is  organized  into  seven  chapters. 
Chapter  2  presents  the  development  of  pressure  bar 
experimentation  from  its  inception  by  Hopkinson  (1914), 
through  the  contributions  of  Davies  (1948),  and  the 
introduction  of  the  SHPB  method  by  Kolsky  (1949).  Also 
presented  is  the  basic  theory  of  measurement  user*  in 
reducing  the  experimental  data,  as  well  as  contributions 
to  the  method  made  by  recent  investigators,  and  the 
methods  application  to  soil  mechanics. 

Chapter  3  describes  the  SHPB  apparatus,  the 
characteristics  of  the  soil,  and  how  the  specimens  were 
prepared.  The  assumptions  of  the  experiment  and  their 
bearing  on  experimental  results  as  well  as  the  data 
reduction  procedure  have  been  addressed  in  chapter  4. 
Chapter  5  presents  the  experimental  results  in  terms  of 
the  soil  stress-strain  response,  and  establishes  that  the 
response  is  nearly  one  of  uniaxial  strain.  In  addition, 
the  ability  to  replicate  experimental  results  is 
demonstrated  and  the  observed  soil  behavior  at  different 
compaction  conditions  to  a  range  of  applied  stresses  are 
presented.  Chapter  6  establishes  that  the  experimental 
assumptions  can  indeed  be  satisfied  when  using  soil 
specimens  in  a  SHPB  experiment.  The  results  presented  in 
chapter  5  are  also  analyzed  to  evaluate  the  dominant 
parameter  governing  the  observed  stress-strain  response 
and  a  strain-rate  independence  is  established  for  specimen 
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strains  less  than  the  initial  gas  porosity.  In  addition, 
specimen  conditions  that  may  lead  to  erroneous 
stress-strain  response  and  an  evaluation  o£  the  SHPB 
experiment  are  presented.  Chapter  7  presents  the 
conclusions  drawn  £rom  the  analyses  and  provides 
recommendations  for  further  research.  Ten  appendices  have 
been  included  to  supplement  and  support  the  discussions  in 
the  main  text. 


2.  DEVELOPMENT  OP  THE  SPL1T-HOPKINSON 
PRESSURE  BAR  EXPERIMENT 

2.1  The  Hopkinson  Pressure  Bar 

In  1914,  Bertram  Hopkinson  devised  a  method  to 
experimentally  study  the  mechanical  behavior  of  solids  in 
response  to  short  duration  impulse  loads.  The  main 
components  of  the  apparatus  constructed  by  Hopkinson  were 
a  long  cylindrical  steel  bar,  a  time  piece,  and  a 
ballistic  pendulum  (see  figure  2.1).  As  shown  in  figure 
2.1,  the  cylindrical  steel  bar  was  suspended  by  threads, 
such  that  it  was  free  to  swing  in  a  vertical  plane  and 
remain  parallel  to  its  original  position.  The  time  piece 
was  the  same  diameter  as  the  cylindrical  bar  and  attached 
to  it  by  magnetic  attraction. 

The  experiment  was  initiated  by  the  impact  of  a  lead 
bullet  or  the  detonation  of  gun  cotton,  near  the  end  of 
the  cylindrical  bar  opposite  the  end  to  which  the  time 
piece  was  attached.  The  impact  produced  a  compressive 
stress  wave  that  was  assumed  to  propagate  down  the  length 
of  the  bar  and  through  the  joint  between  the  bar  and  time 
piece  without  distortion.  The  stress  distribution  over  the 
cross  section  of  the  bar  was  assumed  to  be  uniform.  Upon 
reaching  the  end  of  the  time  piece,  the  compressive  stress 


wave  waa  assumed  to  reflect  as  a  tensile  wave.  The  joint 
between  the  bar  and  the  time  piece  could  not  sustain  a 
tensile  force  greater  than  the  magnetic  force  connecting 
them.  Hence,  when  the  tensile  force  across  the  joint 
exceeded  the  magnetic  force,  the  time  piece  would  separate 
from  the  bar  with  a  given  amount  of  momentum.  This 
momentum  was  measured  by  the  ballistic  pendulum.  The 
momentum  trapped  in  a  given  length  time  piece  corresponded 
to  the  area  under  the  stress-time  curve  between  points  of 
equal  pressure,  separated  by  the  time  required  for  the 
wave  to  travel  the  length  of  the  time  piece  twice.  By 
measuring  the  momentum  delivered  to  different  length  time 
pieces,  the  area  under  the  stress-time  curve  for  different 
time  intervals  could  be  obtained  allowing  a  complete 
stress-time  curve  to  be  constructed.  However,  the  precise 
form  of  the  stress-time  curve  could  not  be  ascertained. 
This  is  because  the  commencement  of  the  different 
intervals  was  not  known.  Although  the  precise  form  of  the 
stress-time  curve  could  not  be  determined,  the  maximum 
stress  and  the  total  duration  of  impact  could  always  be 
evaluated. 

Despite  its  innovative  approach,  Hopkinson's  method 
suffered  from  several  experimental,  as  well  as 
theoretical,  limitations.  Experimentally,  the  force 
necessary  to  separate  the  bar  and  time  piece  was  unknown. 
Also,  the  available  intrumentation  could  not  detect  the 
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commencement  of  the  time  interval  over  which  the  momentum 
was  measured.  Theoretically ,  it  was  assumed  that  the 
applied  stress  would  be  uniform  over  the  cross  section  of 
the  bar  and  that  the  stress  wave  would  propagate  down  the 
bar  without  distortion.  These  limitations  were  addressed 
and  overcome  through  an  extensive  study  of  the  Hopkioson 
bar  method  by  Davies  (1948). 

2.2  The  Davies  Bar 

Davies  (1948)  made  several  changes  to  the  original 
experimental  method  developed  by  Hopkinson.  The  generated 
stress  wave  produced  both  longitudinal  and  radial 
displacements  as  it  propagated  down  the  bar.  Davies 
developed  a  means  of  measuring  the  displacements 
electrically  through  a  bar  condenser  unit  mounted  on  the 
free  end  of  the  pressure  bar  (see  figure  2.2).  The 
condenser  unit  was  initially  charged  to  a  high  voltage 
that  would  be  constant  for  small  time  intervals.  The 
displacements  caused  a  change  in  capacitance  of  the 
condenser  unit,  producing  an  electrical  signal  that  was 
displayed  as  a  function  of  time  on  a  cathode  ray 
oscillograph.  Using  elastic  wave  theory,  Davies  translated 
the  electrical  records  into  a  precise  stress-time  curve. 
The  electrical  recording  system  avoided  the  use  of  a  time 
piece,  and  alleviated  the  associated  experimental 
problems.  Therefore,  with  the  introduction  of  electrical 
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data  recording,  the  experimental  limitations  associated 
with  Hopkinson's  original  experimental  method  are  avoided. 

As  part  of  the  same  stud/,  Davies  (1948)  described 

the  phenomenon  of  wave  dispersion  and  established  the 

accuracy  of  the  experimental  results  when  assuming 

one-dimensional  wave  propagation  in  the  pressure  bar. 

One-dimensional  elastic  wave  theory  assumes  that  a 

propagating  wave  will  travel  without  change  of  form,  at  a 

constant  velocity  (Cq) ,  that  is  independent  of  wavelength 

( A ) •  The  constant  longitudinal  velocity  (0o)  can  be 
2 

t  jsed  as  Co  *  E/p  ,  where  E  is  Young's  modulus  and  p 
is  mass  density.  This  assumption  is  only  true  if  A  of  the 
propagating  wave  is  long  compared  to  the  radius  (R)  of  the 
bar.  The  equation  governing  the  longitudinal  vibrations  of 
a  infinitely  long  circular  cylinder  (referred  to  in  the 
literature  as  the  dispersion  equation),  developed 
independently  by  Pochhammer  (1876)  and  Chree  (1889),  shows 
that  the  one-dimensional  assumption  will  only  be 
approximate,  as  the  velocity  of  wave  propagation  will 
depend  on  A  . 

The  compressive  stress  wave  in  a  pressure  bar 
experiment  is  composed  of  many  frequencies.  Evidence  of 
this  will  be  shown  in  section  4.1.  Because  the  wave 
velocity  (Cn)  will  depend  on  A,  each  frequency  component 
will  travel  at  a  different  Cn»  With  some  frequency 
components  travelling  faster  than  others,  the  wave  will 
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change  form  or  disperse  as  it  progagates.  The  amount  a 
wave  disperses  will  effect  the  uniformity  of  stress 
distribution  over  the  cross  section  of  the  bar.  The  more  a 
wave  disperses,  the  less  uniform  the  stress  distribution. 
Measurements  in  the  experiment  are  made  at  the  surface  of 
the  pressure  bar  (see  figure  2.2).  Therefore,  if  the 
surface  measurements  do  not  represent  the  behavior  over 
the  cross  section  of  the  bar,  the  accuracy  of  the 
experiment  is  reduced. 

By  using  the  dispersion  equation,  Davies  (1948) 
showed  that  the  oscillations  in  the  recorded  displacement 
history  were  attributable  to  dispersion  and  not 
experimental  inaccuracies.  He  also  established  the  error 
in  measuring  the  displacement  history  at  the  surface  of 
the  pressure  bar  to  be  +  3  percent,  provided  that 
R/A  <  0.1. 


2.3  The  Kolskj 


With  the  experimental  and  theoretical  foundations  of 
the  method  established,  Kolsky  (1949)  modified  the 
apparatus  to  permit  dynamic  material  properties  to  be 
determined  by  indirect  measurments .  Kolsky  sandwiched  a 
thin  cylindrical  specimen  (i.e.,  having  an  aspect  ratio, 
1/d  =  0.1,  where  1  is  the  specimen  length,  and  d  is  the 
specimen  diameter '  between  two  pressure  bars.  The  pressure 
bars  were  fitted  with  condenser  microphones  for  data 
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recording.  He  then  applied  an  impact  load  such  that  a 
compressive  stress  wave  propagated  toward  the  specimen 
(see  figure  2.3).  When  the  propagating  stress  wave  reaches 
the  specimen,  a  portion  of  it  will  be  reflected  as  a 
'  tensile  wave  and  a  portion  will  be  transmitted  through  the 

specimen.  The  amplitudes  of  the  reflected  and  transmitted 
waves  will  depend  on  the  physical  properties  of  the 
I  specimen  as  well  as  those  of  the  pressure  bars.  By 

assuming  a  uniform  distribution  of  stress  and  strain  along 
the  longitudinal  axis  of  the  specimen,  Kolsky  was  able  to 

i 

I  develop  relationships  between  the  incident  and  transmitted 

I  displacements  (recorded  at  the  respective  condenser 

microphones  mounted  on  the  bars)  and  the  average  stress, 

I  strain,  and  strain-rate  in  the  specimen.  These  equations 

i 

will  be  derived  in  section  2.4.  This  experimental  method 
is  now  known  as  the  Kolsky  method,  or  the  split-Hopkinson 

I  pressure  bar  (SHPB)  method. 

• 

■  By  using  thin  specimens,  Kolsky  attempted  to  avoid 

i 

I 

;  the  complication  of  axial  inertia.  Radial  inertia  effects 

i  were  accounted  for  through  a  numerical  correction  that 

assumed  a  frictionless  specimen-bar  interface.  The  effects 
of  axial  and  radial  inertia  in  a  SHPB  experiment  are 
i  important  because  they  act  to  oppose  the  equilibration  of 

I 

j  stress  within  the  specimen.  Kolsky  determined  that  by 

i 

J  using  a  thin  layer  of  lubrication  between  the  specimen, 

I  and  the  bars,  the  specimen-bar  interfaces  could  be  assumed 

i 


13 


to  be  frictionless. 


2.4  Theory  of  Measurement 


Assuming  that  a  wave  generated  from  a  disturbance 
produced  in  an  isotropic  elastic  bar  propagates  parallel 
to  the  longitudinal  axis  of  the  bar  such  that  it  can  be 
considered  to  be  a  plane  wave,  the  general  form  of  the 
wave  equation  will  be  one-dimensional.  In  the 
one-dimensional  case,  the  elements  of  the  bar  will  extend 
and  contract  with  no  lateral  displacement  of  the  axis  of 
the  bar  (Kolsky,  1963).  A  one-dimensional  condition  is 
idealized  in  figure  2.4.  The  assumptions  of  the 
one-dimensional  condition  are:  (1)  each  plane  cross 
section  will  remain  plane  while  in  motion,  (2)  the  stress 
over  the  cross  section  will  be  uniform,  and  (3)  that  the 
wavelength  of  the  wave  will  be  long  compared  to  the  cross 
sectional  dimension  of  the  bar  (Kolsky,  1963). 

Isolating  a  small  element  of  the  bar,  dx,  with  cross 
sectional  area  A,  the  stress  at  section  A-A  will  be  o 
and  the  stress  at  section  B-B  will  be  a  +  (  do  /  d x  )  dx. 
Talcing  the  x-axis  as  the  direction  of  wave  propagation, 
the  components  of  Newton's  second  law  of  motion  (F  ■  m  a, 
where  F  is  force,  m  is  mass,  and  a  is  acceleration)  can  be 
written  as; 


F 


(2.1) 


V, 
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; 


P 


dx  A 


» 


m  ■  p  A  dx  , 


where  p  ia  the  mass  density  of  the  bar  and  U  is 
displacement.  Combining  terms  yields; 


dx  A  a 


Equation  (2.5)  reduces  to; 


(2.2) 


(2.3) 

(2.4) 


(2.5) 


(2.6) 


Using  Hooke's  law  for  an  isotropic  elastic  solid  in 
uniaxial  stress; 


O  ■  EC  ,  (2.7) 

and  rewriting  C  in  terms  of  displacement  as; 

-au.  ,  (2.8) 

dx 

equation  (2.7)  can  be  written  as; 

O  »  E_&i.  •  (2.9) 

dx 

Taking  the  partial  derivative  of  equation  (2.9)  with 
respect  to  x  yields; 
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o  -  e  a  u  _ 

x  3 


(2.10) 


Substituting  the  expression  fordo/ 3  x  from  equation 


(2.10)  into  equation  (2.6)  yields; 


aV-  JL  S2U 
3 12  p  d  x2 


(2.11) 


Substituting  Co  “  E/p  into  equation  (2.11)  provides  an 
expression  for  the  propagation  of  longitudinal  waves  of 


infinite  wavelength  along  a  bar; 


2  2  -2 

u  -  c0  3  u  . 

it3  3x2 


(2.12) 


The  solution  of  equation  (2.12)  is  attributed  to 
d'Alembert,  and  may  be  given  as  (Timoshenko  and  Goodier, 


1970) ; 


U  -  f(x  -  C0t)  +  g(x  +  Cot) 


(2.13) 


where  f  and  g  are  arbitrary  functions,  with  f  representing 
a  wave  travelling  in  the  positive  x-direction  and  g 
representing  a  wave  travelling  in  the  negative 


x-direction . 


If  a  wave  is  considered  to  travel  only  in  the 


negative  x-direction,  equation  (2.13)  can  be  written  as; 


U  -  g(x  +  Cot)  . 


(2.14) 


Differentiating  equation  (2.14)  with  respect  to  x  yields; 
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(2.15) 


where  the  prime  denotes  the  derivative  with  respect  to  the 
argument.  Differentiating  equation  (2.14)  with  resect  to  t 
yields; 


Co  g 


(2.16) 


Combining  equations  (2.15)  and  (2.16)  yields; 


au  -  Co  All  .  (2.17) 

dt  ax 


Substituting  Hooke's  law  (equation  2.9)  into  equation 

(2.17)  yields; 

du  -  CnO  .  (2.18) 

at  E 

2 

Substituting  pCo  for  E,  equation  (2.18)  can  be  rewritten 

as: 

3u  »  a  .  (2.19) 

(Ft  p  co 

By  substituting  v  for  the  particle  velocity  (  3u/  dt) , 
equation  (2.19)  simplifies  to; 

0  »  p  Co  v  .  (2.20) 

The  product  pCo  is  commonly  referred  to  as  the 
characteristic  impedance  of  the  material  (Rinehart,  1975). 
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The  procedure  followed  in  deriving  the  equation  for 
particle  velocity  is  similar  to  that  of  other  authors 
(e.g.,  Kolsky,  1963;  Rinehart,  1975;  Zukus ,  Nicholas, 
Greszczuk,  and  Curran,  1982). 

A  diagram  of  the  SHPB  near  the  specimen  is  shown  in 
figure  2.5.  The  stress  waves  Of  (incident  wave)  and  Or 
(reflected  wave)  act  on  interface  1  and  ot  (transmitted 
wave)  acts  on  interface  2.  Before  any  reflections  occur, 
the  particle  velocity  of  the  incident  bar  is  given  as; 

Vj  -  °i  .  (2.21) 

P  C0 

If  the  characteristic  impedance,  or  area  of  the  specimen, 
is  different  than  that  of  the  pressure  bars,  a  portion  of 
the  compressive  stress  wave  at  interface  1  will  be 
reflected  as  a  tensile  wave  and  that  portion  of  the  stress 
wave  which  the  specimen  is  able  to  support  is  transmitted 
through  the  specimen.  When  the  portion  of  the  stress  wave 
propagating  through  the  specimen  reaches  interface  2,  the 
wave  is  once  again  partitioned,  with  a  portion  being 
reflected  back  into  the  specimen  and  a  portion  being 
transmitted  into  the  transmitter  bar.  The  reflected  wave 
at  interface  2  is  compressive;  hence,  it  will  continue  to 
traverse  the  specimen,  increasing  in  amplitude  with  each 
transit.  The  characteristic  impedance  of  the  specimen, 
relative  to  the  characteristic  impedance  of  the  pressure 
bars,  and  any  difference  in  area  will  govern  the  increase 


21 


in  amplitude.  The  result  of  these  internal  reflections  is 
that  the  net  particle  velocity  of  interface  1  will 
increase,  and  then  decrease  with  time,  while  the  net 
stress  on  interface  1  will  decrease,  and  then  increase 
with  time.  The  net  particle  velocity  of  interface  2  will 
increase  with  time,  as  will  the  net  stress  on  interface  2. 

The  equations  for  the  reflection  and  transmission 
coefficients  for  a  wave  impinging  on  an  interface  are 
derived  in  figure  2.6.  Figure  2.7  illustrates  the  process 
of  elastic  wave  propagation  exercising  the  repeated 
application  of  the  reflection  and  transmission 
coefficients.  The  process  is  terminated  at  the  point  of 
stress  equilibrium  with  the  incident  wave.  This  figure  was 
prepared  with  the  assumption  that  the  wavelength  of  the 
incident  wave  is  infinite.  For  illustrative  purposes,  it 
was  also  assumed  that  the  density  of  the  specimen  was 
one-half  that  of  the  pressure  bars  and  the  wave  velocity 
(Co  )  in  the  specimen  and  pressure  bars  is  equal,  thus  the 
characteristic  impedance  of  the  specimen  is  one-half  that 
of  the  pressure  bars.  The  pressure  bars  are  considered  to 
be  of  a  sufficient  length  to  avoid  the  necessity  of 
considering  reflections  from  their  ends.  An  implicit 
assumption  for  using  equations  (4)  and  (5)  derived  in 
figure  2.6  is  that  da/de  is  constant  for  the  specimen 
material.  However,  when  a  specimen  deforms  plastically,  as 
does  the  soil  specimens  in  this  research,  da  /dc  is  not 
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At  interface  A  -  A, 

ir  -  o.  <ot  «■  V  Aj  -  «»tA2( 

Continuity  vi  “  wr  "  tft‘ 

tzom  aquation  (2.20) ,  aquation  (2)  can  ba  writtan  aa, 

°i  °r  °t 

P1C1  '  P1C1  *  B2C2  ‘ 


(1) 

(2) 

(3) 


Solving  aquation  (3)  for  a ,  and  aubatituting  tha  raault  into  aquation  (1), 
a  solution  for  of  in  tarns  of  can  ba  axpraasad  as  follows. 


o 

r 


A2°2C2  *  A1P1C1 

A2P2C2  *  VlCl 


o  . 
1 


ly  following  tha  sans  procadura  and  solving  aquation  (3)  for  ag,  and 
substituting  tha  rasult  into  aquation  (1) ,  a  solution  for  in  tarns  of 
can  ba  asprassad  as. 


(4) 
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t 


iV2C2 


A2P2C2  *  AlPlCi 


(3) 


Figure  2.6 


Derivation  of  reflection  and 
transmission  coefficients. 
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constant  but  will  be  dependent  on  the  velocity  at  which 
the  wave  propagates  in  the  specimen. 

The  net  particle  velocity  of  interface  1  will  be; 

Vj  ■  vj  -  (  -  vr)  ,  (2.22) 


or; 


V;  -  g;  -  (  -Or  )  ,  (2.25) 

pco 

and  the  particle  velocity  of  interface  2; 

v2  "  <*t  -  vt  .  (2.24) 

P  Co 

Representative  particle  velocities  for  interfaces  1  and  ? 
with  soil  as  the  specimen  are  shown  in  figure  2.8.  By 
taking  the  difference  of  the  particle  velocities  at  each 
interface,  the  rate  at  which  the  specimen  is  straining  can 
be  computed  as; 


C  *  (  V;  -  (  -^Vr  )  -  v»  )  ,  (2.25) 


or; 


B  -  (  Q\  -  (  -  Or  )  -  Ot  )  . 
P  Col 


(2.26) 


The  strain  experienced  by  the  specimen  at  any  time  t,  can 
be  computed  by  taking  the  integral  of  the  strain-rate; 


Figure  2.8.  Particle  velocity  history  at  the  specimen-bar 
interfaces,  experiment  115. 


Figures  2.9  and  2.10  show  typical  results  fc  r  e  and  €• 
respectively,  as  computed  from  equations  (2.26)  and 
(2.27).  If  the  specimen  is  straining  uniformly  then  the 
values  obtained  from  equations  (2.26)  and  (2.27)  will  be 
exact.  However,  the  specimen  will  not  usually  experience 
uniform  straining  over  its  length.  Therefore,  the  computed 
strain-rate  and  strain  can  only  be  considered  average 
values  for  the  specimen. 

The  stresses  at  interfaces  1  and  2  can  be  determined 
in  a  similar  fashion.  The  force  at  interfaces  1  and  2  will 
be; 

Ft  -  (  Oj  +  (  -  Or  ))  Aj  -  o'x  A2  ,  (2.28) 

F2  -  (  at  +  (  -  Ox  )  )  A2  -  Ot  Aj  ,  (2.29) 

where  Ot  is  the  stress  transmitted  into  the  specimen.  Or 
is  the  stress  reflected  at  interface  2,  Ai  is  the  area  of 
the  pressure  bars,  and  A2  is  the  area  of  the  specimen.  The 
stress  at  each  interface  is  then; 

01  -  (  Oj  +  (  -  Or ) )  At  ,  (2.30) 

*2 

02  -  Ot Ai  ,  (2.31) 

A2 

and  the  average  stress  in  the  specimen  will  be; 
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(  O:  +  (  -  Or  )  +  Ot  >  A1 


(2.33) 


Figure  2.11  shows  a  graphic  representation  of  the 
averaging  procedure  and  its  result. 


2.5  Contributions  of  Recent  Investigators 


Since  the  modification  made  by  Kolsky  to  the  original 
Hopkinson  bar  experimental  configuration,  the  only 
significant  change  has  been  the  introduction  of  strain 
gauges  for  data  recording.  The  use  of  strain  gauges  was 
first  reported  by  Hauser,  Simmons,  and  Dorn  (1961). 
Replacement  of  the  condenser  microphones  with  strain 
gauges  avoids  the  requirement  of  differentiating  the 
displacement  records  as  a  continuous  strain  history  can 
now  be  recorded. 

To  simplify  equations  (2.26),  (2.27),  and  (2.33)  some 
investigators  have  assumed  the  stress  along  the  axis  of 
the  specimen  to  be  uniform  (e.g.,  Lindholm,  1964).  From 
this  assumption  it  follows  that; 


al  “  a2 


Ot  -  Oj  +  (  -  Or  )  , 


(2.34) 


(2.35) 


Substitution  of  equation  (2.35),  into  equation  (2.33) 


Figure  2.11.  Interface  stress  averaging. 
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yields ; 


°avg"  *i  •  (2.36) 

Similarly,  substitution  of  equation  (2.35),  into  equation 

(2.26) ,  and  (2.27)  yields; 

C  ■  -  2  (  -  gr )  ,  (2.37) 

IpCo 

e  -  -  2  f(  -  CTr)  dt  .  (2.38) 

lpCo  J0 

Because  equations  (2.36),  (2.37),  and  (2.38)  require  an 
additional  assumption  beyond  those  necessary  for  the 
experiment,  they  were  not  used  in  this  investigation. 
Instead,  all  computations  of  specimen  stress,  strain,  and 
strain-rate  have  been  performed  using  equations  (2.26), 

(2.27) ,  and  (2.33).  A  discussion  of  the  experimental 
assumptions  is  presented  in  section  4.1. 

A  critical  analysis  was  made  of  the  Kolsky  technique, 
and  its  assumptions  by  Davies  and  Hunter  (1963).  Their 
experimental  technique  was  essentially  the  same  as  that 
used  by  Kolsky.  However,  by  using  an  analytic  method,  a 
criterion  was  developed  to  estimate  if  the  assumptions  of 
the  experiment  are  satisfied,  principally  the  uniform 
stress  condition  (see  section  4.1). 

Davies  and  Hunter  (1963)  determined  that  to  minimize 
the  effects  of  friction,  the  specimen  aspect  ratio  should 
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be  approximately  unity.  The  criterion  to  estimate  when 

stress  equilibrium  in  the  specimen  has  been  achieved  is 

based  on  the  duration  of  the  applied  stress  wave  (T)  and 

the  characteristic  velocity  at  which  the  wave  propagates 

2 

through  the  specimen  (Cs  -  Es/Ps  ,  where  Es  is  the  slope 
of  the  tangent  to  the  stress-strain  curve  at  a  point 
(do  /de  )  and  is  the  density  of  the  specimen).  The 
characteristic  velocity  (Cs)  is  that  given  by  the 
I’aylor-von  Karmen  theory  for  a  plastically  deforming 
specimen  (von  Karmen,  1942;  Taylor,  1958).  Combining 
energy  principles  with  the  Taylor-von  Karmen  theory,  the 
derived  criterion  for  estimating  when  stress  equilibrium 
in  the  specimen  has  been  achieved  is  given  as  (Davies  and 
Hunter,  1963); 

do  >  TT2 0<  l2  .  (2.39) 

de  t2 

Davies  and  Hunter  (1963)  recommended  that  when  this 
inequality  is  violated,  equilibrium  has  not  been  reached 
and  the  stress-strain  data  may  be  in  error. 

The  SHPB  experimental  environment  was  studied 
numerically  by  Bertholf  and  Karnes  (1974).  They  used  a 
two-dimensional,  finite-difference,  elastic-plastic,  wave 
propagation  computer  code  to  investigate  the  response  of 
an  aluminum  specimen.  The  aluminum  specimen  was  modeled  as 
a  strain-rate  independent  material.  The  code,  TOODY,  was 
developed  by  Bertholf  and  Benzley  (1968).  The  main 
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objectives  of  their  research  were;  (1)  to  examine  the 
validity  of  the  SHPB  experimental  assignations  (i.e.,  in 
particular,  inertia  and  friction),  (2)  to  establish  bounds 
for  the  experimental  parameters,  (3)  to  evaluate  the 
corrections  proposed  by  previous  investigators,  and  (4)  to 
determine  how  the  material  response  is  influenced  when  the 
experimental  assumptions  are  violated  (see  section  4.1) 
(Bertholf  and  Karnes,  1974). 

Bertholf  and  Karnes  (1974)  determined  that  the  SHPB 
experiment  could  be  used  to  accurately,  and  reliably 
determine  material  response  at  high  rates  of  loading, 
provided  that  certain  precautions  were  taken  to  minimize 
the  effects  of  inertia  and  friction.  It  was  shown  that  the 
effects  of  inertia  and  friction  between  the  specimen  and 
the  bars  could  lead  to  the  conclusion  that  the  material 
response  was  rate  dependent  (i.e.,  inertia  and  friction 
affected  the  degree  to  which  stress  and  strain  uniformity 
in  a  specimen  could  be  achieved),  especially  for  the  case 
of  relatively  thin  specimens  (i.e.,  small  aspect  ratios). 
Their  numerical  computations  showed  that  the  influence  of 
inertia  on  experimental  data  could  be  minimized  by 
bounding  the  maximum  specimen  strain-rate  such  that 
(Bertholf  and  Karnes,  1974); 

D *5x10  cm  s  ,  (2.40) 


where  the  loading  wave  is  a  ramp  and; 
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Tr  >  16  ps  cm1  ,  (2.41) 

D 

where  Tr  is  the  rise-tine  of  the  loading  wave  and  D  is  the 
diameter  of  the  pressure  bars.  In  addition  they  showed 
that  if  sufficient  care  is  taken  to  lubricate  the  ends  of 
the  specimen-bar  interfaces/  the  influence  of  friction  on 
experimental  results  can  be  minimized.  Their  investigation 
also  determined  that  the  criterion  proposed  by  Davies  and 
Hunter  to  determine  if  the  reduced  data  may  be  in  error, 
was  indeed  reasonable. 

Since  Davies'  pioneering  work  several  investigators 
have  throughly  examined  the  dispersive  nature  of  wave 
propagation  in  elastic  bars  (e.g.,  Curtis,  1960;  Yeung  Wye 
Kong,  Parsons,  and  Cole,  1974).  However,  only  recently  has 
a  numerical  procedure  to  account  for  wave  dispersion  been 
developed  that  can  easily  be  incorporated  into  the 
standard  SHPB  data  reduction  technique  (Follansbee  and 
Frantz,  1983;  Follansbee  and  Frantz,  1984).  This  technique 
has  been  adapted  for  use  in  this  research. 

As  mentioned  in  section  2.2,  the  generated  stress 
wave  is  not  composed  of  a  single  frequency  but,  instead,  a 
spectrum  of  frequencies,  with  each  frequency  travelling  at 
its  own  respective  phase  velocity  (Cn)  and  wavelength 
(A).  Due  to  this  frequency  variation,  oscillations 
develop  in  the  propagating  stress  wave  which  have  no 
counterpart  in  the  applied  stress  wave  (Curtis,  1960).  The 


nature  of  the  wave  oscillations  can  be  described 
mathematically  through  the  use  of  the  dispersion  equation 
(see  section  2.2).  The  dispersion  equation  can  be  written 
in  the  form  (Wasley,  1973); 

2g(  h2  +  k2  )Jj  (gR)Jj  (hR)  -  (  h2  -  k2  )2  J0  (gR)  J|  (hR)  (2.42) 

R 

-  4  k2ghJ0 (hR)Jx (gR)»  0  , 

where 

92  -  P  Pi  -  k2  , 

A  *2M 

h2-  pp2  -  k2  , 

M 

J0  is  a  Bessel  function  of  the  first  kind  of  order  zero,  Jj 
is  a  Bessel  function  of  the  first  kind  of  order  one,  R  is 
the  radius  of  the  bar,  k  is  the  wave  number  (27T/A  ),  A 
and  |j  are  Lame's  constants,  p  is  the  circular  frequency, 
and  p  is  the  bar  density.  The  solution  to  the  first  mode 
of  vibration  of  the  dispersion  equation  is  shown  in  figure 
2.12  for  a  material  with  a  Poisson's  ratio  of  0.29. 

The  plot  of  the  fundamental  mode  of  the  dispersion 
equation  shows  that  a  high  frequency  wave  will  travel 
slower  than  a  low  frequency  wave  (see  figure  2.12).  Hence, 
as  a  wave  propagates,  the  higher  frequency  components  of 
the  wave  will  lag  behind  the  lower  frequency  components 
causing  a  change  in  the  original  shape  of  the  wave.  This 
wave  dispersion  as  related  to  the  SHPB  experiment 


on  to  the  fundamental  mode  of  vibration  of 
spersion  equation  (adapted  from  Davies,  1948) 
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increases  the  difficulty  in  interpreting  specimen  behavior 
as  deduced  from  the  standard  data  reduction  technique 
(Follansbee  and  Frantz,  1983). 

An  assumption  of  the  dispersion  correction  procedure 
is  that  only  the  fundamental  mode  of  vibration  is  excited 
during  the  SHPB  experiment.  This  assumption  is  supported 
by  the  work  of  Davies  (1948)  and  Follansbee  and  Frantz 
(1983).  Based  on  this  assumption,  the  fundamental  mode  of 
vibration  of  the  dispersion  equation  has  been  approximated 
by  fitting  the  data  with  an  equation  of  the  form 
(Follansbee  and  Frantz,  1983); 

SR  *  0.5764  + _ Q  .  4236 _  (2.43) 
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Figure  2.13  compares  the  approximation  to  the  solution  of 
the  dispersion  equation  for  the  fundamental  mode  of 
vibration.  Also,  the  applied  stress  wave  has  been 
represented  as  a  Fourier  cosine  series; 

oo 

f  (t)  ■  Ao  +  I  Dn  cos (n W0t  -  )  ,  (2.44) 

2  n=1 

where  A0  is  the  amplitude  of  the  largest  frequency 
component,  <jUq  is  the  lowest  frequency  component,  Dn  is  the 
amplitude  of  frequency  component  na>o  and  <J>  is  the  phase 
angle.  The  phase  angle  is  the  component  used  to  correct 
the  waveform  for  dispersion.  A  complete  mathematical 
description  of  the  dispersion  correction  procedure  is 


satisfied  (see  section  4.1>.  The  equations  they  present  to 
determine  the  specimen  stress  aud  strain  assume  that  the 
stress  is  uniform  across  the  specimen.  This  assumption 
needs  to  be  questioned  for  soil  specimens. 

A  more  recent  investigation  into  the  response  of 
soils  using  the  SHPB  has  been  conducted  by  Gaffney,  Brown, 
and  Felice  (1985).  The  soil  was  a  clayey  sand  obtained 
from  the  CARES-Dry  site,  located  on  Luke  Air  Force  Base, 
Arizona.  The  specimens  were  prepared  by  static  compaction 
in  thick-walled  confining  cylinders.  Four  specimens  were 
prepared  at  a  moisture  content  of  3.5  percent  and  a  wet 
density  of  1.77  g/cm3  .  Two  nominal  specimen  lengths  were 
used,  13  mm,  and  25  mm.  The  initial  volume  of  air 
contained  in  the  specimen  was  approximately  29  percent  of 
the  total  specimen  volume.  The  applied  stresses  for  the 
experiments  ranged  up  to  500  MPa. 

In  addition  to  the  SHPB  experiments,  gas  gun 
experiments  at  high  strain-rates  (  >  5000  s1)  and 
quasi-static  experiments  (strain-rates  on  the  order  of 
(5  X  IQ3  s'1)  were  also  conducted.  Through  comparison  of  the 
results  obtained  over  a  range  of  strain-rate  regimes,  it 
was  concluded  that  for  strain-rates  below  5000  s1  the 
response  of  the  soil  with  a  moisture  content  of  3.5 
percent  was  independent  of  strain-rate.  It  was  noted  that 
in  the  SHPB  experiments,  the  specimens  retained  some 
volume  of  air  and  that  the  soil  response  may  be 
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significantly  different  if  all  the  air  voids  had  been 
closed  (Gaffney  et  al.,  1985). 

The  same  SHPB  apparatus  used  by  Gaffney  et  al.  (1985) 
was  also  used  for  this  investigation. 


3.  EXPERIMENTAL  APPARATUS  AND  MATERIAL 


3.1  Experimental  Apparatus 

A  photograph  of  the  SHPB  apparatus  used  in  this 
research  is  shown  in  figure  3.1.  Pigure  3.2  presents  a 
diagrammatic  representation  of  the  apparatus.  Pigure  3.3 
presents  an  enlarged  schematic  of  the  incident  and 
transmitter  pressure  bars  showing  the  specimen  positioning 
and  the  location  of  the  strain  gauges.  The  apparatus  is 
the  property  of  the  Geophysics  Group  at  the  Los  Alamos 
National  Laboratory,  Los  Alamos,  New  Mexico.  The  main 
components  of  the  system  are  the  gas  gun,  the  reaction 
frame,  and  the  incident  and  transmitter  pressure  bars 
(additional  components  are  noted  in  figure  3.2). 

The  incident  and  transmitter  pressure  bars  are 
constructed  of  Vascomax  3S0  CVM  maraging  steel  that  has 
been  heat  treated  to  sustain  a  yield  stress  of  about  2 
GPa.  Each  pressure  bar  is  60.3  mm  in  diameter  and  1.22  m 
in  length.  The  bars  ride  in  adjustable  teflon  bearings 
that  allow  unrestricted  motion  in  the  horizontal  plane. 

The  teflon  bearings  do  not  restrict  the  passage  of  the 
stress  wave  (Lindholm,  1964).  Adjustments  of  the  teflon 
bearings  for  system  alignment  are  made  through  four-wedge 
supported  mounting  rings.  The  entire  system  is  supported 
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by  a  massive  four-column  reaction  frame  mounted  on  a  steel 
I-beam. 

The  applied  stress  wave  is  initiated  by  the  impact  of 
the  striker  bar  on  the  incident  bar.  The  striker  bar  is 
constructed  of  the  same  material  and  has  a  slightly  larger 
diameter  (60.5  mm)  than  the  pressure  bars.  Three  striker 
bar  lengths  were  available/  0.127,  0.254,  and  0.508  ro.  To 
minimize  misalignment  during  impact  and  increase  the  rise 
time  of  the  applied  stress  wave,  the  impact  end  of  the 
striker  bar  is  slightly  rounded  (Follansbee  and  Frantz, 
1983).  The  striker  bar  is  accelerated  down  the  launch  tube 
to  the  appropriate  impact  velocity  by  the  gas  gun.  The 
amplitude  of  the  applied  stress  wave  is  proportional  to 
the  impact  velocity  of  the  striker  bar  which  is  controlled 
by  the  gas  gun  breech  pressure  (Nagy  and  Muelenhaupt, 
1983).  A  description  of  the  gas  gun  is  presented  in 
appendix  B  along  with  the  correspondence  between  breech 
pressure  and  the  striker  bar  impact  velocity.  In  appendix 
C  the  relationship  between  the  impact  velocity  of  the 
striker  bar  and  the  amplitude  of  the  applied  stress  wave 
is  derived.  The  end  of  the  launch  tube  is  vented  so  that 
the  gas  driving  the  striker  bar  is  expelled  before  impact. 
This  is  done  to  avoid  multiple  impacts  by  the  striker  bar 
on  the  incident  bar.  Just  prior  to  impact,  the  striker  bar 
velocity  is  measured  by  three  sets  of  diode  lasers  and 
photodetectors  mounted  in  the  end  of  the  launch  tube.  The 
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laser  velocity  measuring  system  is  described  in  appendix 
B . 

The  applied  stress  wave  in  the  bars  is  monitored  by 
resistance  strain  gauges  mounted  on  the  radial  surface  of 
the  pressure  bars  (see  figure  3.3).  The  strain  gauges 
(Micro-Measurement  No.  EA-06-125AV-12Q  with  extra  long 
foil  tabs)  are  mounted  in  pairs  on  opposite  sides  of  the 
respective  pressure  bars  and  connected  in  a  half-bridge 
configuration  to  nullify  bending  strains.  The  strain 
gauges  are  bonded  (bonding  agent  is  CA120SL  Bean  adhesive) 
to  the  pressure  bars  at  mid-length.  The  positioning  of  the 
strain  gauges  on  the  incident  bar  is  particularly  critical 
because  it  is  necessary  to  obtain  continuous  records  of 
the  incident  and  reflected  waves  without  overlap. 

The  data  recorded  from  strain  gauge  bridges  are 
filtered  and  preamplified  and  then  routed  to  a  data 
acquisition  system.  The  data  are  recorded  by  CAMAC 
(IEEE-583)  based  waveform  digitizers.  The  incident  and 
transmitter  bar  records  were  recorded  on  separate  channels 
with  8-bit  resolution  at  a  rate  of  0.5  microseconds  per 
data  point  for  all  but  experiments  31  and  32  where  the 
data  collection  rate  was  0.25  microseconds  per  data  point. 
The  duration  of  data  recording  is  8  milliseconds  per 
channel.  Positioning  of  this  data  collection  window  such 
that  only  the  data  of  interest  is  gathered  is  controlled 
by  a  timing  signal  generated  by  the  laser  velocity 
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measuring  system  (see  appendix  B) .  The  data  are  read  by  a 
Digital  Equipment  Corporation  LSI-11/02  microcomputer  and 
stored  on  a  flexible  disk.  The  data  are  then  transferred 
by  magnetic  tape  to  a  CRAY-1S  mainframe  computer  located 
at  the  Air  Force  Weapons  Laboratory,  Kirtland  Air  Force 
Base,  New  Mexico  for  complete  data  processing.  Additional 
information  on  the  SHPB  instrumentation  is  presented  in  in 
appendix  B. 

An  additional  feature  of  this  particular  SHPB  is  that 
a  static  preload  of  up  to  200  MPa  can  be  applied  to  the 
specimen.  This  is  accomplished  by  a  hydraulic  ram  located 
at  the  downstream  end  of  the  transmitter  bar  (see  figure 
3.2).  A  tapered  transition  plug  is  used  to  transfer  the 
load  from  the  hydraulic  ram  system.  Another  transition 
plug  is  located  at  the  impact  end  of  the  incident  bar 
which  is  fabricated  with  a  large  diameter  shoulder  to 
support  the  static  preload  that  is  applied  to  the  specimen 
(see  figure  3.2).  It  is  tapered  to  match  the  diameter  of 
the  incident  bar  where  they  contact.  The  transition  plugs 
are  constructed  of  the  same  material  as  the  pressure  bars 
and  the  striker  bar.  A  static  preload  was  not  applied  to 
the  specimens  in  this  research. 

Appendix  D  includes  a  detailed  outline  of  the 
procedure  followed  when  conducting  a  SHPB  experiment. 
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3.2  Soil  Characteristics 

The  soil  used  for  the  experiments  described  herein 
was  sampled  in  bulk  quantities  from  a  test  pit  in  an 
undisturbed  area  of  the  McCormick  Ranch  test  site  located 
on  Kirtland  Air  Force  Base,  New  Mexico.  The  site  has  been 
used  for  explosive  testing  by  the  Department  of  Defense 
since  the  early  1960s.  The  geology  of  the  site  consists  of 
some  surface  playa  underlain  by  unconsolidated  eolian 
sand,  alluvium,  and  lacustrine  deposits  15  to  30  m  thick 
(Bedsun,  1983).  The  water  table  is  estimated  to  lie  at  a 
depth  of  60  to  90  m  (Bedsun,  1983). 

In  order  that  the  soil  be  as  free  of  organic  material 
as  possible,  the  surface  vegetation  was  scraped  away  and 
the  sample  taken  from  a  depth  of  1  to  2  m.  After  arriving 
at  the  University  of  Utah  soils  laboratory,  the  samples 
were  mixed  to  oreak  apart  large  clumps  and  to  achieved  a 
uniform  r-.xxture.  The  following  tests  were  then  performed 
to  determine  the  soil  properties:  (1)  mechanical  sieve  and 
hydrometer  analyses  to  determine  grain-size  distribution, 
(2)  Atterberg  limits,  (3)  hygroscopic  moisture 
determination,  and  (4)  specific  gravity.  These  tests  were 
conducted  in  accordance  with  the  American  Society  for 
Testing  and  Materials  standards  (ASTM,  1982).  The 
moisture/densi ty  relationship  was  determined  by  the 
Harvard  miniature  compaction  procedure  using  three  layers 
with  25,  40  pound  tamps  per  layer.  An  x-ray  diffraction 
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study  was  also  performed  to  determine  the  soil 
constitutive  minerals.  This  study  conducted  at  the  Los 
Alamos  National  Laboratory,  Los  Alamos,  New  Mexico. 

The  results  of  these  tests  showed  the  soil  to  be  a 
nonplastic,  clayey  silty  sand  (see  figure  3.4)  with  a 
specific  gravity  of  2.67  and  hygroscopic  moisture  content 
of  2.9  percent.  The  constitutive  minerals  and  their 
percentages  are,  quartz,  60  percent,  feldspar,  30  percent, 
with  traces  of  hornblende,  illite,  montmoril Ionite, 
cristobalite,  and  calcite.  A  carbonate  test  showed  the 
percentage  of  calcite  to  be  approximately  2  percent. 

3 ♦ 3  Specimen  Preparation 

The  majority  of  specimens  were  prepared  near  the 
optimum  moisture  content  (13.3  %)  and  dry  density  (1.87 
g/cm3 )  as  determined  by  the  Harvard  miniature  compaction 
procedure  (see  figure  3.5).  A  complete  tabulation  of  the 
initial  specimen  parameters  for  each  specimen  is  presented 
in  table  E.l  of  appendix  E.  To  achieve  as  uniform 
conditions  among  specimens  as  possible,  the  soil  was  mixed 
in  batches  sufficient  to  prepare  a  minimum  of  five 
experimental  specimens.  Before  the  addition  of  water,  the 
soil  was  passed  through  a  No.  4  sieve  (4.75  mm  opening), 
weighed  and  placed  in  a  large  flat  pan.  The  correct  amount 
of  moisture  was  added  by  using  a  spray  bottle  so  that  an 
even  distribution  could  be  obtained.  The  soil  was  then 
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mixed  throughly.  After  mixing,  a  damp  cover  was  placed 
over  the  soil  for  a  period  of  20  minutes  to  allow  the 
soil-water  mixture  to  equilibrate.  Poliowing  the 
equilibration  period,  the  soil  was  again  mixed  to  assure 
that  an  even  mixture  was  obtained. 

The  individual  experimental  specimens  were  then 
prepared  by  removing  the  appropriate  soil  mass  from  the 
batch  mix  to  yield  a  given  volume  when  compacted.  At  this 
time  a  small  specimen  of  the  batch  mix  was  taken  for  a 
moisture  content  determination.  The  standard  laboratory 
techniqv  r  moisture  content  determination  was  followed 
<  ASTM-i/  -.*16-80).  The  specimens  were  statically  compacted 
in  the  confining  cylinders  by  using  a  hydraulic  press  (see 
figure  3.6)  with  spacer  rings  (see  figure  3.7)  to  control 
the  specimen  length  (i.e.,  density).  A  description  of  the 
specimen  confining  system  is  presented  section  3.4.  Two 
nominal  specimen  lengths  were  used  12.7  mm  and  6.35  mm. 
These  lengths  were  chosen  based  on  the  work  of  Gaffney  and 
Brown  (1984).  They  reported  that  for  specimen  lengths 
greater  than  13  mm  difficulty  was  encountered  in  obtaining 
a  transmitted  stress  signal. 

Bach  specimen  was  then  sealed  in  a  plastic  bag  to 
minimize  any  moisture  loss  that  might  occur  prior  to  the 
experiment.  Eight  tests  were  conducted  to  quantify  the 
amount  of  moisture  loss  to  be  expected  between  the  time 
the  specimen  was  prepared  and  the  time  the  experiment  was 
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Figure  3.7.  Spacer  rings,  confining  cylinder,  and  plungers. 
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performed.  The  results  showed  an  average  loss  was  5.4 
percent  of  the  total  moisture  with  a  standard  deviation  of 
1.8.  These  results  are  believed  to  be  representative  of 
the  pre*experiment  moisture  losses  for  the  specimens 
prepared  near  optimum  compaction  conditions.  The  result 
for  each  test  is  presented  in  table  E.2  of  appendix  E. 

The  same  preparation  procedure  was  followed  for 
specimens  prepared  on  the  wet  and  dry  side  of  optimum 
conditions . 

Prior  to  performing  the  experiment ,  the  specimen  was 
removed  from  the  plastic  bag  and  the  length  of  the 
specimen  was  measured  inside  the  confining  cylinder.  This 
was  accomplished  by  lowering  a  depth  gauge  from  the  lip  of 
the  cylinder  to  the  surface  of  the  specimen  on  both  ends 
of  the  cylinder.  Hence,  by  knowing  the  depth  to  the 
specimen  surface  on  each  end  of  the  cylinder,  and  the 
length  of  the  cylinder,  the  specimen  length  could  be 
computed.  At  the  completion  of  each  experiment,  the 
specimen  length  was  again  measured.  The  specimen  was  then 
extruded  from  the  confining  cylinder  and  a  portion  of  the 
specimen  used  to  determine  the  postexperiment  moisture 
content . 

Whenever  a  specimen  of  soil  is  compacted  there  will 
exist  moisture  and  density  variations  over  its  length.  A 
limited  experimental  study  was  conducted  to  evaluate  the 
extent  of  these  variations  for  specimens  compacted  using 
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the  static  compaction  technique.  The  nominal  length  of  the 
specimens  prepared  for  this  study  was  12  cm.  The 
experimental  procedure  and  results  are  presented  in 
appendix  F.  The  results  indicate  that  in  general,  a 
specimen  will  have  the  highest  density  at  its  ends  and  the 
lowest  density  at  its  center.  The  moisture  content  will 
vary  from  a  high  at  the  center  to  lower  at  the  ends. 

Although  the  specimen  lengths  used  in  the  moisture 
and  density  variation  study  greatly  exceeded  the  specimen 
lengths  used  for  the  SHPB  experiments,  the  results  can 
readily  be  interpolated.  In  short  sections  near  the  center 
of  the  specimen,  the  variations  in  density  and  moisture 
were  less  than  2  and  4  percent,  respectively.  During 
compaction,  these  variations  are  affected  by  friction 
between  the  soil  and  the  confining  cylinder.  For  long 
specimens  a  greater  area  will  exist  over  which  friction 
forces  can  act.  Hence,  friction  will  influence 
moisture/density  variations  to  a  greater  degree  in  longer 
specimens  chan  it  will  in  short.  Therefore,  these  values 
can  be  considered  upper  bounds  for  the  short  specimens 
used  in  the  SHPB  experiments. 

3.4  Specimen  Confining  System 

To  achieve  a  nearly  uniaxial  strain  environment  for 
the  experiments,  the  soil  specimen  was  compacted  in  a 
thick-walled  confining  cylinder.  The  condition  of  uniaxial 
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strain  will  be  confirmed  in  section  5.3.  The  concept  of 
confining  a  SHPB  specimen  such  that  it  will  experience  a 
condition  of  nearly  uniaxial  strain  was  first  presented  by 
Bhushan  and  Jashman  (1978). 

The  nominal  dimensions  of  the  confining  cylinder  are 
60.3  mm  inside  diameter/  102  mm  outside  diameter/  and  44.5 
mm  in  length.  The  confining  cylinder  served  several 
purposes;  first,  to  contain  the  soil  specimen  itself,  and 
second,  because  the  specimen  will  experience  a  state  of 
nearly  uniaxial  strain,  the  effects  of  radial  inertia  are 
avoided.  The  confinement  also  prevented  specimen 
distortion  or  barrelling  during  the  experiment. 

When  the  confining  cylinder  containing  the  specimen 
was  placed  between  the  pressure  bars  approximately  19.0  mm 
of  the  cylinder  overlapped  the  bars  on  each  end  (see 
figure  3.8).  To  determine  if  the  confining  cylinder  was 
transferring  any  stress  to  the  transmitter  bar,  a  test  was 
conducted  with  the  bars  separated  a  distance  of  3.0  mm  and 
the  confining  cylinder  placed  over  the  air  gap.  The  3.0  mm 
distance  was  greater  than  the  anticipated  displacement  of 
the  incident  bar.  If  the  confining  cylinder  did  indeed 
transfer  stress  to  the  transmitter  bar,  a  signal  would  be 
recorded  at  the  transmitter  bar  strain  gauge;  if  not,  the 
strain  gauge  record  should  be  flat.  The  results  indicated 
that  the  confining  cylinder  did  not  transfer  any 
measurable  stress  to  the  transmitter  bar. 
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stress  wave,  many  reflections  can  take  place  within  the 
specimen.  The  stress  gradient  will  diminish,  resulting  in  a 
nearly  uniform  distribution  of  stress  over  the  length  of 
the  specimen.  If  there  is  insufficient  time  for  the  stress 
in  the  specimen  to  equilibrate,  the  experiment  is  one  of 
wave  propagation  rather  than  one  of  uniform  stress.  Should 
this  be  the  case,  a  constitutive  equation  would  be  required 
to  study  the  phenomena.  This  defeats  the  purpose  of  the 
experiment  as  it  is  the  the  constitutive  equation  that  is 
sought  in  the  first  place.  Hence,  if  care  is  exercised  in 
selecting  specimen  geometry  and  the  duration  of  the  applied 
stress  wave,  the  SHPB  experiment  provides  a  method  of 
obtaining  a  material  relationship  without  the  disadvantages 
of  having  to  undertake  a  wave  propagation  analysis. 

This  assumption  is  perhaps  the  most  difficult  to 
satisfy  when  using  soil  specimens  in  a  SHPB  experiment. 
Because  soils  exhibit  nonlinear  hysteretic  behavior  and 
have  low  wave  speeds  (e.g.,  =  300  m/s)  there  may  be 
insufficient  time  for  the  stress  gradient  between  the  two 
faces  of  the  specimen  to  diminish.  Therefore,  to  satisfy 
the  uniform  stress  assumption  the  specimen  aspect  ratio 
will  have  to  be  small  and  the  duration  of  the  applied 
stress  wave  long. 
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4.1.2  Wave  Dispersion 


The  stress  wave  initiated  by  the  impact  of  the 
striker  bar  on  the  incident  pressure  bar  is  initially  very 
complex.  Thij  is  due  in  part  to  friction  at  the  striker 
bar-incident  bar  interface ,  as  well  as  the  appearance  of 
spherical  and  dilatational  waves  (Follansbee  and  Frantz, 
1983).  These  end  effects  disappear  at  about  10  bar 
diameters  from  the  point  of  impact  (Wasley,  1973,  and  Yeung 
Wye  Kong  et  al.,  1974). 

Figure  4.1  shows  the  waves  as  recorded  by  the  incident 
bar  strain  gauge  in  a  test  in  which  no  specimen  was  placed 
between  the  incident  and  transmitter  bars.  The  incident  and 
transmitter  bars  were  separated  such  that  the 
specimen-incident  bar  interface  was  a  free  end.  In  this 
configuration,  the  applied  stress  wave  upon  reaching  the 
free  end  is  reflected  as  a  tensile  wave  with  amplitude 
equal,  but  opposite  in  sign  to  the  applied  wave.  If  the 
wave  does  not  disperse  during  propagation,  the  two  waves 
should  be  identical.  Clearly,  wave  dispersion  does  occur  as 
the  period  of  the  oscillations  of  the  reflected  wave  are 
greater  than  those  of  the  incident  wave. 

Figures  4.2  and  4.3  demonstrate  the  application  of  the 
dispersion  correction  (see  appendix  A)  to  the  waves 
generated  in  the  SHPB  test  with  the  specimen- incident  bar 
interface  as  a  free  end  (see  section  2.5).  Figure  4.2 
compares  the  reflected  stress  wave  as  recorded  at  the 
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Figure  4.3.  Corrected  incident  and  reflected  (inverted)  stress 
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strain  gauge  location  with  the  dispersion  corrected  wave  a? 
reconstructed  at  the  specimen- incident  bar  interface.  It  is 
evident  that  the  dispersion  correction  has  effectively 
compressed  the  reflected  wave  during  the  reconstruction. 
This  is  to  be  expected  as  the  phase  angle  will  decrease  for 
negative  values  of  Ax  (see  appendix  A,  equation  A. 26).  In 
figure  4.3  the  dispersi.  .orrected  incident  and  reflected 
stress  waves  as  reconstructed  at  the  specimen-incident  bar 
interface  are  shown.  From  figure  4.3  it  can  be  seen  there 
is  a  correspondence  of  the  peaks  and  valleys  of  the  two 
waves.  This  correspondence  was  not  evident  between  the 
uncorrected  stress  waves  at  the  recorded  strain  gauge 
position  (see  figure  4.1).  This  correspondence  supports  the 
assumption  that  only  the  fundamental  mode  of  vibration  is 
excited  in  the  SHPB  experiment  (Pollansbee  and  Frantz, 
1983).  If  other  modes  of  vibration  had  been  present,  the 
procedure  of  performing  the  dispersion  correction  based  on 
obtaining  the  waves'  component  phase  velocity  and 
wavelength  from  the  fundamental  mode  would  not  have  been 
successful  (Follansbee  and  Frantz,  1983). 

Figure  4.1  clearly  shows  that  wave  dispersion  does 
occur,  while  figure  4.3  shows  that  the  Follansbee  and 
Frantz  (1983)  correction  procedure  can  be  used  to  account 
for  the  phenomenon.  If  dispersion  is  not  accounted  for  in 
the  reduc*^  experimental  lata,  large  oscillations  will 
appear  ir.  t'?».  stress-stra  -i  curv  whi.cn  can  mask  trends  and 


increase  the  difficulty  in  interpreting  specimen  response 
(see  section  6.1.2). 


4.1.3  Stress  Variation  Over  the  Cross 


Section  of  the  Bar 


The  information  required  to  evaluate  the 
strain-rate-time,  strain-time,  and  stress-time  response  of 
a  specimen  in  a  SHPB  experiment  is  obtained  from  strain 
gauges  mounted  on  the  surface  of  the  pressure  bars  (see 
figure  3.3).  If  this  recorded  information  is  not  consistent 
with  the  stress-time  response  the  pressure  bar  is 
experiencing  over  its  cross  sectional  area,  the  resulting 
computations  may  be  seriously  in  error.  Hence,  a 
significant  condition  that  must  be  present  in  any  SHPB 
experiment  is  that  the  longitudinal  stress  and  displacement 
be  nearly  constant  over  the  cross  sectional  area  of  the 
bar. 

A  prediction  of  the  dispersion  equation  is  that  the 
longitudinal  stress  and  displacement  will  vary  over  the 
cross  sectional  area  of  the  pressure  bars,  and  that  the 
variation  will  be  dependent  on  the  ratio  R  /A  (Davies, 
1948).  Therefore,  errors  are  unavoidable  if  one-dimensional 
wave  analysis  is  used  to  reduce  the  data.  The  significance 
of  the  error  and  its  impact  on  interpreting  experimental 
results  was  investigated  by  Davies  (1948).  Using  the 
limiting  condition  of  R  /A  •  0  as  an  error  baseline  (i.e. , 
theoretically  at  R  /A  *  0  the  longitudinal  displacement 


and  stress  are  constant  over  the  cross  sectional  area  of 
the  bar),  Davies  was  able  to  establish  that  for  R  /A  <  0.1, 
the  displacement  recorded  at  the  surface  of  the  pressure 
bar  would  differ  from  that  at  the  bars  longitudinal  axis  by 
±  3  percent  (Davies,  1948). 

In  this  research  the  diameter  of  the  pressure  bars  are 
larger  than  those  used  by  past  investigators  by  at  least  a 
factor  of  2  and  in  some  cases  a  factor  of  3.  For  example/ 
in  the  experiments  performed  by  Lindholm  (1964)  the 
diameter  of  the  pressure  bars  were  1.9  cm.  A  concern  in 
using  the  larger  bars  is  that  the  behavior  over  the  cross 
section  will  deviate  from  one-dimensional  condition  that  is 
assumed;  hence  increasing  the  error  in  the  computations 
beyond  that  established  by  Davies.  Any  additional  error  and 
its  affect  on  the  computed  specimen  stress-strain  response 
has  been  addressed  in  section  6.1.3. 

4.1.4  Specimen-Bar  Interface  Friction 

In  the  traditional  5HPB  experiment,  a  specimen  with  a 
diameter  slightly  less  than  that  of  the  bars  is  placed 
between  them  (Lindholm,  1964).  This  is  to  allow  the 
specimen  to  expand  radially  during  the  experiment  while  not 
exceeding  the  diameter  of  the  pressure  bars.  As  stress  is 
applied  to  the  specimen,  radial  shear  stresses  are  created 
between  the  pressure  bars  and  the  specimen.  This  has  been 
commonly  referred  to  as  "end  effects"  or  "friction 
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ef facts."  A  result  of  these  effects  is  that  the  specimen 
tends  to  be  clamped  at  the  specimen-bar  interfaces  (Green 
and  Perkins ,  1969).  This  prevents  the  specimen  from 
expanding  uniformly;  hence,  barreling  of  the  specimen  is 
often  observed.  This  may  result  in  a  strengthing  effect 
being  observed  in  the  stress-strain  response  (Betholf  and 
Karnes,  1974).  Zn  addit.o-^,  these  effects  may  greatly 
increase  the  degree  of  nonuniformity  of  stress  and  strain 
in  the  specimen  (Bertholf  and  Karnes,  1974). 

To  reduce  the  effects  of  friction  it  is  common 
practice  to  apply  a  film  of  lubrication  to  the  ends  of  the 
pressure  bars  in  contact  with  the  specimen.  The  results  of 
the  numerical  study  conducted  by  Bertholf  and  Karnes  (1974) 
have  shown  this  to  be  an  acceptable  method  to  miminize  the 
effects  of  friction  on  stress-strain  response.  However,  the 
effects  of  friction  cannot  be  considered  independently  of 
specimen  geometry,  i.e.,  aspect  ratio  (Rand,  1967).  Rand 
(1967)  has  shown  that  decreasing  the  specimen  aspect  ratio 
will  have  the  same  effect  as  increasing  the:  coefficient  of 
friction  between  the  specimen  and  the  pressure  bars.  Again 
as  a  result  of  the  work  of  Bertholf  and  Karnes  (1974),  a 
specimen  aspect  ratio  of  0.5  is  thought  to  be  an  optimum. 

However,  for  soil  specimens  an  aspect  ratio  of  0,5  may 
be  too  restrictive  due  to  the  same  concerns  mentioned  in 
section  4.1.1  (i.e.,  low  wave  speed  and  nonlinear 
hysteretlc  behavior).  Therefore  the  aspect  ratio  for  a  soil 
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specimen  must  be  large  enough  so  that  friction  effects  are 
minimized  but  small  enough  to  permit  stress  equilibration 
within  the  duration  of  the  applied  stress  wave.  A  method  to 
retain  a  small  aspect  ratio  while  avoiding  friction  effects 
is  to  contain  a  soil  specimen  in  a  confining  cylinder.  The 
confining  cylinder  prevents  barreling  and  excessive  radial 
strain  so  that  the  specimen  is  in  a  state  of  nearly 
uniaxial  strain  (see  section  5.4).  Hence,  friction  effects 
on  specimen  stress-strain  response  is  minimized  (see 
section  6.1.4). 


This  section  will  present  the  steps  used  in  the  data 
reduction  procedure  to  compute  the  average  specimen 
strain-rate-time,  strain-time,  and  stress-time  response. 
The  steps  are: 

(1)  data  input, 

(2)  correct  each  wave  for  dispersion, 

(3)  compute  the  average  specimen  strain-rate-time, 

strain-time,  and  stress-time  response  using  equations 

(2.26),  (2.27),  and  (2.33),  and 

(4)  output  processd  data. 

To  perform  the  computations  required  in  the  data 
reduction  procedure  a  computer  program  (REDUCE)  was 
developed.  A  flowchart  and  listing  of  REDUCE  is  presented 
in  appendix  G  along  with  an  example  of  the  output  for 
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experiment  134. 

4.2.1  Raw  Data 

Examples  of  the  raw  data  recorded  by  the  strain  gauges 
for  the  two  nominal  specimen  lengths  are  shown  in  figures 
4.4  and  4.5.  Figures  4.6  and  4.7  compare  the  raw  data  of 
figures  4.4  and  4.5  with  the  data  after  the  dispersion 
correction  procedure  has  been  performed.  Because  the  raw 
data  are  recorded  at  some  distance  away  from  the 
specimen-bar  interfaces,  and  the  equations  used  in  the 
computations  are  only  valid  at  the  specimen-bar  interfaces, 
each  wave  must  be  corrected  for  the  dispersion  that  occurs 
between  the  strain  gauge  and  the  specimen-bar  interface. 

4.2.2  Processed  Data 

Figures  4.8  and  4.9  illustrate  the  waves  of  figures 
4.4  and  4.5  after  correction  for  dispersion  and  time 
shifting  to  the  specimen-bar  interfaces.  The  incident  and 
reflected  waves  are  shifted  to  the  specimen-incident  bar 
interface  so  that  they  will  start  at  the  same  moment.  The 
transmitted  wave  will  be  shifted  to  the 

specimen-transmitter  bar  interface  but  with  a  delay  which 
depends  on  the  wave  velocity  in  the  specimen.  As  the 
distance  between  the  specimen-bar  interface  and  the  strain 
gauge  is  the  same  for  both  pressure  bars,  Ax,  each  wave  is 
time  shifted  by  Ax/C q,  where  Co  (4900  m/s)  is  the  wave 
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Figure  4.4.  Uncorrected  stress  waves  for  a  6.35  mm  specimen 
experiment  115. 
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Figure  4.8.  Stress  waves  for  a  6.35  mm  specimen  corrected  and 
time  shifted  to  the  specimen-bar  interfaces, 
experiment  IT. 5. 


Figure  4.9.  Stress  waves  for  a  12.7  mm  specimen  corrected  and 
time  shifted  to  the  specimen-bar  interfaces, 
experiment  134. 


velocity  in  the  pressure  bars. 

Examples  of  the  average  s train-rate-time,  strain-time 
stress-time,  and  stress-strain  response  as  computed  by 
REDUCE  for  the  two  nominal  specimen  lengths  are  shown  in 
figures  4.10  through  4.17. 


0) 

w 

3 

O' 

■H 

ttt 


7S 


* 

[* 


i 

J. 

,1 

.1 


A 

i 


'i 

4 


Figure  4.11.  Strain-rate-tirae  response  for  a  12.7  mm  specimen,  experiment 


Figure  4.14.  Stress-time  response  for  a  6.35  mm  specimen 
experiment  115. 


Figure  4.15.  Stress-time  response  for  a  12.7  mm  specimen 


115. 


Figure  4.17.  Stress-strain  response  for  a  12.7  nun  specimen, 
experiment  134. 


5.  PRESENTATION  OF  EXPERIMENTAL  RESULTS 


5.1  Experimental  Results 

In  this  research  50  experiments  were  conducted  on 
compacted  sand  specimens.  The  compaction  moisture  and 
density  combinations  ranged  from  dry  of  optimum  to  wet  of 
optimum  conditions  as  determined  by  the  Harvard  miniature 
compaction  procedure  (see  figure  3.5).  The  majority  of 
experiments  were  performed  with  specimens  prepared  at  or 
near  optimum  conditions.  In  most  cases,  a  minimum  of  two 
experiments  were  conducted  at  each  combination  of  specimen 
moisture  content  and  density.  The  moisture/density  data 
points  for  each  group  (i.e.,  specimens  prepared  dry,  near, 
and  wet  of  optimum)  and  their  relationship  to  the  Harvard 
miniature  compaction  curve  are  presented  in  appendix  H. 

The  applied  stresses  and  maximum  strain-rates  for  the 
experiments  ranged  from  130  MPa  to  760  MPa,  and  from  500  s1 
to  5000  s'1,  respectively.  The  applied  stress,  peak  stress, 
strain  at  peak  stress,  and  strain-rate  at  peak  stress  for 
each  experiment  are  tabulated  in  table  E.3  of  appendix  E. 

Three  different  length  striker  bars  were  available  to 
generate  the  applied  stress  wave,  0.127,  0.254,  and 
0.508  m.  Several  experiments  were  performed  with  the 
0.508  m  striker  bar;  however,  the  incident  wave  and  the 
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reflected  wave  were  found  to  overlap.  Hence,  the 
termination  of  the  incident  wave  and  the  commencement  of 
the  reflected  wave  could  not  be  resolved.  Dispersion  of 
the  incident  wave  was  determined  to  be  the  cause  of  the 
overlap,  and  not  the  incorrect  placement  of  the  incident 
bar  strain  gauge.  It  was  shown  in  section  4.1.2  that  for 
positive  values  of  Ax,  the  phase  angle  will  increase 
causing  the  propagating  wave  to  be  stretched.  The  data 
from  these  experiments  were  not  reduced.  No  experiments 
were  attempted  with  the  0.127  m  striker  bar  as  the 
specimen  would  not  be  able  to  reach  equilibrium  within  the 
duration  of  the  stress  wave  (£50  microseconds).  Hence, 
the  0.254  m  striker  was  selected  to  initiate  the  stress 
wave  for  all  the  experiments  reported  herein. 

Although  50  experiments  were  performed,  only  the 
results  of  27  are  used  for  the  purpose  of  comparison.  In 
the  process  of  preparing  the  SHPB  apparatus  for  firing,  a 
seating  strain  was  applied  to  the  specimen.  The  amount  of 
seating  strain  sustained  by  a  specimen  varied  from 
experiment  to  experiment.  This  was  an  artifact  of  the  SHPB 
apparatus  and  required  a  change  in  the  experimental 
procedure  so  that  the  inherent  variability  was  accounted 
for  in  the  data  reduction  procedure.  The  cause  of  this 
difficulty  and  the  procedure  that  was  developed  to  measure 
the  variable  seating  strain  is  discussed  in  section  6.2. 

The  procedure  to  measure  the  seating  strain  was 
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applied  to  experiment  112  and  every  experiment  thereafter 
(8ee  section  6.2).  Accurate  seating  strains  are  necessary 
to  determine  the  change  in  specimen  length  and  phase 
relationships  (e.g.,  initial  gas  porosity).  A  knowledge  of 
these  parameters  is  essential  as  they  will  influence  a 
specimens  response  to  a  given  stress.  If  accurate 
assessments  are  to  made  concerning  for  example 
experimental  replication/  the  initial  parameters  of  the 
specimens  being  compared  must  be  similar.  This  is 
particularly  critical  for  soil  specimens.  Hence/  in  the 
experiments  where  seating  strains  were  not  measured  (i.e., 
31  through  76),  the  results  are  not  suitable  for  the 
purpose  of  comparison. 

The  average  seating  strain  experienced  by  the 
specimens  in  experiment  112  through  167  was  7.5  percent 
with  a  standard  deviation  of  3.8.  The  results  are 
tabulated  in  table  G.4  of  appendix  E.  For  the  experiments 
where  the  seating  strain  was  not  measured,  the  initial 
specimen  length  was  adjusted  by  7.5  percent  for  use  in  the 
computations.  Figure  5.1  shows  the  stress-strain  response 
for  experiment  55  using  the  specimen  length  that  was 
measured  just  before  it  was  placed  between  the  bars  and 
the  adjusted  length.  It  is  evident  that  the  overall 
reponse  is  not  affected  by  adjusting  the  specimen  length 
by  7.5  percent.  However,  the  correct  response  to  use  for 
comparison  is  not  known.  It  will  be  shown  throughout  this 


chapter  the  response  is  not  inconsistent  with  the  response 
obtained  in  experiments  where  the  seating  was  measured. 
Hence,  although  the  specimen  stress-strain  response  for 
experiments  31  through  76  is  consistent,  they  are  not 
suitable  for  the  purpose  of  comparison.  The  computed 
stress-strain  response  and  summary  information  for  each  of 
the  50  experiments  are  presented  in  appendix  I. 

5.2  Stress-Strain  Response 

Figure  5.2  shows  the  stress-strain  response  for  a 
12.7  mm  specimen  to  an  applied  stress  of  400  MPa.  The 
response  is  representative  for  both  specimen  lengths 
subjected  to  the  applied  stresses  used  in  this  research. 

The  stress-strain  curve  can  be  divided  into  three 
regions  which  describe  the  overall  specimen  response.  The 
first,  region  0  to  A,  is  where  inertia  effects  act  to 
oppose  the  equilibration  of  stress.  The  extent  of  this 
region  can  be  estimated  by  using  the  criterion  of  Davies 
and  Hunter  (1963)  as  well  as  experimentally  (see  section 
6.1.1).  At  A  the  specimen  is  considered  to  be  in 
quasi-equilibrium  and  the  deformation  nearly  homogeneous. 
Region  A  to  B  is  characterized  by  particle  rearrangment 
into  a  denser  packing  and  closing  of  the  air  voids.  The 
strain  at  B  is  equal  to  the  initial  gas  porosity  of  the 
specimen.  At  this  point  there  are  theoretically  no  air 
voids  remaining. 


The  initial  gae  porosity  of  the  specimen  is  computed 
as  the  volume  of  air  contained  in  the  specimen  divided  by 
the  total  volume  of  the  specimen.  In  dynamic  uniaxial 
strain  experiments  on  partially  saturated  sand  where  the 
stress  changes  are  large,  the  pore  air  will  be  highly 
compressed,  and  then  dissolved  in  the  pore  water  (Whitman, 
1970).  This  behavior  is  characteristic  of  the  response 
observed  in  region  A  to  B.  With  no  pore  air  remaining 
(zero  gas  porosity),  the  specimen  will  be  fully  saturated 
(point  B).  Under  these  conditions  the  compressive 
resistance  of  the  water  will  greatly  exceed  the 
compressive  resistance  of  the  soil  skeleton  (Whitman, 
1970).  When  this  occurs,  the  specimen  will  become  strongly 
resistant  to  additional  deformation.  This  behavior  governs 
the  response  in  the  region  B  to  C. 

In  the  unloading  region  of  the  curve  (i.e.,  beyond 
point  C) ,  the  specimen  in  many  of  the  experiments 
continued  to  acculumate  strain  (e.g.,  see  figures  4.16  and 
4.17).  The  cause  of  the  additional  strain  accumulation  may 
be  that  the  specimen  can  retain  some  memory  of  its 
loading.  Another  cause  might  be  the  extrusion  of  material. 
The  potential  for  loss  of  material  during  the  experiment 
is  discussed  in  section  6.4.  The  unloading  region  of  the 
curve  will  not  be  used  in  the  analyses. 


5.3  Uniaxial  strain  Path 


A  nearly  uniaxial  strain  state  is  forced  upon  the 
specimen  by  the  fact  that  it  is  contained  in  a  confining 
cylinder.  This  condition  may  be  verified  by  comparing  the 
radial  strain  experienced  by  the  confining  cylinder  and 
maximum  longitudinal  strain  experienced  by  the  specimen. 
The  radial  strain  was  estimated  by  computing  a  maximum 
radial  deflection  of  the  confining  cylinder  using 
thick-walled  cylinder  theory.  The  equation  used  to  compute 
the  radial  deflection  of  a  thick-walled  cylinder  subject 
to  internal  pressure  was  (Seely  and  Smith,  1952); 


AR  -  P  Rt 


R22 


Rl 


+  V 


2  2 
R2  -  Rj 


(5.1) 


where  Ri  is  the  inside  radius  of  the  cylinder,  R2  is  the 
outside  radius  of  the  cylinder,  P  is  the  internal 
pressure,  E  (10  GPa)  is  Young's  modulus  of  the  confining 
cylinder,  and  v  (0.3)  is  Poisson's  ratio  for  the  confining 
cylinder.  The  internal  pressure  was  taken  as  the  peak 
stress  experienced  by  the  specimen.  This  assumes  that 
Poisson's  ratio  for  the  soil  is  0.5.  The  assumption  is 
reasonable  if  the  soil  is  near  saturation  as  the  pore 
water  will  be  highly  resistant  to  volume  change  (Whitman, 
1970)  In  nearly  all  the  experiments  conducted  the 
specimens  were  completely  saturated  at  peak  stress.  The 
computed  radial  strain  is  thought  to  be  an  upper  bound  as 


the  internal  pressure  is  not  distributed  over  the  entire 
length  of  the  confining  cylinder  due  to  some  overlap  on 
the  pressure  bars  (see  figure  3.8).  The  results  of  the 
computations  are  presented  in  table  E.5  in  appendix  E. 

The  computed  radial  strain  of  the  confining  cylinder 
for  all  experiments  did  not  exceed  1  percent.  Comparing 
the  radial  strain  of  the  confining  cylinder  (  cr )  to  the 
maximum  longitudinal  strain  (  Gj  )  experienced  by  the 
specimen  (see  table  E.3  in  appendix  B) ,  it  was  found  that; 


Gr<  0.06Cj  ,  (5.2) 

for  all  experiments  except  one  (experiment  137),  where  the 
radial  strain  was  9  percent  of  the  longitudinal  strain. 
Hence,  by  comparison,  the  specimen  can  be  considered  to  be 
in  a  state  of  nearly  uniaxial  strain  during  the 
experiment. 


5.4  Experimental  Replication 


If  the  SHPB  technique  is  to  prove  useful  for 
evaluating  dynamic  soil  response,  the  reproducibility  of 
experimental  results  must  be  established.  This  section 
will  present  the  results  of  replicate  experiments  on 
different  length  specimens  subjected  to  a  range  of  applied 
stresses . 

Replicate  experiments  were  conducted  at  the  nominal 
applied  stresses  of  250,  400,  and  520  MPa  and  specimen 
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lengths  of  6.35,  and  12.7  mm.  The  variation  in  the  applied 
stresses  did  not  exceed  ±  8  percent.  The  results  are 
presented  in  figures  5.3  through  5.8.  The  applied  stress 
is  given  in  parenthesis  next  to  the  experiment 
identification. 

Prom  figures  5.3  through  5.8  the  experimental  results 
can  be  considered  to  be  reproducible.  It  can  be  seen  that 
in  nearly  all  cases  the  slopes  of  the  stress-strain 
curves,  the  peak  stress,  and  the  strain  at  peak  stress  for 
the  replicated  experiments  are  directly  comparable  (the 
computed  values  of  peak  stress  are  tabulated  in  table  E.3 
in  appendix  E) .  However,  there  are  several  discrepancies 
to  be  noted.  In  figure  5.5  both  specimens  behaved 
similarily  up  to  275  MPa.  At  that  point,  the  specimen  in 
experiment  162  began  to  accumulate  strain  with  very  small 
changes  in  stress,  whereas  the  specimen  in  experiment  163 
continued  to  accumulate  stress  with  very  small  changes  in 
strain.  The  initial  specimen  lengths,  phase  relationships, 
and  applied  stresses  were  virtually  identical  for  both 
specimens  (see  table  E.l  in  appendix  E  for  phase 
relationships  and  table  E.3  in  appendix  E  for  applied 
stresses).  In  figure  5.7,  the  stress-strain  response  for 
each  specimen  is  very  similar  except  that  each  begins  to 
stiffen  at  different  strains. 

These  discrepancies  may  be  attributed  to  several 
factors.  One  is  the  difficulty  in  preparing  specimens  with 


116  (241  MPa) 


Figure  5.3.  Replicate  experiments  for  6.35  mm  specimen 
an  applied  stress  of  approximately  250  MPa 


112  (425  MPa) 


Figure  5.4.  Replicate  experiments  for  6.35  nun  specimens  at  an 
applied  stress  of  approximately  400  MPa. 
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of  approximately  400  MPa 
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identical  properties  (i.e.,  moisture,  density,  etc.).  A 
second  is  that  as  the  seating  strain  cannot  be  controlled 
with  great  precision  it  is  difficult  to  have  specimens 
with  identical  phase  relationships.  Third,  it  is  difficult 
to  achieve  duplicate  impact  velocities  due  to  a  variable 
friction  force  between  the  launch  tube  and  the  striker 
bar.  These  potential  sources  of  difficulty  in  replicating 
experimental  results  are  discussed  in  section  6.2.  Another 
factor  that  might  hamper  the  replication  of  experimental 
results  is  that  soil  and  or  moisture  might  be  escaping 
through  the  annular  region  between  the  confining  cylinder 
and  the  pressure  bars  during  the  experiment  (see  figure 
3.8).  A  discussion  of  soil  and  moisture  loss  is  presented 
in  section  6.4. 


5.5  Soil  Behavior 

Three  methods  of  analysis  will  be  used  to  describe 
the  observed  soil  behavior:  (1)  by  comparison  of  the 
stress-strain  response  of  specimens  with  similar  lengths 
to  a  range  of  applied  stresses,  (2)  by  comparing  the 
response  of  specimens  with  different  lengths  to  the  same 
applied  stress,  and  (3)  by  comparing  the  response  of 
specimens  with  similar  lengths  but  prepared  at  different 
initial  moisture  contents  to  the  same  applied  stress. 

Figures  5.9  and  5.10  show  the  typical  stress-strain 
response  for  the  nominal  specimen  lengths  of  6.35  mm  and 


Figure  5.9.  Stress-strain  response  for  6.35  mm  specimens 
to  a  range  of  applied  stresses. 
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12.7  am  respectively ,  to  a  range  of  applied  stresses.  It 
can  be  observed  that  the  average  stress  experienced  by  the 
specimen  increased  with  increasing  applied  stress 
independent  of  the  specimen  length.  For  both  specimen 
lengths ,  the  stress-strain  response  is  very  similar  for 
applied  stresses  up  to  400  MPa  with  some  increase  in 
stiffness  observed  at  the  higher  applied  stress.  For  all 
applied  stress  levels,  the  specimens  began  to  stiffen  at 
strains  approximately  equal  to  the  initial  gas  porosity. 
For  v  jpecimen  lengths  and  at  all  applied  stresses,  the 
strain  at  peak  stress  experienced  by  the  specimen  exceeded 
the  initial  gas  porosity  of  the  material. 

Figures  5.11,  5.12,  and  5.13  compare  the 
stress-strain  response  of  specimens  to  the  same 
approximate  applied  stress  but  with  different  specimen 
lengths.  The  nominal  applied  stresses  are  250,  400,  and 
520  MPa  respectively.  From  figures  5.11  and  5.12  it  can  be 
seen  that  the  strain  experienced  by  the  shorter  specimens 
exceeded  those  of  the  longer  specimens.  Also  the  longer 
specimens  showed  greater  stiffness  at  lower  strains  than 
the  shorter  specimens.  This  apparent  discrepancy  in 
response  between  the  two  specimen  lengths  may  be  explained 
by  examining  the  initial  gas  porosity  for  each  specimen. 

To  observe  how  moisture  content  variations  affect 
stress-strain  response,  specimens  were  prepared  at  the 
following  nominal  moisture  contents;  7,  13,  and  15 
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percent.  Figure  5.14  shows  the  effect  of  moisture  content 
on  the  stress-strain  response  of  6.35  nun  specimens  to  an 
applied  stress  of  250  MPa,  and  figure  5.15  shows  the 
response  of  12.7  mm  specimens  to  an  applied  stress  of  400 
MPa. 

These  figures  show  that  the  average  stress 
experienced  by  the  specimens  increased  while  the  strain  at 
peak  stress  decreased  with  increasing  moisture  content.  As 
with  the  other  stress-strain  curves  shown,  there  is  a 
marked  break  in  slope  near  a  strain  equal  to  the  initial 
gas  porosity.  This  change  in  slope  is  not  observed  for  the 
12.7  mm  specimen  with  the  lowest  moisture  content 
(experiment  135);  however,  the  maximum  strain  (16  percent) 
did  not  approach  the  initial  gas-filled  porosity  of  23.4 
percent.  Also,  the  specimens  became  somewhat  stiff er  with 
increasing  moisture  content,  at  least  at  strains  in  excess 
of  the  initial  gas  porosity.  At  strains  less  than  the 
initial  gas  porosity  this  observation  is  less  clear. 
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6.  DISCUSSION  OF  EXPERIMENTAL  RESULTS 


In  chapter  4  the  experimental  assumptions  that  need 
to  be  satisfied  in  order  to  obtain  meaningful  data  were 
presented.  Chapter  5  presented  the  stress-strain  response 
for  sand  specimens  compacted  to  different  moisture/density 
conditions  and  subjected  to  a  range  of  applied  stresses. 
The  stress-strain  response  wa3  computed  using  the  theory 
established  in  chapter  2  and  the  reduction  procedure 
outlined  in  chapter  3.  This  chapter  will  establish  that 
the  experimental  assumptions  stated  in  chapter  4  can  be 
satisfied  when  using  soil  specimens  and  will  address  the 
discrepancies  noted  in  section  5.5  concerning 
stress-strain  response.  In  addition,  the  specimen 
conditions  that  were  necessary  to  satisfy  the  assumptions 
and  those  that  could  contribute  to  erroneous  soil  response 
will  be  discussed.  Also,  the  advantages  of  the  SHPB 
experiment  over  similar  methods  to  determine 
one-dimensional  dynamic  soil  properties  will  be  addressed. 


6.1  Experimental  Assumptions 


From  the  strain  gauge  data  collected  during  a  SHPB 
experiment  it  is  assumed  that  the  stress,  strain,  and 
strain-rate  response  in  the  specimen  at  any  time  can  be 


computed  by  averaging  the  forces  and  particle  velocities 
at  the  specimen-bar  interfaces.  For  the  computation*  to  be 
meaningful  the  assumptions  of  the  experiment  must  be 
satisfied  (see  section  4.1).  The  degree  to  which  the 
assumptions  are  satisfied  using  soil  specimens  and  large 
diameter  pressure  bars  will  be  analyzed. 

6.1.1  Uniform  Distribution  of  Stress 

Using  the  criterion  developed  by  Davies  and  Hunter 
(1963); 

do  >  r2 Ps  _12  ,  (2.39) 

de  T  2 

an  estimate  can  be  established  as  to  when  stress 
uniformity  might  be  achieved  in  a  SHPB  specimen.  The 
inequality  is  conservative  because  its  derivation  is  based 
on  a  lower  permissible  value  for  the  propagation  velocity 
of  waves  in  the  specimen  (Davies  and  Hunter,  1963). 

Because  unloading  waves  generated  in  the  experiment  will 
travel  faster  than  the  plastic  waves  in  the  specimen, 
quasi-equilibrium  will  be  reached  sooner  than  indicated  by 
the  inequality  (Davies  and  Hunter,  1963). 

Figures  6.1,  6.2,  and  6.3  show  results  of  experiment 
134.  Using  equation  (2.39),  the  stress-strain  curve 
(figure  6.1),  the  appropriate  specimen  length  (1.259  cm), 
density  (2.12  q/crt?  )  ,  and  T  (2:130  microseconds)  stress 
equilibrium  is  estimated  to  be  reached  at  42  microseconds. 
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Figure  6.2.  Stress-time  response  of  the  specimen-bar  interfaces 
experiment  134. 
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Figure  6.3.  Average  specimen  stress-time  response  and  the 
difference  between  the  stress-time  response  o 
the  specimen-bar  interfaces,  experiment  134. 


The  values  for  specimen  length  and  density  are  given  in 
table  E.l  of  appendix  E  and  T  is  approximated  from  figure 
4.5.  Numerical  values  for  the  slope  of  the  stress-strain 
curve  for  experiment  134  are  given  in  the  example  computer 
output  in  appendix  G.  Figure  6.2  shows  the  stress 
histories  at  the  specimen-bar  interfaces.  From  figure  6.2 
it  can  be  seen  that  after  about  25  microseconds  the 
stresses  at  the  two  interfaces  begin  to  converge  to 
approximately  the  same  value  and  are  indeed  equal  at  50 
microseconds.  The  stress  difference  across  the  specimen  as 
a  function  of  time  and  the  average  stress-time  response  is 
shown  in  figure  6.3.  If  inertia  forces  were  absent,  the 
stress  difference,  would  be  zero.  It  can  be  seen  that 
after  25  microseconds  the  stress  difference  is  very  small 
relative  to  the  average  stress  from  which  it  can  be 
concluded  that  the  stress  gradient  in  the  specimen  has 
diminished. 

The  experimental  evidence  presented  in  figures  6.2 
and  6.3  support  the  estimate  computed  from  the  Davies  and 
Hunter  (1963)  criterion  as  to  when  stress  uniformity  is 
achieved  in  a  SHPB  soil  specimen.  Hence,  the  Davies  and. 
Hunter  criterion  can  be  used  to  estimate  when  stress 
uniformity  is  likely  to  be  achieved  in  a  soil  specimen. 

In  contrast  to  the  criterion  of  Davies  and  Hunter 
(1963)  which  provided  a  reasonable  estimate  to  the  time 
required  for  the  specimen  to  reach  equilibrium,  the 


criterion  of  Bertholf  and  Karnes  (1974)  is  violated  in 
every  experiment  by  at  least  a  factor  of  2  (see  section 
2.5).  For  example ,  a  maximum  strain-rate  of  1800  s1  was 
computed  for  experiment  134.  Substituting  1800  s*into 
equation  (2.40)  along  with  the  diameter  of  the  bar  (6.03 
cm)  yields; 

•  -l 

Demax*  10800  cm  s  ,  (6.3 


and  subtituting  the  rise-time  (S45  microseconds)  of  the 
loading  wave  into  equation  (2.41)  yields; 


-l 


7.5  ps  cm 


(6.4 


The  rise-time  of  the  loading  wave  can  be  approximated  from 
figure  4.5.  To  satisfy  the  Bertholf  and  Karnes  (1974) 
criterion  the  specimen  strain-rate  and  loading  wave 
rise-time  would  have  to  be  less  than  800  s*and  greater 
than  90  microseconds,  respectively . 

The  purpose  of  the  Pertholf  and  Karnes  (1974) 
criterion  was  to  bound  the  experimental  parameters  such 
that  the  influence  of  inertia  would  be  minimized.  Although 
the  experiments  in  this  research  did  not  satisfy  the 
criterion,  it  has  been  shown  that  stress  equilibrium  will 
be  achieved,  satisfying  the  uniform  stress  assumption. 


The  influence  of  inertia  and  its  affect  on  the 
ability  of  a  soil  specimen  to  reach  stress  equilibrium  has 
been  investigated  experimentally,  and  the  results  compared 
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with  the  estimates  computed  from  the  criterion  of  Davies 
and  Hunter  (1963)  and  Bertholf  and  Karnes  (1974). 
Experimentally  it  has  been  demonstrated  that  the  stress  at 
the  specimen-bar  interfaces  will  be  approximately  equal 
within  the  rise-time  of  the  incident  stress  wave.  This 
means  that  the  first  1  to  2  percent  strain  the  specimen 
sustains  must  be  disregarded  as  stress  equilibrium  within 
the  specimen  has  not  been  achieved.  It  has  been  shown  that 
the  criterion  of  Davies  and  Hunter  provides  a  reasonable 
estimate  as  to  when  stress  uniformity  is  achieved  in  the 
specimen  by  direct  comparsion  to  the  experi mental  results. 
In  contrast  to  the  criterion  of  Davies  and  Hunter,  the 
bounds  on  the  experimental  parameters  established  by 
Bertholf  and  Karnes  are  too  restrictive  when  applied  to 
soil  specimens. 

6.1.2  Wave  Dispersion 

Davies  (1948)  showed  that  oscillations  in  the 
recorded  stress  waves  were  attributable  to  dispersion  and 
could  be  predicted  (see  sections  2.4  and  4.1.2).  However, 
until  recently  there  has  been  no  technique  that  could 
easily  be  incorporated  into  the  standard  SHPB  data 
reduction  procedure  to  account  for  this  phenomenon 
(Follansbee  and  Frantz,  1983).  Consequently,  the 
oscillations  in  the  raw  data  were  retained  and  a  smoothing 
technique  applied  to  the  reduced  data.  Although  the 
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smoothing  of  reduced  SHPB  data  is  an  accepted  practice 
(e.g.,  Rand#  1967)#  a  systematic  method  of  accounting  for 
the  dispersive  nature  of  a  propagating  wave  is  preferable. 

The  effects  of  dispersion  on  the  experimental  results 
obtained  in  this  research  was  investigated  by  conducting  a 
test  with  the  end  of  the  incident  bar  being  a  free  end 
(see  section  4.1.2).  The  correction  procedure  was  shown  to 
account  for  the  wave  dispersion  occuring  between  the 
strain  gauge  location  and  the  free  end  of  the  incident 
bar. 

Figure  6.4  illustrates  the  effect  of  dispersion  on 
the  appearance  of  the  stress-strain  curve  for  experiment 

Here  it  can  be  seen  that  the  dispersion  corrected 
curve  is  smoother  and  easier  to  interpret.  For  instance# 
it  would  have  been  difficult  to  correlate  the  sudden 
increase  in  stiffness  to  the  condition  of  zero  air  voids 
had  the  dispersion  correction  procedure  not  been  employed. 
Hence#  an  effect  of  dispersion  on  experimental  results  is 
to  mask  trends. 

6.1.3  Stress  Variation  Over  the  Cross 
Section  of  the  Bar 

The  diameter  of  the  pressure  bars  used  in  this 
research  was  significantly  larger  than  that  used  in  past 
SHPB  experimentation.  A  consequence  of  using  a  larger 
diameter  bar  is  that  the  stress-time  history  over  the 
cross  sectional  area  of  the  pressure  bar  may  not  be 
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constant.  Hence,  measurements  made  at  the  surface  of  the 
bar  might  not  be  representative  of  the  behavior  over  its 
cross  section.  In  the  work  of  Davies  (1948)  he  showed  that 
if  R/A  <  0.1  the  displacements  at  the  bar  surface  would 
differ  from  those  along  the  axis  of  the  bar  by  ±  3 
percent. 

In  order  to  investigate  the  degree  to  which  the 
criterion  established  by  Davies  is  satisfied  for  the  SHPB 
used  in  this  research,  the  incident  wave  for  experiment 
134  was  decomposed  into  its  respective  Fourier  components. 
Table  6.1  lists  the  values  n,  R  /A  ,  Cn/Co#  end  the  ratios 
of  the  amplitudes  of  the  Fourier  component  nu>0  (Dn)  to  the 
largest  magnitude  Fourier  component  (A<>). 

For  n  <  9,  the  value  R/A  is  less  than  0.1,  and  the 
magnitude  of  the  Fourier  component  (Dn)  is  approximately  6 
percent  of  the  largest  magnitude  Fourier  component  (A<>). 
These  computations  illustrate  that  the  high  frequency 
Fourier  components  of  the  wave  are  damped  out  quickly 
which  tends  to  support  a  nearly  uniform  stress-time 
distribution  over  the  cross  sectional  area  of  the  bar. 
However,  n  •  17  is  required  to  model  the  experimentally 
observed  wave;  therefore,  the  error  in  the  computations 
will  be  in  excess  of  +,  3  percent. 

Figure  6.5  shows  the  incident  wave  modeled  with 
n  ■  17,  and  n  ■  9  along  with  their  difference.  The 
difference  was  taken  as  a  qualitative  indicator  of  the 
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Table  6.1 


Fourier  Coefficients,  Values  for  R/A  ,  and  Cn/Co 
for  the  Incident  Wave  of  Experiment  134 


n 

R/A 

Cn/Co 

Dn/A0  (a) 

1 

0.01170 

0.99966 

0.89275 

2 

0.02340 

0.99913 

0.65881 

3 

0.03500 

0.99864 

0.41980 

4 

0.04670 

0.99804 

0.18439 

5 

0.05840 

0.99733 

0.02868 

6 

0.07020 

0.99646 

0.12833 

7 

0.08200 

0.99538 

0.13652 

8 

0.09380 

0.99402 

0.11074 

9 

0.10570 

0.99229 

0.05668 

10 

0.11770 

0.99012 

0.02368 

11 

0.12980 

0.98743 

0.05596 

12 

0.14210 

0.98408 

0.05409 

13 

0.15450 

0.98002 

0.04164 

14 

0.16730 

0.97505 

0.02027 

15 

0.18030 

0.96911 

0.01707 

16 

0.19380 

0.96196 

0.02908 

17 

0.20770 

0.95353 

0.02016 

18 

0.22220 

0.94358 

0.00962 

19 

0.23750 

0.93186 

0.01329 

20 

0.25370 

0.91817 

0.00912 

a)  Ao  -  25.32032 


Figure  6.5.  Incident  stress  wave  modeled  with  n  =  17  and 
n  =  9  Fourier  terms. 
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error  .In  excess  of  3  percent.  A  quantitative  measure  can 
be  computed  by  summing  the  absolute  values  of  the  Fourier 
coefficients  for  n  equal  10  to  17  and  dividing  by  the 
summation  of  the  Fourier  coefficients  for  n  equal  0  to  17. 


\Dn 

nslo 


.  0.068 


17 

+  1 
n=j 


Dn 


(6.5) 


This  relative  error  will  be  an  upper  bound  as  Dn  is 
computed  from  An  and  Bn  which  can  be  opposite  in  sign. 
Hence,  as  a  worst  case  the  computational  error  will  be 
approximately  ±  10  percent. 

Figure  6.6  shows  the  difference  in  the  computed 
average  stress-strain  response  between  using  n  *  17  and 
n  »  9.  it  can  be  seen  that  the  stress-strain  response  is 
not  significantly  altered;  hence  ±  10  percent  can  be 
considered  to  be  an  upper  bound  on  the  error  when  assuming 
one-dimensional  wave  analysis  to  reduce  the  data  obtained 
using  60.3  mm  diameter  pressure  bars. 


6.1.4  Specimen-par  Interface  Friction 


The  numerical  work  of  Bertholf  and  Karnes  (1974)  has 
shown  that  friction  between  the  specimen  and  the  pressure 
bars  can  affect  specimen  response,  particularly  for 
relatively  thin  specimens.  Interface  friction  affects 
specimen  response  in  two  ways;  (1)  a  greater  stress  will 
be  required  to  produce  a  given  strain,  and  (2)  the 
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specimen  will  barrel.  From  their  work  they  concluded  that 
the  effects  of  friction  could  be  minimized  if  the 
specimen-bar  interfaces  were  lubricated  and  a  specimen 
aspect  ratio  of  approximately  0.5  selected. 

All  the  specimens  in  this  research  were  the  same 
diameter  as  the  pressure  bars  and  placed  in  a  confining 
cylinder  which  prevented  excessive  barreling  and  hence 
negligible  shear  stresses  between  the  specimen  and  the 
pressure  bars.  This  was  verified  in  section  5.3  where  it 
was  shown  that  the  specimen  did  sustain  some  radial  strain 
but  for  all  experiments  was  less  than  1  percent. 

The  results  of  section  5.5  can  be  used  to  examine  the 
effect  of  friction  on  specimen  stress-strain  response.  In 
figures  5.11/  5.12/  and  5.13  the  stress-strain  response  is 
compared  for  different  length  specimens  to  the  nominal 
applied  stresses  of  250,  400,  and  520  MPa,  respectively. 
These  figures  are  useful  because  the  specimen  lengths  are 
different  by  a  factor  of  two;  hence  if  friction  effects 
exist  their  stress-strain  response  should  not  be  the  same. 

In  all  three  figures  the  stress-strain  response  for 
the  different  length  specimens  is  nearly  identical  up  to 
their  initial  gas  porosities  (see  also  section  6.4).  For 
strains  in  excess  of  the  initial  gas  porosity  the  specimen 
response  is  governed  principally  by  the  pore  water.  This 
is  evident  from  the  results  presented  throughout  chapter 
5.  In  all  three  figures  the  peak  stress  of  the  shorter 
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specimens  exceeded  that  o£  the  lon9er  by  at  least  25 
percent.  As  the  specimen  does  sustain  some  radial  strain 
(see  section  5.3)  this  phenomenon  may  be  due  to  friction. 

However,  to  be  conclusive,  experiments  need  to  be 
conducted  on  specimens  with  as  similar  initial  gas 
porosities  as  possible. 

Hence,  for  soil  specimens  with  lengths  different  by  a 
factor  of  two  and  aspect  ratios  less  than  0.2,  friction  is 
not  observed  to  influence  the  stress-strain  response  at 

i 

strains  less  than  the  initial  gas  porosity.  Prom  this 
conclusion  it  can  be  determined  that  the  optimum  aspect 
ratio  of  0.5  is  too  restrictive  when  applied  to  soil 
specimens  and  that  by  maintaining  a  condition  of  nearly 
uniaxial  strain,  aspect  ratios  less  than  0.5  can  be  used 
with  the  effects  friction  on  stress-strain  response  held 
to  a  minimum. 

6.2  Experimental  Replication 

The  replicate  experiments  shown  in  section  5.4  were 
conducted  to  evaluate  the  degree  to  which  the  results  of 
an  experiment  could  be  reproduced.  In  this  research, 
several  aspects  of  the  experimental  process  have  been 
identified  as  hindering  complete  experimental 
reproducibility  with  compacted  sand  specimens. 

In  attempting  to  achieve  replicate  experimental 
results,  each  specimen  was  prepared  according  to  the  same 
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procedure.  The  procedure  for  preparing  each  specimen  has 
been  presented  in  section  3.3.  However,  regardless  of 
following  a  consistent  procedure  for  preparation, 
variations  in  specimen  parameters  oxisted  from  specimen  to 
specimen.  For  example,  the  optimum  moisture  content  for 
the  soil  was  13.3  percent  (see  section  3.3),  the  moisture 
content  for  the  specimens  prepared  near  optimum  ranged 
from  10.4  to  13.9  percent.  The  results  shown  in  section 
5.4  have  demonstrated  that  the  initial  parameters  will 
govern  the  specimen  stress-strain  response.  Additional 
evidence  will  be  presented  in  section  6.3.  If  complete 
experimental  replication  is  to  be  achieved,  then  each 
specimen  must  have  the  same  initial  parameters.  This 
requirement  has  been  difficult  to  attain  for  the  compacted 
specimens  prepared  in  this  research.  Therefore,  some 
discrepancies  in  the  results  are  inevitable. 

Due  to  the  nature  of  the  specific  SHPB  used  for  this 
research  and  the  compressibility  of  the  specimens,  the 
seating  strains  produced  by  positioning  the  specimen 
between  the  bars  could  not  be  controlled  with  great 
precision.  After  the  specimen  was  positioned  between  the 
bars,  and  a  tight  contact  between  system  components 
checked,  the  hydraulic  actuator  was  used  to  position  the 
momentum  trap  (see  figure  3.2).  During  this  process,  a 
variable  stress  was  applied  through  the  transmitter  bar  to 
the  specimen  causing  it  to  compress.  The  range  of  this 
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stress  was  from  2  to  8  MPa.  A  measure  of  the  seating 
strain  was  obtained  by  recording  the  displacement  the 
momentum  trap  caused  by  the  movement  of  the  hydraulic 
actuator  (see  section  5.1).  However  the  inability  to 
control  the  amount  of  seating  strain  made  it  difficult  to 
have  two  specimens  with  similar  lengths  and  phase 
relationships  (e.g.,  initial  gas  porosity).  This  increased 
the  difficulty  in  achieving  complete  experimental 
replication  and  contributed  to  the  observed  discrepancies, 
(e.g.,  see  figures  5.5  and  5.7). 

How  well  the  experimental  results  can  be  replicated 
are  also  dependent  on  how  consistently  the  load  is 
applied.  In  section  3.1  it  was  noted  that  the  magnitude  of 
the  applied  load  was  proportional  to  the  impact  velocity 
of  the  striker  bar  which  is  controlled  by  the  breech 
pressure  in  the  gas  gun.  To  determine  the  required  gas  gun 
breech  pressure  to  produce  a  given  striker  bar  impact 
velocity,  a  theoretical  velocity-breech  pressure 
calibration  curve  was  developed.  The  velocity-breech 
pressure  calibration  curve  for  the  0.254  ra  striker  bar  and 
its  construction  are  presented  in  section  B.2  of  appendix 
B.  Over  the  period  of  time  the  experiments  were  conducted, 
it  was  found  that  the  breech  pressure  to  produce  a  given 
impact  velocity  differed  from  that  which  was  predicted 
from  the  theoretical  curve.  The  variation  has  been 
attributed  to  a  variable  friction  force  between  the  launch 


tube  and  the  striker  bar.  To  minimize  the  effects  of 
friction ,  the  launch  tube  was  periodically  swabbed  out  and 
treated  with  a  light  coating  of  lubricating  oil.  Also  a 
dust  cover  was  fabricated  for  the  vent  holes  on  the  launch 
tube  so  that  contamination  could  be  reduced  when  the 
system  was  inactive.  Regardless  of  these  efforts,  it  was 
difficult  to  achieve  duplicate  impact  velocities  with  a 
given  breech  pressure. 

In  spite  of  the  above  mentioned  difficulties, 
experimental  replication  can  be  achieved  when  using  soil 
specimens  in  a  SHPB  experiment.  This  has  been  demonstrated 
convincingly  by  figures  5.3,  5.6  and  5.8. 

6.3  Soil  Behavior 

In  one-dimensional  compression  the  general 
stress-strain  response  exhibited  by  soil  is  S-shaped.  For 
small  stresses  changes,  yielding  is  observed  with  the 
stress-strain  curve  concave  to  the  strain  axis.  For  large 
stress  changes,  the  behavior  is  characterized  by 
stiffening  with  the  stress-strain  curve  being  concave  to 
the  stress  axis.  This  behavior  is  shown  in  figure  6.7.  The 
soil  is  a  dry  desert  alluvium  taken  from  a  site  near  Yuma, 
Arizona  (Oass  and  Bratton,  1983).  The  general  specimen 
stress-strain  response  observed  in  this  research  (see 
section  5.5)  is  consistent  with  this  description  and  is 
similar  to  that  found  by  other  investigators  who  performed 


\ 


4 

4 


I 


i 


134 


uniaxial  strain  experiments  on  similar  soils  but  at  lower 
loading  rates  (Jackson,  1968;  Calhoun  and  Kraft,  1966). 

In  section  5.5  three  observations  were  made 
concerning  figures  5.9  and  5.10:  (1)  the  peak  stress 
increases  with  increasing  applied  stress,  (2)  the 
specimens  began  to  stiffen  significantly  at  a  strain  near 
the  initial  gas  porosity,  and  (3)  the  strain  at  peak 
stress  consistently  exceeded  the  initial  gas  porosity.  The 
first  observation  is  consistent  with  anticipated  behavior, 
i.e.,  as  the  applied  stress  is  increased  the  peak  stress 
sustained  by  the  specimen  should  also  increase.  The  second 
observation  was  alt anticipated  and  will  be  analyzed 
below.  The  third  observation  was  not  anticipated  and  will 
be  analyzed  in  section  6.4. 

i 

In  figures  5.9  and  5.10,  at  strains  less  than  the 
initial  gas  porosity  the  compressibility  of  the  specimens 
is  nearly  constant.  In  figure  5.9  the  stress-strain 
response  for  the  specimen  of  experiment  163  does  not 
adhere  to  this  behavior.  However,  the  initial  gas  porosity 
of  the  specimen  was  50  percent  less  than  that  of  the  other 
two  so  that  a  stiffer  response  at  a  lower  strain  is  to  be  j 

expected.  At  strains  in  excess  of  the  initial  gas 
porosity,  the  compressibility  of  the  specimens  is  again 
nearly  constant  but  with  a  value  greater  than  the  initial 
compressibility.  For  example  in  figure  5.10  the  tangent 
modulus  for  experiment  134  at  4  percent  strain  is  0.2  GPa 
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whereas  at  11  percent  it  is  1.5  GPa.  It  should  be  noted 
that  the  tangent  modulus  at  11  percent  strain  is  nearly 
equal  to  the  bulk  modulus  of  water  which  is  approximately 
2  GPa.  Between  these  zones  of  nearly  constant 
compressibility  lies  a  transition  zone  that  can  be 
identified  through  the  initial  gas  porosity.  In  this  zone 
the  response  of  the  specimen  is  shifted  from  being 
governed  by  the  soil  mass  to  being  governed  by  the  pore 
water.  Hence,  the  initial  gas  porosity  of  a  specimen  can 
be  used  as  an  indicator  to  determine  a  change  in  the 
response  governing  mechanism.  It  is  also  of  interest  to 
note  that  a  similar  response  has  been  observed  for 
McCormick  Ranch  sand  subjected  to  high  hydrostatic 
compression  (Mazanti  and  Holland,  1970). 

Figures  5.11,  5.12,  and  5.13  compared  the  response  of 
different  length  specimens  to  a  range  of  applied  stresses. 
In  section  5.5  the  discrepancies  between  the  response  for 
the  two  specimen  lengths  were  noted  and  an  examination  of 
the  specimens  initial  gas  porosities  were  proposed  as  a 
means  to  explain  them.  In  figures  5.11  and  5.12  the 
initial  gas  porosity  of  the  longer  specimens  were 
approximately  50  percent  less  than  that  of  the  shorter 
specimens  (i.e.,  6.2  pecent  as  compared  to  9.0  percent, 
respectively  for  figure  5.10  and  5.8  percent  as  compared 
to  10.6  percent,  respectively  for  figure  5.11).  Hence, 
greater  stiffness  at  lower  strains  is  to  be  expected  for 


the  longer  specimens.  In  addition,  smaller  strains  should 
be  expected  for  specimens  with  lower  initial  gas 
porosities.  In  figure  5.13  the  situation  is  reversed  with 
the  initial  gas  porosity  of  the  shorter  specimen  less  than 
of  the  longer  (i.e.,  4.8  percent  as  compared  to  7.5 
percent,  respectively).  Hence,  the  behavior  shown  in 
figure  5.13  is  consistent  with  the  conclusion  that  the 
specimen  stress-strain  response  is  governed  by  the  initial 
gas  porosity.  In  addition,  it  can  be  concluded  that  the 
stress-strain  response  of  specimens  with  different  lengths 
to  the  same  applied  stress  will  be  similar  if  their 
initial  gas  porosities  are  similar. 

Figures  5.14  and  5.15  compare  the  stress-strain 
response  of  specimens  with  similar  lengths  but  prepared  at 
different  moisture  contents  to  the  same  nominal  applied 
stress.  The  response  shown  in  these  figures  for  both 
specimen  lengths  is  consistent  with  the  exception  of  the 
specimen  of  experiment  167.  However,  as  was  the  case  for 
experiment  163  discussed  above,  the  initial  gas  porosity 
of  the  specimen  is  approximately  50  percent  lower  than  the 
other  two  specimens.  Hence  the  response  is  not 
inconsistent  with  what  has  already  been  presented.  For  the 
specimens  prepared  near  and  wet  of  optimum,  the 
stress-strain  response  is  consistent  with  the  results 
presented  throughout  chapter  5;  hence  their  analysis  will 
follow  what  has  already  been  discussed.  Of  interest  in 
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figures  5.14  and  5.15  is  the  response  of  the  specimens 
that  were  prepared  at  moisture/density  conditions  dry  of 
optimum. 

In  figures  5.14  and  5.15  the  experiments  where  the 
specimens  were  prepared  dry  of  optimum  were  135  and  139. 

It  was  mentioned  in  section  5.5  that  a  change  in  slope  was 
not  observed  for  the  specimen  of  experiment  135  and  that 
the  strain  did  not  approach  the  initial  gas  porosity.  It 
can  also  be  seen  that  the  slope  of  the  stress-strain  curve 
is  nearly  constant  which  indicates  that  the 
compressibility  of  the  soil  skeleton  can  be  considered  to 
be  nearly  constant  and  that  the  compression  of  the  air 
voids  is  directly  proportional  to  changes  in  stress. 

A  difference  is  observed  in  the  response  of  the 
specimens  prepared  at  the  moisture/density  conditions  dry 
of  optimum  for  the  two  specimen  lengths.  In  the  12.7  mm 
specimen  (experiment  135),  which  had  a  gas  porosity  of 
23.4  percent,  very  little  build  up  of  stress  was 
experienced  at  the  maximum  strain  of  IS  percent.  For  the 
6,35  mm  specimen  (experiment  139),  which  had  a  gas 
porosity  of  17.7  percent,  a  substantially  higher  stress 
was  built  up  as  the  strain  approached  this  value.  Two 
factors  may  be  responsible  for  this  difference  in 
behavior.  First,  the  density  of  the  6.35  mm  specimen  (1.98 
g/cm3)  is  greater  than  that  of  the  12.7  mm  specimen  (1.87 
g/cm3 ) .  Second,  the  initial  gas  porosity  of  the  of  the 


6.35  mm  specimen  is  lower  than  that  of  the  12.7  mm 
specimen.  As  the  effect  of  density  variations  on  specimen 
response  was  not  investigated  it  can  not  be  concluded  that 
one  factor  is  more  dominate  than  the  other.  However ,  the 
respone  of  the  6.35  mm  specimen  is  consistent  with  the 
response  of  other  specimens  where  the  strain  approached  or 
slightly  exceeded  the  initial  gas  porosity  (e.g., 
experiment  119);  hence  supporting  the  conclusion  that  the 
initial  gas  porosity  is  a  governing  parameter  in 
determining  specimen  stress-strain  response. 

6.4  Strain  in  Excess  of  Initial  Gas  Porosity 


As  discussed  in  section  5.2  the  initial  gas  porosity 
is  an  approximate  upper  bound  on  the  strain  that  a 
specimen  can  sustain.  However,  the  results  presented  in 
section  5.5  showed  that  for  nearly  every  experiment,  the 
strain  sustained  by  the  specimen  exceeded  the  initial  gas 
porosity  (see  figure  6.8).  There  are  several  factors  which 
may  be  working  together  to  account  for  this  discrepancy: 

(1)  loss  of  soil  and  or  moisture, 

(2)  compression  of  the  pore  water  as  the  specimen 
becomes  saturated,  and 

(3)  radial  expansion  of  the  confining  cylinder  and 
specimen . 


initial  gas  porosity. 
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6.4.1  Loss,  of  soil  and  Moiature 

Soil  and  moiature  loss  are  possible  because  the 
diameter  of  the  bars  and  the  inside  diameter  of  the 
confining  cylinder  are  not  equal.  Although  the  tolerance  is 
very  small  (see  figure  3.8),  soil  and  moisture  loss  can 
occur  through  the  annular  region  between  them.  The 
mechanism  for  the  loss  is  that  a  pressure  gradient  is  set 
up  between  specimen-bar  interfaces  and  the  atmosphere. 

At  the  specimen-incident  bar  interface  a  portion  of 
the  initial  compressive  stress  wave  is  reflected  as  a 
tensile  wave  due  to  the  lower  impedance  of  the  specimen 
relative  to  the  bar.  This  tensile  wave  travels  back  down 
the  incident  bar  toward  the  end  at  which  the  impact 
occurred.  As  the  impact  end  of  the  incident  bar  is  now  a 
free  end,  the  tensile  wave  is  reflected  as  a  compressive 
wave  travelling  once  again  toward  the  specimen,  hence, 
reloading  the  specimen.  For  the  SHPB  experiment  the 
specimen  response  of  interest  is  that  due  to  the  initial 
stress  wave. 

This  multiple  impact  situation,  complicates  the 
ability  to  quantify  the  amount  of  soil  and  moisture  lost, 
as  the  quantity  lost  during  each  impact  is  indeterminate. 

In  addition,  the  specimen  may  sustain  additional  strain 
during  the  second  impact,  making  it  meaningless  to  compare 
a  measured  strain  with  the  computed  strain  which  results 
from  the  initial  stress  wave. 
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An  attempt  to  measure  the  mass  of  soil  lost  was  made 
in  several  experiments.  The  mass  of  soil  lost  was 
determined  by  recording  the  combined  mass  of  the  confining 
cylinder  and  specimen  before  it  was  positioned  between  the 
bars  and  just  after  it  was  removed,  at  the  completion  of 
the  experiment,  and  then  taking  the  difference.  The  maximum 
and  minimum  soil  loss  was  computed  to  be  7.9  and  1.7  g, 
respectively.  The  results  are  tabulated  in  table  E.6  of 
appendix  E.  These  measurements  overestimate  the  amount  of 
soil  lost  because  the  specimen  is  subjected  to  multiple 
impacts . 

Because  a  moisture  content  was  determined  at  the  time 
the  specimen  was  prepared  and  immediately  after  the 
experiment,  a  measure  of  the  amount  of  moisture  loss 
sustained  by  the  specimen  can  be  computed.  The  maximum  and 
minimum  moisture  losses  were  computed  to  be  27.6  and  1.7 
percent  of  the  initial  moisture  content,  respectively.  The 
moisture  loss  for  each  experiment  is  tabulated  in  table  E.7 
of  appendix  E.  If  all  the  moisture  loss  is  assumed  to  occur 
during  the  first  impact,  the  specimen  can  sustain  an 
additional  strain  proportional  to  the  volume  of  moisture 
lost.  Using  this  assumption,  the  strain  contribution  from 
moisture  loss  for  each  experiment  has  been  computed  and  the 
results  tabulated  in  table  6.2.  These  strains  will  be 
overestimates  for  the  same  reason  the  measured  soil  loss 


was  an  overestimate. 
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6.4.2  Compression  of  Pore  Water 


In  most  engineering  applications,  water  is  considered 
to  be  an  incompressible  fluid.  At  the  stresses  applied  in 
this  research  (>  100  MPa),  this  assumption  will  be  invalid. 
The  compression  of  the  pore  water  will  contribute  to  the 
strain  sustained  by  the  specimen. 

The  specific  volume  of  water  at  a  variety  of 
temperatures  subjected  to  a  range  of  applied  stresses  has 
been  tabulated  in  the  Handbook  of  Physical  Constants 
(Sydney,  1966).  These  values  have  been  used  to  compute  the 
strain  contribution  due  to  the  compression  of  the  pore 
water.  The  computation  results  are  presented  in  table  6.2. 
The  stress  used  in  the  computations  was  taken  as  the  peak 
stress  sustained  by  the  specimen. 


6.4.3  Radial  Expansion  of  the  Confininc 
Cylinder  and  Specimen 


It  has  been  shown  in  section  5.3  that  the  specimens 
are  constrained  to  a  state  of  nearly  uniaxial  strain 
because  of  the  small  amount  of  radial  strain  that  occurs. 
However  small,  the  radial  strain  will  make  some 
contribution  to  the  overall  specimen  strain.  The  method  by 
which  the  radial  strain  of  the  confining  cylinder  and 
specimen  has  been  computed  and  the  assumption  used  in  the 
computations  have  already  been  presented  in  section  5.3. 
The  strain  contribution  is  computed  by  accounting  for  the 
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change  in  specimen  area  resulting  from  the  radial 
displacement.  The  radial  strain  contribution  to  specimen 
strain  in  each  experiment  is  tabulated  in  table  6.2. 

6.4.4  Strain  Correction 

Table  6.2  sums  the  strain  contributions  from  water 
compression,  moisture  loss,  and  radial  strain  and  compares 
it  to  the  specimen  strain  at  peak  stress.  These  data  are 
plotted  in  figure  6.9.  It  can  be  seen  that  in  all  but  a  few 
experiments  a  balance  of  strain  still  remains.  Had 
sufficient  information  been  obtained  to  compute  the  strain 
contribution  from  soil  loss  the  remaining  balance  for  each 
experiment  would  be  less. 

Due  to  the  uncertainty  in  the  evaluation  of  the  strain 
contributions  from  the  above  mentioned  factors,  a  gross 
strain  adjustment,  to  account  for  the  excess  between  strain 
at  peak  stress,  and  the  initial  gas  porosity  has  been 
computed.  The  details  of  the  computations  are  presented  in 
appendix  J.  The  gross  strain  correction  requires  none  of 
the  assumptions  made  in  the  computations  of  the  strain 
contribution  from  moisture  loss,  pore  water  compression,  or 
radial  expansion.  It  is  based  solely  on  the  peak  stress 
sustained  by  the  specimen.  In  addition,  the  strain 
contribution  of  soil  loss  is  also  included. 

The  results  show  that  the  strain  balance  from  the 
gross  correction  is  less  than  that  computed  by  summing  the 
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individual  components  indicating  that  the  strain 
contribution  due  to  soil  loss  can  be  if  leant  (see 
figure  6.10).  Regardless  of  the  method  used,  an  accounting 
of  the  excess  between  the  strain  at  peak  stress  and  the 
initial  gas  porosity  can  be  made. 

6.5  Strain-Rate  Dependence 


s 


To  determine  the  strain-rate  sensitivity  of  the 
compacted  sand,  stress-strain-rate  curves  at  constant 
strains  were  constructed  for  each  nominal  specimen  length. 
The  data  were  taken  from  experiments  conducted  at  the 
nominal  applied  stresses  of  250  MPa  and  400  MPa.  The 
results  are  presented  in  figures  6.11  and  6.12  for  the 
6.35  mm  and  the  12.7  mm  long  specimens  compacted  at 
moisture  and  density  conditions  near  optimum,  respectively. 
The  details  of  the  construction  are  as  follows. 

For  each  experiment,  the  stress  and  strain-rate  for  a 
given  strain  was  recorded.  For  the  6 . 35  mm  specimens  the 
range  of  strains  was  from  5  percent  to  14  percent,  and  for 
the  12.7  mm  specimens  the  range  of  strains  was  from  5 
percent  to  10  percent.  The  stress-strain-rate  data  points 
from  the  experiments  for  a  given  strain  are  then  grouped 
according  to  the  applied  stress  at  which  the  experiment  was 
conducted  and  then  plotted.  From  figures  6.11  and  6.12  it 
can  be  seen  that  considerable  scatter  results.  To  simplify 
the  presentation,  the  group  of  stress-strain-rate  data 
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Figure  6.11.  Stress-strain-rate  plot  for  6.35  mm  specimens 
compacted  at  optimum  conditions. 


Figure  6.12.  Stress-strain-rate  plot  for  12.7  mm  specimens 
compacted  at  optimum  conditions. 


points  for  a  given  strain  and  nominal  applied  stress  were 
averaged.  In  figures  €.11  and  €.12  ,  the  dashed  lines  are 
the  average  stress-strain-rate  trajectories  for  a  given 
applied  stress.  The  solid  lines  connect  points  of  constant 
strain  between  the  stress-strain-rate  trajectories. 

In  this  type  of  plot,  the  degree  to  which  the 
compacted  sand  is  strain-rate  dependent  can  be  judged  by 
the  slope  of  the  constant  strain  curves  connecting  the 
stress-strain-rate  trajectories.  If  the  slope  is  zero,  it 
can  be  concluded  that  the  compacted  sand  response  is  not 
strain-rate  dependent.  The  results  shown  in  figures  6.11 
and  €.12  indicate  the  response  of  the  compacted  sand  is  not 
dependent  on  strain-rate  for  strains  below  the  initial  gas 
porosity.  This  is  consistent  with  the  results  reported  by 
Gaffney  et  al.  (1985)  who  found  no  strain-rate  dependence 
in  a  drier  but  otherwise  similar  soil  for  total  strains 
less  than  the  initial  gas  porosity. 

It  was  observed  in  figures  6.11  and  6.12  that  the 
stress  which  produced  a  given  strain  did  not  increase  with 
increased  strain-rate,  indicating  that  the  specimen 
response  was  not  dependent  on  strain-rate  for  strains  less 
than  the  initial  gas  porosity.  However,  this  apparent 
strain-rate  independence  should  be  viewed  with  caution  for 
two  reasons,  first,  the  constant  strain  curves  used  to  show 
this  apparent  strain-rate  independence  are  only  rough 
averages  developed  from  a  few  data  points,  and  second,  the 


factors  identified  to  account  for  the  discrepancy  between 
the  strain  at  peak  stress  and  initial  gas  porosity  cannot 
be  quantified  with  the  necessary  accuracy  to  determine  what 
their  effect  on  the  apparent  strain-rate  dependency  might 
be. 


6,6  Conditions  of  Experiment 

In  order  to  satisfy  the  assumptions  of  the  experiment 
(see  sections  4.1  and  6.1)  the  soil  specimens  were  prepared 
with  aspect  ratios  less  than  or  equal  to  0.2  and  contained 
in  a  confining  cylinder.  The  small  aspect  ratios  were 
necessary  to  satisfy  the  uniform  stress  assumption  while 
the  containment  facilitated  the  control  of  boundary 
conditions . 

If  the  specimens  had  been  longer  stress  equilibrium 
would  not  have  been  achieved  during  the  early  portion  of 
the  experiment.  Had  this  been  the  case,  a  greater  portion 
of  the  stress-strain  curve  would  have  had  to  be  discarded 
as  inertial  forces  would  have  dominated  the  response. 
Therefore,  soil  specimens  in  the  SHPB  experiments  are 
restricted  to  small  aspect  ratios.  It  may  be  possible  to 
use  greater  aspe_c  ratios  if  the  duration  of  the  applied 
stress  wave  is  extended  (e.g.,  by  using  a  longer  striker 
bar ) . 

In  the  traditional  SHPB  experiment  the  specimen  is  in 
a  state  of  nearly  uniaxial  stress.  Because  the  unconfined 
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strength  of  soil  (e.g.,  <  0.1  MPa)  is  insignificant 
compared  to  the  applied  stresses  (e.g.,  >  100  MPa), 
reproducing  the  traditional  SHPB  experiment  with  soil 
specimens  is  not  feasible.  Hence,  some  means  of  radial 
confinement  is  required.  Three  degrees  of  radial  support 
can  be  considered:  minimal,  intermediate,  and  stiff. 
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Providing  minimal  support  would  most  closely  approximate 
the  uniaxial  stress  state  in  the  traditional  experiment. 
However,  even  the  flimsiest  of  containment  systems  would 
provide  a  radial  restraint  that  would  be  significant 
relative  to  the  unconfined  strength  of  the  soil.  If  an 
intermediate  confinement  is  used,  radial  stress  and  strain 
histories  would  have  to  be  known  to  evaluate  specimen 
response.  Therefore,  the  use  of  minimal  or  intermediate 
containment  systems  would  complicate  the  data  analysis.  In 
contrast  to  the  first  two,  a  stiff  confinement  system  is 
easy  to  implement  and  has  the  advantage  that  it  constrains 
the  specimen  to  a  nearly  uniaxial  strain  path  (see  section 
5.3)  which  can  be  easily  duplicated  at  both  higher  and 
lower  strain-rates . 


6.7  Applications  of  the  SHPB  Technique 


In  recent  years  complex  constitutive  models  have  been 
developed  for  use  in  numerical  computations  to  predict  the 
response  of  soil  in  the  region  close  to  a  high  yield  weapon 
explosion  (e.g.,  Rubin  and  Sandler,  1977).  These  models  are 
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used  In  large  computer  codes  to  study  ground-shock  problems 
of  both  the  free-field  and  soil-structure  interaction  types 
(Nelson,  Baron,  and  Sandler,  1971).  If  these  models  are  to 
accurately  predict  field  response  they  must  be  fit  to  data 
which  reflect  the  type  of  confinement,  magnitude  of  stress 
change,  and  time  scale  of  stress  application  to  be  expected 
(Whitman,  1970).  These  requirements  have  placed  demands  on 
experimentalists  to  develop  laboratory  procedures  that  will 
reflect  field  conditions  and  provide  data  for  the 
evaluation  of  the  material  constants  which  appear  in  the 
models  (Nelson  et  al.,  1971). 

The  soil  test  most  often  selected  to  duplicate  the 
loading  condition  in  the  overpressure  region  associated 
with  weapons  effects  is  the  uniaxial  strain  test.  For  the 
past  30  years  the  dynamic  uniaxial  strain  device  has  been 
used  to  investigate  the  fundamental  nature  of  soil  behavior 
to  rapidly  applied  loads.  Although  current  devices  can 
apply  stresses  up  to  400  MPa  with  loading  times  of  0.3 
milliseconds  che  need  for  a  laboratory  test  that  more 
closely  simulates  the  environment  near  a  weapon  explosion 
is  clear.  The  SHPB  technique  is  a  means  to  satisfy  the 
demand  for  more  reliable  data  in  support  of  a  wide  variety 
of  weapon  effects  problems. 

It  is  common  practice  to  extrapolate  laboratory  data  to 
evaluate  the  material  constants  which  appear  in  a 
constitutive  model  that  will  be  used  to  represent  the  soil 
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behavior  to  an  explosive  event.  As  the  SHPB  technique  can 
be  used  to  evaluate  material  response  to  stresses  and 
strain-rates  that  are  closer  to  in  situ  conditons, 
extrapolation  can  be  avoided  and  the  material  constants 
determined  directly.  This  can  lead  to  a  decrease  in  the 
time  required  to  evaluate  the  adequacy  of  a  model  to  a 
particular  problem  as  well  as  a  better  model  because  the 
constants  can  be  evaluated  from  data  that  more  closely 
duplicates  field  conditions.  In  addition,  a  reduction  in 
field  costs  may  also  be  realized  as  extra  in  situ 
measurements  can  be  avoided  if  the  laboratory  data  can 
provide  a  closer  approximation  of  the  anticiapted  field 
conditions . 

The  flyer  plate  experiment  is  a  particular  example 
where  the  data  from  a  SHPB  experiment  would  prove  useful. 
The  purpose  the  flyer  plate  experiment  is  to  develop  a 
technique  for  field  calibrating  soil  stress  and  motion 
instrumentation  in  the  pressure  range  of  100  to  1000  MPa. 
In  this  experiment  a  planar  impact  is  applied  to  an 
instrumented  soil  test  bed  by  an  explosively  driven  metal 
plate  such  that  the  loading  path  is  approximately  uniaxial 
strain.  To  calibrate  the  instrumentation  the  soil  must  be 
properly  characterized  such  that  the  material  model  will 
accurately  predict  the  stress  and  motion  fields.  The  SHPB 
technique  is  the  most  appropriate  means  to  evaluate  the 
required  material  properties  as  precise  and  accurate 
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measurements  can  be  made  over  the  strews  rings  of  interest. 

The  SHPB  technique  can  also  be  used  to  extend  the 
knowledge  of  dynamic  soil  response  in  several  other  areas. 
In  section  6.4  it  was  described  that  during  the  experiment 
the  specimen  was  subjected  to  multiple  inpacts.  This 
phenomenon  may  be  exploited  to  investigate  soil  response  to 
multiple  loadings  in  rapid  succession  without  altering  the 
specimen  properties  in  between  loadings.  Hence,  providing  a 
laboratory  means  to  investigate  soil  response  to  a 
successive  impulse  loading  environment.  Another  area  of 
interest  is  the  effect  of  confinment  on  soil  response.  By 
varying  the  degree  of  radial,  support  to  a  soil  specimen  the 
effect  of  confinment  on  the  soil  stress-strain  response  can 
be  evaluated.  Such  information  would  be  valuable  to  the 
understanding  of  soil  response  to  airblast  loading  during 
the  transition  from  the  superseismic  region  to  the 
outrunning  region.  Also  the  SHPB  technique  is  capable  of 
providing  information  on  stress  attenuation  and  changes  in 
waveform  with  distance  and  their  influence  on  stress-strain 
response  as  a  function  of  stress  level. 
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7.  CONCLUSIONS  AND  RECOMMENDATIONS 


7.1  Conclusions 

The  analyses  contained  in  this  research  permit  the 
following  conclusions  regarding  the  use  of  the  SHPB 
experimental  method  to  measure  the  dynamic  response  of 
soil. 

(1)  For  compacted  sand  specimens  constrained  to  a  nearly 
uniaxial  strain  state/  the  assumptions  necessary  to 
obtain  meaningful  data  from  a  SHPB  experiment  can  be 
satisfied.  It  has  been  shown  that  stress  uniformity  over 
the  length  of  the  specimen  is  achieved  and  that  the 
criterion  of  Davies  and  Hunter  (1963)  can  be  used  to 
estimate  the  time  it  occurs.  The  criterion  of  Bertholf 
and  Karnes  (1974)  was  violated  for  nearly  every 
experiment  and  determined  to  be  to  restrictive  when 
applied  to  soil  specimens.  Wave  dispersion  was  shown  to 
occur  and  increase  the  difficulty  in  interpreting 
experimental  results.  However/  the  effects  of  wave 
dispersion  can  be  minimized  through  the  correction 
procedure  developed  by  Follansbee  and  Frantz  (1983).  The 
error  in  assuming  stress  uniformity  over  the  cross 
sectional  area  of  the  bar  was  determined  to  be  in  excess 
of  the  ±  3  percent  determined  by  Davies  (1948) 
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when  R/A  <0.1  but  less  than  ±  10  percent.  Friction 
effects  were  minimized  by  confining  the  specimen  and 
were  found  not  to  influence  stress-strain  response  for 
specimens  different  in  length  by  a  factor  of  two. 

(2)  Experimental  replication  can  be  achieved.  However, 
every  precaution  must  be  taken  to  prepare  specimens  with 
similar  phase  relationships  (e.g.,  initial  gas  porosity) 
and  to  conduct  the  experiments  following  the  same 
procedure  if  discrepancies  in  the  results  are  to  be 
avoided.  The  factors  that  were  found  to  complicate 
experimental  replication  are:  preparing  specimens  with 
similar  initial  parameters  (e.g.,  initial  gas  porosity), 
the  inablity  to  control  seating  strain  which  altered  a 
specimens  initial  parameters,  difficulty  in  achieving 
duplicate  impact  velocities,  and  soil  and  moisture  loss 
through  the  annular  region  between  the  confining 
cylinder  and  pressure  bars. 

(3)  The  stress-strain  response  is  governed  principally 
by  the  initial  gas  porosity  of  the  specimen.  At  strains 
less  than  the  initial  gas  porosity  the  specimen 
compresses,  closing  the  air  voids.  At  a  strain  near  the 
initial  gas  porosity  the  specimen  begins  to  stiffen 
significantly.  The  initial  gas  porosity  was  taken  as  an 
upper  bound  on  the  strain  a  specimen  could  sustain. 
However,  in  some  experiments,  the  specimen  continued  to 
strain  beyond  the  initial  gas  porosity.  The  factors 
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accounting  for  the  excess  strain  were  identified  as: 
soil  and  moisture  loss,  compression  of  the  pore  water, 
and  radial  expansion  of  the  confining  cylinder  and 
specimen.  This  was  verified  by  an  overall  strain 
correction  and  by  summing  the  strain  attributed  to  each 
factor. 

(4)  The  specimens  prepared  near  optimum  conditions  were 
insensitive  to  strain-rate  at  strains  less  than  the 
initial  gas  porosity.  At  strains  in  excess  of  the 
initial  gas  porosity  a  strain-rate  dependency  is-. 
apparent.  Pore  water  pressure  is  likely  to  play  a  role 
in  this  change,  but  how  is  not  clear. 

It  is  concluded  that  the  SHPB  experimental  method  can 
be  used  to  determine  the  dynamic  response  of  soil. 

Although  a  few  limitations  are  associated  with  using  soil 
specimens  in  a  SHPB  experiment,  there  are  no  serious 
obstacles  that  will  prevent  them  from  being  overcome.  The 
SHPB  method  is  a  relatively  simple  experiment  and  can 
significantly  extend  the  range  of  stresses  and 
strain-rates  that  can  be  applied  beyond  the  capabilities 
of  current  equipment  used  for  dynamic  soil  investigations. 

7.2  Recommendations 

The  success  of  this  research  effort  provides  a 
justification  for  the  expanded  use  of  the  SHPB  technique 
for  evaluating  dynamic  soil  response.  Consequently,  it  is 
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recommended  that  study  be  initiated  in  three  areas t  (1) 
alternative  measurements,  (2)  containment  systems,  and  (3) 
modification  to  equipment. 

As  a  supplement  to  the  traditional  measurement 
techniques  which  have  been  described  in  chapters  2  and  3, 
the  feasibility  of  sensing  the  stress  waves  directly  on 
the  specimen-bar  interfaces  should  be  explored.  For  foam 
specimens  such  a  technique  has  been  employed  that  uses 
quartz  crystals  mounted  on  the  ends  of  the  pressure  bars 
in  contact  with  the  specimen  (Hodge  and  Wasley,  1969). 

The  soil  containment  system  used  in  this  research  has 
presented  several  difficulties  the  most  important  of  which 
is  that  loss  of  soil  and  moisture  cannot  be  prevented. 
Therefore  it  would  be  advantageous  to  have  a  system  that 
would  prevent  such  losses.  Also  such  a  system  would  allow 
unloading  behavior  to  be  investigated.  In  addition,  if  a 
confining  cell  similar  to  those  used  for  triaxial  testing 
could  be  built  that  would  be  able  to  ply  confining 
stresses  of  several  hundred  megapascals,  a  means  would  be 
available  to  investigate  dynamic  shear  response.  Such  a 
containment  system  would  extend  the  capabilities  of  the 
method  when  using  soil  specimens  considerably.  As  a 
mimimum  however,  measurements  of  radial  stresses  and 
strain  should  be  attempted.  This  information  alone  would 
be  very  valuable  in  deriving  constitutive  models  for 
general  deformation  paths. 
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The  final  area  in  which  recommendation  will  be  made 
concerns  experimental  equipment.  Two  main  problems  were 
encountered  in  using  the  equipment  described  in  chapter  3 
when  using  soil  specimens:  (1)  the  inability  to  accurately 
control  seating  strains,  and  (2)  only  stresses  in  excess 
of  200  MPa  could  be  applied  with  any  degree  of 
consistency. 

The  first  difficulty  can  be  eliminated  by  the 
introduction  of  controls  sensitive  to  small  changes  in 
stress.  These  would  be  beneficial  for  two  reasons.  First, 
an  increased  control  on  seating  strains  could  be 
established  and  second,  a  preload  could  be  precisely 
controlled  to  simulate  geostatic  overburden  stresses. 

The  second  problem  may  be  overcome  by  using  a  striker 
bar  with  a  characteristic  impedance  less  than  that  of  the 
steel  pressure  bars.  As  the  applied  stress  is  proportional 
to  the  characteristic  impedance  of  the  striker  bar  (see 
appendix  C)  a  reduction  in  the  impedance  will  cause  a 
smaller  stress  to  be  applied  for  the  same  impact  velocity. 
For  example,  if  the  striker  bar  is  constructed  of  aluminum 
which  has  a  characteristic  Impedance  nearly  three  times 
less  than  that  of  steel,  for  the  same  impact  velocity  the 
applied  stress  would  be  reduced  by  nearly  a  factor  of 
three. 

A  final  recommendation  concerning  equipment  pertains 
to  the  pressure  bars.  The  ratio  of  the  characteristic 


161 


Impedances  for  the  soil  specimens  and  the  pressure  bars 
used  in  this  research  was  approximately  60  to  1.  Although 
the  uniform  stress  assumption  was  not  seriously 
compromised,  stress  equilibrium  in  the  specimen  would  have 
been  achieved  sooner  if  the  impedance  mismatch  had  been 
smaller.  Also  the  constancy  of  strain- rate  could  have  been 
inproved  if  the  characteristic  impedance  of  the  pressure 
bars  had  been  closer  to  that  of  the  soil  (Rand,  1967). 
Therefore  it  is  recommended  that  pressure  bars  be 
constructed  of  material  other  than  steel.  Potential 
candidates  may  be  aluminum  and  possibly  some  type  of 
plastic . 
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APPENDIX  A 


MATHEMATICS  OF  THE  DISPERSION 
CORRECTION  PROCEDURE 


A.l  General 


l 

The  procedure  for  correcting  a  waveform  for  j 

dispersion  can  be  broken  down  into  four  basic  operations 
(Follansbee  and  Frantz,  1983): 

(1)  transformation  of  the  original  waveform  to  a  Fourier 
series, 

(2)  computation  of  the  phase  velocity  for  each  frequency 
component, 

(3)  computation  of  the  adjustment  to  the  phase  angle 
based  on  position, 

(4)  reconstruction  of  the  original  waveform  with 
adjusted  phase  angles. 

A. 2  Transformation  to  a  Fourier  Series 

The  general  expression  for  a  Fourier  series  can  be 
written  as  follows  (Wiley  and  Barrett,  1982); 

oo 

f(t)  ■  An  +  I  An  cos  nu>ot  +  Bn  sin  nw0t  ,  (A.l) 

2  n=i 

where  f(t)  is  a  periodic  function,  uJo  is  the  wave 
frequency  defined  as; 


U)0*  2  Tt  ,  (A. 2) 

T 

where  T  is  the  period  of  the  wave,  and  A0/  An,  and  Bn  are 
the  Fourier  coefficients  which  can  be  expressed  as; 
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ln 


Bn 


fT 

Aq  -  2_  /f(t)  dt  , 

T  Jq 

rT 

2_  Jf(t)  cos  nui0t  dt, 
T  Jo 

rT 

2_  |f(t)  sin  n(u0t  dt. 
T  J0 


(A.  3) 

(A. 4) 

(A. 5) 


Upon  evaluation  of  the  Fourier  coefficients  A0,  An,  and 
Bn ,  the  general  Fourier  series  (equation  A.l)  is 
transformed  to  a  pure  Fourier  cosine  series.  The  reduction 
of  the  general  Fourier  series  involving  the  sum  of  a  sine 
and  cosine  term  is  accomplished  by  first  multiplying  the 
sum; 


An  cos  n  U)ot  +  Bn  sin  n  u>0t  , 


(A. 6) 


*V  An  *  Bn  i 
VAn  +  Bn 


(A. 7) 


to  yield. 


VAn  +  Bn 


An  cos  nuuo^  +  Bn  sin  n  cu0t 


LV  An  +  Bn 


VAn  +  Bn 


(A. 8) 


By  establishing  the  new  variables; 


Dn  *  VAn  +  Bn  , 


(A. 9) 


which  is  the  resultant  amplitude  of  the  frequency 
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component  nujot'  and; 

cos<p  ■  An  •  (A. 10) 

VAn  +  Bn 

sirup  ■  Bn  ,  (A. 11) 

"V  An  +  Bn 

where  <J>  is  the  phase  angle  which  is  a  measure  of  the  lead 
or  lag  of  the  nth  harmonic  in  reference  to  the  sine  or 
cosine  wave  of  the  same  frequency  (Wiley  and  Barrett, 

1982).  Expression  (A. 8)  can  be  rewritten  as; 

Dn  (cos<J>coa  nto0t  +  sin<psin  nu>0t).  (A. 12) 

Applying  the  following  trigonometric  substitution  to 
expression  (A. 12)  (Selby,  1975); 

cos  a  cos  b  +  sin  a  sin  b  -  cos  (a  -  b)  ,  (A. 13) 

yields. 


Dn  cos  (nuJot  -<p).  (A. 14) 

Substituting  expression  (A. 14)  into  equation  (A.l),  yields 
a  pure  Fourier  cosine  series  of  the  form; 

oo 

f(t)  -  Ao  +  Z  Dn  cos  (ncu0t  -<p)  .  (A. 15) 

2  n=1 

There  is  one  note  on  the  computation  <p  .  Instead  of 
computing  <p  as; 
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<p  -  tanT1^  ,  (A. 16) 

An 


<J>  will  be  computed  as; 

<p-  cos*1  An  .  (A.  17) 

Dn 


This  is  because  the  range  of  the  arctangent  is  ±ir /2, 
hence  for  some  phase  angles  (e.g.,  3  7T/4)  an  incorrect 
value  would  be  returned.  The  range  of  the  arccosine  is  0 
to  7r  .  To  maintain  the  proper  sign  of  the  phase  angle  in 
the  computations;  if  Bn  <  0/  the  resultant  phase  angle 
will  be  the  inverse  sign  of  the  value  obtained  in  equation 
(A. 17) . 


A. 3  Computation  of  Phase  Velocity 

The  propagation  velocity  of  a  wave  at  a  given 
frequency  is  called  its  phase  velocity  (Cn) .  Cn  can  be 
written  in  terms  of  cyclic  frequency  (f)  and  wavelength 
( A  )  as; 


Cn  -  fnA  ,  (A. 18) 

which  can  be  related  to  the  natural  frequency  (U)o) 
through; 


as, 


2  7T  fn  »  nu>0 


/ 


(A. 19) 
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Cn  *  nWftA  .  (A. 20) 

£  7T 

Equation  (A. 20)  can  be  written  in  dimensionless  form  by 
rearranging  and  dividing  both  sides  of  the  equation  by  C<>R 
to  yield. 

27T  Cn  R  “  R  n  u>n  .  (A. 21) 

A  Co  Co 

The  terms  on  the  right  side  of  equation  (A. 21)  are  all 
known  quantities.  The  only  unknown  in  equation  (A. 21)  is 
Cn>  A  value  for  cn  can  be  computed  using  the  polynomial 
approximation  (equation  2.43  to  the  fundamental  mode  of 
the  dispersion  equation  (equation  2.42).  The  two  variables 
in  the  polynomial  approximation  are  Cn/Co ,  and  R/A  .  By 
computing  Cn/C0  for  a  given  value  of  R/A  ,  the  quantities 
on  each  side  of  equation  (A. 21)  can  be  compared.  Hence,  by 
iterating  through  values  of  R/A  ,  a  value  for  Cn/C0  can  be 
computed  such  that  the  two  sides  of  equation  (A. 21)  will 
be  equal.  Once  equality  has  been  established,  Cn  can  be 
computed  directly. 


From  the  dispersion  equation  (equation  2.42)  it  has 
been  determined  that  Cn  is  dependent  on  A  (Davies,  1948). 
This  causes  a  wave  to  disperse  as  it  propagates  because 
the  high  frequency  components  will  lag  behind  the  low 
frequency  components.  Pollansbee  and  Frantz  (1983)  have 
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derived  a  relationship  based  on  the  phase  angle  «J»  to 
account  for  the  dispersion  that  occurs  as  a  wave 
propagates  a  distance  Ax.  The  derivation  of  the  phase 
angle  correction  is  shown  below. 

The  term  inside  the  brackets  of  equation  (A. 15)  can 
be  rewritten  in  terms  of  position  as; 

nuu0t  -  0  ■  nu>o  ft  +  Ax  -  Axl  ,  (A. 22) 

L  Co  CnJ 


where  Ax/Cn  is  the  time  for  a  given  component  to  travel  a 
distance  Ax.  Rewriting  the  terms  of  equation  (A. 22) 
yields ; 


n  tu0  t  -  <J> 


nw0 


m 

t  +  Ax 

'  i  -  c0' 

^7 

>  « 

(A. 23) 


nu>0  t  -  <p 


nuu0  t 


-  ng>n  A  x 

C0 


LCjl- 

Cn 


ll 


(A. 24) 


The  phase  angle  can  now  be  computed  as  a  function  of 
position  according  to; 


n  (juft  A  x  Cn  -  1  . 
C0  [Cn 


(A. 25) 


A. 5  Wave  Reconstruction 


The  last  operation  in  the  dispersion  correction 
procedure  is  wave  reconstruction.  Essentially  this  entails 
adding  or  subtracting  the  phase  angle  adjustment  computed 
for  a  given  position  Ax,  from  equation  (A. 25)  to  the 
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phaaa  angle  computed  at  the  original  position  x,  computed 
from  equation  (A. 17); 

<p  -  <p  ♦  <J>  .  (A.  26) 

The  value  for  <p  is  then  substituted  into  equation  (A. 15) 
for<f>,  and  the  wave  reconstituted  with  the  effects  of 
dispersion  accounted  for. 
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B.l  Gas  Gun  Operation 

A  schematic  o£  the  gas  gun  used  to  propel  the  striker 
bar  down  the  launch  tube  is  shown  in  figure  B.l.  The  three 
chamber  arrangement  is  prepared  for  firing  by  pressurizing 
with  nitrogen  gas  according  to  the  following  procedure. 

(1)  The  set  pressure  chamber  is  filled  to  set  the  piston 
assembly  and  prevent  any  leakage  of  gas  to  the  launch 
tube. 

(2)  The  firing  pressure  chamber  is  filled  to  a  pressure 
not  to  exceed  the  pressure  in  the  set  pressure  chamber 
so  that  the  risk  of  accidental  firing  is  reduced. 

(3)  The  breech  is  filled  to  a  pressure  such  that  when 
the  nitrogen  gas  is  released,  the  striker  bar  will  be 
propelled  down  the  launch  tube  at  a  predetermined 
velocity.  The  required  pressure  to  produce  a  given 
velocity  is  determined  from  the  velocity-breech  pressure 
calibration  curve  for  the  striker  bar  used.  To  assist  in 
experimental  replication,  the  breech  is  fitted  with  a 
pressure  transducer  that  is  connected  to  a  digital  volt 
meter  to  accurately  establish  pressurization  levels. 

(4)  Piring  is  accomplished  by  venting  the  set  pressure 
chamber  which  allows  the  gas  in  the  firing  pressure 
chamber  to  release  the  piston  assembly.  This  breaks  the 

j  seal  between  the  gun  and  launch  tube  causing  the  breech 

i 

|  pressure  to  be  released  behind  the  striker  bar 

!  propelling  it  down  the  launch  tube. 

» 

t 
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Figure  B.l.  Gas  gun  (Nagy  and  Muelenhaupt ,  1983). 
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B . 2  Velocity-Breech  pressure  Calibration 


The  impact  velocity  of  the  striker  bar  as  a  function 
of  gas  gun  breech  pressure  can  be  determined  from  the 
following  equation  (Nagy  and  Muelenhaupt,  1983); 


V(x)  -  -2PftVn 

m( 7  “  1  ) 


(B.l) 


where 


V ( x )  is  the  projectile  velocity  at  distance  x, 

X  is  the  distance  traveled  while  the  pressure  is 
applied, 

P0  is  the  initial  breech  pressure, 

V(j  is  the  breech  volume  (65,500  cm3), 

m  is  the  mass  of  the  projectile,  and 

7  is  the  ratio  of  specific  heats  for  the  driving  gas. 

The  calibration  curve  is  shown  in  figure  B.2 

B. 3  Velocity  Measuring  System 

The  velocity  measuring  system  consists  of  three  diode 
lasers  with  focusing  lenses  and  three  photodetectors.  The 
lasers  are  model  OL5100  3mW  continuous  wave  lasers  that 
operate  at  8000  to  8200  angstroms  in  the  infrared  and  are 
manufactured  by  Optical  Information  Systems.  The  focusing 
lenses  are  SELFOC-Micro  Lens  (SML)  made  by  Nippon  Sheet 
Glass  Company  Incorporated.  The  photodetectors  are 
manufactured  by  Hewlett-Packard.  A  photograph  of  the 
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velocity  measuring  system  detached  from  the  launch  tube  is 
shown  in  figure  B.3. 

The  lasers  with  the  focusing  lenses  and 
photodetectors  are  mounted  in  sets  across  from  each  other. 
Each  set  is  separated  a  distance  of  1.27  cm  from  the 
other.  As  the  striker  bar  passes  the  first  set,  a 
multichannel  clock  is  started.  Each  successive  set  stops  a 
channel  of  the  clock;  hence,  recording  the  time  for  the 
striker  bar  to  travel  across  each  of  the  two  1.27  cm 
intervals.  By  measuring  the  travel  time  across  a  given 
distance,  the  striker  bar  velocity  for  each  interval  is 
computed.  The  striker  bar  velocity  for  the  experiment  is 
then  obtained  by  averaging  the  striker  bar  velocity  of  the 
two  intervals. 

In  addition  to  measuring  the  striker  bar  velocity, 
the  velocity  measuring  system  is  also  used  to  locate 
within  an  8  millisecond  window  the  recording  of  the  strain 
gauge  data  by  the  digitizers.  This  is  accomplished  by 
initiating  a  time  delay  as  the  striker  bar  passes  the  last 
set  of  lasers  and  photodetectors.  The  time  delay  is  set 
according  to  the  desired  striker  bar  impact  velocity.  As 
the  digitizers  are  continuously  running,  the  time  delay 
establishes  a  termination  point  for  data  collection,  such 
that  only  the  strain  gauge  data  of  interest  are  captured 
within  the  8  millisecond  window. 


Figure  B.3.  Laser  velocity  assuring  system 
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B.4  Electro-Hvdraulic  Control  and  Monitor  Unit 

The  primary  circuits  of  the  electro-hydraulic  control 
and  monitor  unit  are  shown  in  figure  B.4.  The  symbols  Jl 
through  J15  refer  to  ports  on  the  rear  of  the 
electro-hydraulic  control  unit.  Front  and  rear  views  of 
the  electro-hydraulic  control  unit  are  shown  in  figures 
B.5  and  B. 6, respectively .  A  brief  description  of  each 
major  component  will  be  presented. 

B.4.1  Carrier  Oscillator  (Al) 

This  system  component  provides  the  carrier  signal  for 
the  transducer  bridges  and  the  demodulation  in  the 
transducer  preamplifiers  (Nagy  and  Muelenhaupt,  1983). 

Also  incorporated  is  an  adjustable  field-effect  transistor 
amplitude  control  circuit  that  limits  the  variations  in 
carrier  level  to  approximately  0.01  percent  (Nagy  and 
Muelenhaupt,  1983). 

B.4.2  Transducer  Preamplifiers  (A2.A3) 

There  are  two  carrier  preamplifiers,  one  for  axial 
load  (A2),  and  one  for  firing  pressure  (A3).  Each  one 
amplifies  the  signal  from  its  respective  transducer  and 
demodulates  it  to  produce  two  analog  signals  (Nagy  and 
Muelenhaupt,  1983).  Each  amplifier  operates  with  a 
carrier-driven  Wheatstone  bridge  circuit. 
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L- **J14  MASTER  OUT 

Figure  B.4.  Electro-hydraulic  control  block  diagram  and  monitor 
unit  circuits  (Nagy  and  MuelenhauDt,  1983). 


The  error  detector  is  an  integral  part  of  the 
hydraulic  actuator  control  system.  The  error  detector 
continuously  compares  the  feedback  signal  from  the  load 
control  amplifier  with  the  command  signal  and  produces  an 
appropriate  error  signal  to  be  used  by  the  servo-amplifier 
to  drive  the  servo-valve  (Nagy  and  Muelenhaupt,  1983). 

B.4.4  Servo-Amplifier  (A6) 

The  servo-amplifier  controls  the  current  to  the  „ 
servo-valve.  The  signal  from  the  load  control  (A7)  is 
summed  with  the  error  control  signal  and  converted  to  a 
current  by  the  servo-amplifier  to  drive  the  servo-valve. 

B.4.5  Load  Control  (A7) 

This  is  a  manual  control  to  adjust  axial  load 
magnitude.  It  is  recommended  that  some  compressive  force 
be  applied  to  the  system  with  the  hydraulic  actuator  to 
guarantee  contact  between  the  specimen  and  the  pressure 
bar3.  To  avoid  damaging  the  specimen  after  completing  the 
experiment,  a  bypass  circuit  has  been  installed  to 
maintain  control  over  the  actuator  at  all  times  (Nagy  and 
Muelenhaupt,  1983).  The  function  of  the  bypass  circuit  is 
to  prevent  the  actuator  from  applying  any  compressive 
force  to  the  specimen  after  the  experiment  is  completed. 


A  multistation  pushbutton  switch  allows  monitoring  of 
the  data  output  from  A2,  A3,  and  A7  (Nagy  and  Muelenhaupt, 
1983). 

B. 5  Dynamic  Signal  Conditioning  System 

The  dynamic  signal  conditioning  system  consists  of 
four  major  components: 

(1)  the  mainframe, 

(2)  eight  signal  conditioning  amplifiers, 

(3)  supply  voltage  regulation,  and 

(4)  calibration. 

B.5.1  Mainframe 

The  mainframe  consists  of  a  Vector  Model  CMA3B.20 
case  which  houses  eight  modular  amplifiers.  Power  is 
supplied  by  a  LND-Y-152  Lambda  Power  Supply  that  provides 
the  system  with  +  15  volts  regulated  power  (Nagy  and 
Muelenhaupt,  1983). 


B.5.2  Signal  Conditioning  Amplifiers 


operated  in  half  bridge  configuration.  Figure  B.8  shows 
the  wiring  for  the  half  bridge  operation.  The  bandwidth  of 
the  amplifiers  is  1  MHz ,  to  prevent  aliasing  of  the 
digital  records. 

B.5.3  Voltage  Regulation 

As  the  voltage  output  of  the  mainframe  is  +  15  volts 
and  the  requirement  of  the  system  is  10  volts,  voltage 
regulation  is  necessary  to  assure  proper  system 
performance  (Nagy  and  Muelenhaupt 1983). 

B.5.4  Calibration  Control 

The  system  is  designed  so  that  it  is  capable  of 
switching  to  a  calibration  signal  which  can  verify  the 
overall  gain  of  the  data  acquisition  system  (Nagy  and 
Muelenhaupt,  1983). 


APPENDIX  C 

DERIVATION  OF  THE  RELATIONSHIP  BETWEEN 
THE  STRIKER  BAR  IMPACT  VELOCITY  AND 
THE  MAGNITUDE  OF  THE  APPLIED 
STRESS  WAVE 
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The  relationship  between  the  striker  bar  impact 
velocity  and  the  magnitude  o£  the  applied  stress  wave  is 
accomplished  through  momentum  considerations  derived 
through  Kewton's  second  law  of  motion; 

P  ■  m  a  ,  (C.l) 

where  F  is  force,  m  is  mass,  and  a  is  acceleration. 

Equation  (C.l)  can  be  rewritten  by  making  the  following 
substitutions ; 

P  -  O  A  ,  (C . 2 ) 

a  -  _dy.  ,  (C .  3 ) 

dt 


where  O  is  stress,  A  is  the  cross  sectional  area  of  the 
striker  bar,  V  is  velocity,  and  t  is  time,  as; 


A  fa  dt  -  m  jdV  .  (C.4) 

The  limit  of  integration  t,  on  the  left  side  of  equation 
(C.4)  is  the  duration  of  impact  between  the  striker  bar 
and  the  incident  bar.  When  the  striker  bar  impacts  the 
incident  bar,  a  compressive  wave  is  generated  in  the 
incident  bar  and  striker  bar  as  well.  The  compressive  wave 
in  the  incident  bar  propagates  toward  the  sample.  The 
compressive  wave  in  the  striker  bar  propagates  in  the 
opposite  direction,  toward  the  free  end  of  the  striker 
bar.  Upon  reaching  the  free  end,  the  compressive  wave  is 
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reflected  as  a  tensile  wave  which  returns  to  .the  impact 
interface.  As  the  interface  is  unable  to  sustain  teneion, 
unloading  occurs ,  terminating  the  applied  stress  wave. 
Hence,  the  duration  of  impact  is  the  time  required  for  a 
wave  to  traverse  the  length  of  the  striker  bar  twice,  and 
can  be  expressed  as; 

t  «  2  L  ,  (C.5) 

Co 

where  L  is  the  length  of  the  striker  bar,  and  Co  is  the 
bar  velocity  of  the  striker  bar.  As  mentioned  previously, 
a  0.254  m  long  striker  bar  was  used  in  this  research, 
hence  the  duration  of  the  applied  stress  wave  is 
approximately  100  microseconds.  The  limit  of  integration 
V,  on  the  right  side  of  equation  (C.4)  is  the  striker  bar 
impact  velocity.  Carrying  out  the  integration  of  equation 
(C.4)  yields; 


a  A  2  L  -mV  ( C .  6  ) 

Co 

Substituting  PA  L  -  m  reduces  equation  (C.6)  to; 

o  - p  Co  v  .  ( c . 7 ) 

2 

The  magnitude  of  the  applied  stress  is  therefore 
determined  by  the  velocity  at  which  the  striker  bar 
impacts  the  incident  bar. 
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1)  Select  breech  pressure  to  produce  the  required  striker 
bar  impact  velocity. 

2)  Lubricate  specimen-bar  interfaces  with  a  thin  film  of 
silicon  grease. 

3)  Measure  specimen  length. 

4)  Position  specimen  between  pressure  bars. 

5)  Verify  that  the  momentum  trap  is  attached. 

6)  Move  momentum  trap  into  position  with  hydraulic 
actuator . 

7)  Record  specimen  seating  strain. 

8)  Check  for  continuous  system  contact. 

9)  Power  and  set  instrumentation. 

10)  Lock  exclusion  area. 

11)  Turn  on  warning  lights. 

12)  Open  nitrogen  gas  bottle  manual  valve. 

13)  Open  selected  nitrogen  gas  bottle  solenoid  valve. 

14)  Open  system  isolation  valve. 

15)  Retract  striker  bar  projectile. 

a)  Open  firing  chamber  vent. 

b)  Move  projectile  into  launch  tube  past  vent  holes. 

c)  Close  breech  vent. 

d)  Close  vacuum  valve. 

e)  Start  vacuum  pump. 

f)  Open  vacuum  valve  to  indicated  position  and  maintain 
vacuum  until  the  striker  bar  bumps. 

g)  Open  vacuum  to  full  open  position  for  5  to  10 
seconds . 

h)  Turn  off  vacuum  pump. 

i)  Close  vacuum  valve. 

j)  Open  breech  vent  to  eliminate  vacuum. 

k)  Close  breech  vent. 

l)  Close  firing  chamber  vent. 

16)  Pressurize  breech. 

a)  Set  "set  pressure"  to  300  psig. 

b)  Close  set  pressure  valves  and  allow  pressure  to 
stabilize. 
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c)  Pressurize  breech  to  desired  pressure  level. 

d)  Put  on  ear  protection. 

e)  Set  "firing  pressure"  to  200  psig. 

f)  Close  firing  pressure  valves  and  allow  pressure  to 
stablilize. 

17)  Evacuate  and  secure  exclusion  area. 

18)  Verify  instrumentation  ready. 

19)  Fire  gas  gun. 

a)  Power  fire  control  panel. 

b)  Set  standby  switch  to  "ARM"  position. 

c)  Press  "FIRE"  switch. 

d)  Set  switch  to  "STANDBY"  position. 

e)  Turn  fire  control  panel  off. 

20)  Close  nitrogen  gas  bottle  solenoid  valve. 

21)  Vent  nitrogen  gas  bottle  manifold. 

22)  Close  system  isolation  valve. 

23)  Close  nitrogen  gas  bottle  manifold  vent  valve. 

24)  Enter  exclusion  area. 

25)  Close  nitrogen  bottle  manual  valve. 

26)  Open  breech  vent. 

27)  "SLOWLY"  vent  firing  chamber  pressure. 

28)  Turn  off  warning  lights. 

29)  Open  exclusion  area. 

30)  Turn  off  laser  velocity  system  power  supply. 


< 


31)  Remove  specimen  from  between  pressure  bars. 

32)  If  possible  measure  final  specimen  strain. 

33)  Remove  portion  of  specimen  for  moisture  content 
determination . 


it 
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Table  E.l 

Tnitial  Specimen  Parameters 


ftqjariant 

Wo. 

Ungtb 

imi 

Noiitur* 

Oontant 

(t) 

Qm 

Vorosity 

<%> 

Nat 

ttaMity 
(fl/oc  ) 

Void 

tmtlo 

D*gr«*  of 
Saturation 
(%> 

31 

1.283 

U.9 

7.07 

2.11 

0.42 

76.0 

32 

1.216 

11.9 

0.00 

2.31 

0.29 

100.0 

38 

1.230 

13.2 

19.01 

1.81 

0.67 

52.6 

39 

1.232 

13.2 

9.33 

2.03 

0.49 

71.7 

40 

1.201 

13. S 

4.00 

2.14 

0.42 

86.4 

41 

1.235 

13.3 

7.22 

2.07 

0.46 

77.1 

42 

1.198 

13.4 

4.47 

2.13 

0.42 

84.9 

55 

1.235 

12.4 

6.45 

2.U 

0.42 

78.3 

56 

1.207 

12.6 

6.79 

2.10 

0.43 

77.6 

57 

1.252 

12.7 

15.31 

1.90 

0.58 

58.4 

58 

1.261 

13.9 

9.76 

2.00 

0.52 

71.4 

59 

0.584 

12.7 

6.59 

2.10 

0.43 

78.2 

60 

0.703 

12.8 

22.97 

1.73 

0.74 

46.1 

61 

0.549 

13.1 

4.18 

2.14 

0.41 

85.6 

62 

0.552 

12.9 

9.76 

2.02 

0.49 

70.3 

63 

0.537 

13.3 

U.22 

1.98 

0.53 

67.5 

70 

1.216 

12.1 

14.29 

1.94 

0.54 

59.4 

71 

1.228 

11.7 

11.91 

2.00 

0.49 

63.4 

72 

0.549 

15.4 

0.90 

2.16 

0.42 

97.0 

73 

0.611 

12.2 

U.40 

1.98 

0.51 

63.5 

74 

0.552 

12.0 

7.09 

2.10 

0.42 

76.1 

75 

0.555 

12.3 

11.87 

1.99 

0.51 

64.7 

76 

1.240 

U.9 

14.16 

1.95 

0.54 

59.4 

112 

0.644 

U.8 

7.91 

2.09 

0.43 

73.6 

113 

0.645 

12.1 

7.73 

2.09 

0.43 

74.5 

114 

0.655 

12.1 

9.21 

2.05 

0.46 

70.7 

115 

0.645 

U.4 

10.55 

2.04 

0.46 

66.4 

116 

0.637 

U.2 

8.07 

2.10 

0.41 

72.4 

117 

0.635 

10.7 

8.50 

2.10 

0.41 

70.5 

118 

0.645 

10.6 

10.25 

2.07 

0.43 

65.9 

119 

0.645 

10.4 

9.91 

2.08 

0.42 

66.4 

131 

1.313 

U.4 

9.57 

2.04 

0.47 

70.2 

132 

1.269 

U.5 

6.88 

2.10 

0.43 

77.2 

133 

1.291 

U.4 

8.02 

2.07 

0.45 

74.1 

134 

1.259 

U.4 

5.80 

2.U 

0.41 

80.2 

135 

1.307 

7.0 

23.37 

1.84 

0.55 

34.1 

136 

1.223 

15.1 

5.17 

2.08 

0.48 

84.1 

137 

0.665 

15.4 

8.51 

2.00 

0.54 

75.8 

138 

0.596 

13.0 

4.86 

2.13 

0.42 

83.5 

139 

0.609 

7.0 

17.68 

1.98 

0.44 

42.3 

14; 

1.265 

U.8 

6.23 

2.13 

0.40 

78.3 
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Table  E.l  cont. 


toper  iant 

NO. 

^length 

(<a) 

Maietora 

Oartant 

(«) 

Om 

tensity 

<%> 

Ntt 

Density 
(9/ac  ) 

fi 

ngm  of 

Saturation 

(«) 

146 

1.267 

11.9 

6.26 

2.13 

0.41 

78.3 

147 

0.636 

11.9 

6.64 

2.12 

0.41 

77.2 

148. 

0.631 

14.4 

7.41 

2.04 

0.50 

77.6 

162 

0.629 

13.0 

4.44 

2.14 

0.41 

84.7 

163 

0.631 

12.9 

4.84 

2.13 

0.41 

83.4 

164 

1.298 

12.9 

7.49 

2.07 

0.45 

76.0 

165 

1.289 

12.7 

7.05 

2.09 

0.44 

77.0 

166 

0.624 

14.0 

4.42 

2.12 

0.44 

85.5 

167 

0.622 

14.0 

4.15 

2.12 

0.43 

86.3 
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Table  E.2 

Moisture  Loss  Between  the  Time  a  Specimen  is  Prepared 
and  the  Time  it  is  Positioned  Between  the  Bars 


Experiment  No. 

MCI (a) 

(%) 

MC2(b) 

(%) 

Percent  Difference 

(-) 

A 

13.4 

12.3 

8.2 

B 

13.0 

12.2 

6.2 

C 

12.9 

12.2 

5.4 

D 

13.0 

12.7 

2.3 

E 

13.1 

12.6 

3.8 

F 

13.2 

12.6 

4.5 

G 

13.2 

12.4 

6.1 

H 

13.3 

12.4 

6.8 

Arithmetic  mean  *  5.4 

Standard  Deviation  -  1.8 

a)  MCI  is  the  specimen  moisture  content  at  the  time  it  was 
prepared. 

b)  MC2  is  the  specimen  moisture  content  at  the  time  it  was 
positioned  between  the  pressure  bars. 

c)  The  specimens  for  these  tests  were  prepared  following 
the  procedure  that  was  outlined  in  section  3.3.  The 
same  sequence  of  steps  from  specimen  preparation  to 
just  before  firing  of  the  gas  gun  were  followed  as  if 
an  actual  experiment  were  to  be  performed. 


Table  E.3 


Experimental  Results 


feqpariMnt  lie. 

*»U«1 

(Ml) 

tak  8traM 
<Wa> 

Strain  at 

(t) 

Strain-Mat*  at 

Saak  Strata 

(  rt 

31 

135 

10 

5.3 

266 

32 

555 

379 

15.0 

667 

X 

SSI 

273 

18.8 

502 

39 

307 

134 

11.0 

140 

40 

392 

168 

14.0 

536 

41 

463 

205 

15.8 

566 

42 

253 

S3 

9.2 

669 

SS 

457 

248 

14.3 

718 

56 

374 

220 

12.8 

315 

57 

365 

129 

14.0 

442 

58 

301 

89 

11.3 

345 

59 

489 

349 

19.2 

1070 

60 

509 

371 

14.8 

781 

61 

260 

168 

14.4 

537 

62 

270 

247 

10.9 

150 

63(a) 

365 

••• 

■••• 

70(b) 

629 

»** 

•••• 

•••• 

71 

631 

340 

16.2 

1019 

72 

390 

164 

29.7 

904 

73 

416 

285 

18.4 

606 

74 

764 

333 

30.9 

2159 

75 

261 

115 

18.8 

1077 

76 

260 

29 

10.2 

334 

112 

425 

253 

19.8 

1327 

113 

385 

247 

17.4 

1062 

114 

386 

250 

19.6 

678 

115 

395 

260 

18.0 

1051 

116 

243 

153 

14.2 

424 

117 

268 

208 

14.5 

294 

116 

244 

160 

14.1 

313 

119 

246 

153 

14.1 

389 

131 

367 

194 

12.7 

440 

132 

375 

203 

10.7 

673 

133 

368 

172 

12.5 

472 

134 

399 

209 

12.2 

696 

135 

385 

36 

15.6 

382 

136 

397 

283 

11.2 

408 

137 

249 

239 

6.3 

90 

136 

251 

184 

11.6 

579 

139 

249 

114 

17.2 

347 

145 

269 

105 

10.1 

ISO 

146 

523 

313 

13.3 

867 

147 

237 

147 

13.0 

585 

_ _ ’iW b  •£ fr^vl-vv.!  •>  Ir-v-v  ••I.Iva’C*-’! 
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Table  E.3  cont 


teperl— it  No. 

Applied 

Stress 

(We) 

Peek  Stress 
(We) 

Strain  at 

Peek  Stress 

(«> 

Strain  Nets  et 

"wr* 

148 

249 

117 

15.2 

972 

162 

S19 

325 

20.6 

1327 

163 

522 

425 

15.3 

459 

164 

507 

309 

12.4 

696 

165 

254 

76 

9.3 

136 

166 

261 

124 

13.3 

967 

167 

260 

151 

12.6 

779 

a)  In  experiment  63  the  gain  setting  for  the  txanmaitter  bar  strain  gauge  was 
incorrectly  aet.  This  caused  the  digitisers  collecting  the  strain  gauge  data  to  be 
overdriven.  As  a  result,  oily  two  of  the  required  three  data  aeta  for  data  reduction 
ware  collected. 

b>  In  experiment  70  e  tight  fit  between  ell  system  opponents  was  not  maintained.  Aa  a 
result*,  an  irregular  incident  stress  wave  was  initiated  that  was  inauitable  for  data 
reduction. 


i 


i 


i 
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Table  B.4 

Specimen  Seating  Strain 


Experiment 

No. 

Preplacement 

Length 

(cm) 

In- Position 
Length 
(cm) 

Seating 
Stra in 
(%) 

112 

0.673 

0.644 

4.3 

113 

0.674 

0.645 

4.3 

114 

0.723 

0.655 

9.4 

115 

0.677 

0.645 

4.7 

116 

0.715 

0.637 

10.9 

117 

0.697 

0.635 

8.9 

118 

0.694 

0.645 

7.1 

119 

0.721 

0.645 

10.5 

131 

1.349 

1.313 

2.7 

132 

1.338 

1.269 

5.2 

133 

1.367 

1.291 

5.6 

134 

1.382 

1.259 

8.9 

135 

1.337 

1.307 

2.2 

136 

1.330 

1.223 

8.0 

137 

0.670 

0.665 

0.7 

138 

0.688 

0.596 

13.4 

139 

0.728 

0.609 

16.3 

145 

1.360 

1.265 

7.0 

146 

1.316 

1.267 

3.7 

147 

0.759 

0.636 

16.2 

162 

0.646 

0.629 

2.6 

163 

0.684 

0.631 

7.7 

164 

1.362 

1.298 

4.7 

165 

1.345 

1.289 

4.2 

166 

0.687 

0.624 

9.2 

167 

0.654 

0.622 

4.9 

A 

0.727 

0.641 

11.8 

B 

0.711 

0.623 

12.4 

C 

0.715 

0.631 

11.7 

D 

1.395 

1.273 

8.7 

E 

1.364 

1.275 

6.5 

F 

1.369 

1.273 

7.0 

G 

1.335 

1.273 

4.6 

Arithmetic  mean  »  7.5 
Standard  Deviation  ■  3.9 

a)  Experiments  numbers  A  through  G  do  not  represent 
complete  experiments.  These  were  tests  performed  to 
provide  additional  information  on  seating  strains  and 
moisture  loss. 
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Table  E.5 


Computed  Radial  Displacements  Due  to  Applied  Load 


Experiment  No.  Inside  Radius  of  Cylinder  Displacement 

(cm)  (cm) 


31 

3.017 

0.0007 

32 

3.018 

0.0264 

38 

3.025 

0.0191 

39 

3.030 

0.0094 

40 

3.059 

0.0120 

41 

3.055 

0.0147 

42 

3.018 

0.0037 

55 

3.048 

0.0177 

56 

3.030 

0.0154 

57 

3.058 

0.0093 

58 

3.055 

0.0064 

59 

3.058 

0.0251 

60 

3.028 

0.0261 

61 

3.033 

0.0118 

62 

3.020 

0.0172 

63 

3.045 

****** 

70 

3.030 

***  **  * 

71 

3.032 

0.0240 

72 

3.063 

0.0118 

73 

3.035 

0.0201 

74 

3.053 

0.0238 

75 

3.038 

0.0082 

76 

3.018 

0.0020 

112 

3.060 

0.0182 

113 

3.056 

0.0177 

114 

3.026 

0.0175 

115 

3.060 

0.0187 

116 

3.060 

0.0111 

117 

3.056 

0.0149 

118 

3.059 

0.0115 

119 

3.040 

0.0109 

131 

3.045 

0.0138 

132 

3.062 

0.0146 

133 

3.060 

0.0124 

134 

3.055 

0.0150 

135 

3.030 

0.0025 

136 

3.020 

0.0198 

137 

3.018 

0.0166 

138 

3.025 

0.0129 

139 

3.025 

0.0080 

145 

3.025 

0.0073 

146 

3.028 

0.0220 

V 


s 


> 


•I 
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Table  G.5  cont. 


Experiment  Mo. 

Inside  Radius  of  Cylinder 
(cm)‘ 

Displacement 

(cm) 

147 

3.056 

0.0105 

148 

3.035 

0.0083 

162 

3.028 

0.0228 

163 

3.023 

0.0297 

164 

3.025 

0.0217 

165 

3.055 

0.0053 

166 

3.035 

0.0087 

167 

3.028 

0.0106 
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Table  E.6 

Mass  of  Soil  Lost  Daring  the  Experiment 


Experiment  Mo. 

- 

Mass  Lost 

<g> 

138 

1.7 

145 

2.1 

147 

1.9 

148 

5.3 

162 

7.9 

164 

4.8 

165 

2.1 

166 

4.8 

167 

2.3 

Table  E.7 


Measured  Moisture  Content  Changes  Before 
and  After  the  Experiment 


Experiment  No.  Preexperiment  -  Postescperiment  Percent 

Moisture  Content  Moisture  Content  Difference 


(%) 

(%) 

(-) 

31 

11.9 

10.0 

16.0 

32 

11.9 

10.8 

9.2 

38 

13.2 

12.2 

7.6 

39 

13.2 

12.9 

2.3 

40 

13.5 

13.1 

3.0 

41 

13.3 

10.7 

19.5 

42 

13.4 

9.7 

27.6 

55 

12.4 

**** 

56 

12.6 

**** 

**** 

57 

12.7 

**** r 

**** 

■i 

13.9 

**** 

**** 

59 

12.7 

**** 

**** 

60 

12.8 

**** 

61 

13.1 

**** 

it*  it* 

62 

12.9 

**** 

**** 

63 

13.3 

It*  irk 

**** 

70 

12.1 

10.3 

14.9 

71 

11.7 

9.2 

21.4 

72 

15.4 

8.5 

44.8 

73 

12.2 

10.5 

13.9 

74 

12.0 

11.8 

1.7 

75 

12.3 

10.0 

18.7 

76 

11.9 

10.3 

13.4 

112 

11.8 

10.4 

11.9 

113 

12.1 

10.0 

17.4 

114 

12.1 

10.1 

16.5 

115 

11.4 

10.2 

10.5 

116 

11.2 

9.8 

12.5 

117 

10.7 

9.9 

7.5 

118 

10.6 

9.6 

9.4 

119 

10.4 

9.5 

8.7 

131 

12.4 

11.4 

8.1 

132 

12.5 

11.7 

6.4 

133 

12.5 

11.0 

12.0 

134 

12.4 

11.4 

8.1 

135 

7.0 

6.2 

11.4 

136 

15.1 

11.5 

23.8 

137 

15.4 

12.1 

21.4 

138 

13.0 

12.0 

7.7 

139 

7.0 

6.3 

10.0 

145 

11.8 

11. 3< 

4.6) 

4.2 

146 

11.9 

11. 1< 

6.2) 

6.7 

147 

11.9 

10. 7< 

6.4) 

10.1 

v  r_-  ‘.t  *-•  ■ 
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Table  E.7  cont. 


Experiment  No. 

Preexperiment 
Moisture  Content 
(%) 

Postexperiment 
Moisture  Content 

(%) 

Percent 

Difference 

(-) 

148 

14.4 

12.8(11.4) 

11.1 

162 

13.0 

11. 7(  6.6) 

10.0 

163 

12.9 

11.9 

7.8 

164 

12.9 

11. 9(  7.4) 

7.8 

165 

12.7 

12.1(11.1) 

4.7 

166 

14.0 

12.7(11.5) 

9.3 

167 

14.0 

12.7(11.0) 

9.3 

a)  The  numbers  in  parentheses  beside  the  post  experiment 
moisture  content  percentages  are  moisture  contents  ( in 
percent)  of  the  material  scraped  from  the  region  where  the 
bars  and  confining  cylinder  overlap  (see  figure  3.8). 
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APPENDIX  P 


PROCEDURE  AND  RESULTS  OP  MOISTURE/DENSITY 
VARIATION  STUDY 
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The  soil  used  for  this  study  was  the  same  as  that 
used  in  the  SHPB  experiments.  The  specimens  were  prepared 
at  moisture/density  conditions  near  optimum  as  determined 
from  the  Harvard  miniature  compaction  procedure  (see 
figure  3.5).  The  specimens  for  this  study  were  prepared 
following  the  same  procedures  outlined  in  section  3.3. 

The  specimens  were  compacted  in  stainless  steel 
cylinders  with  the  nominal  dimensions  of  12  cm  in  length 
and  6  cm  inside  diameter,  static  compaction  of  the 
specimen  was  performed  from  both  ends  of  the  cylinder. 

Pour  lifts  of  soil  were  used  to  complete  the  compaction 
processs.  After  compaction ,  the  specimen  was  extruded  from 
the  cylinder  and  sliced  into  1  cm  sections. 

The  density  variation  over  the  specimen  length  was 
determined  by  immersing  a  piece  of  each  1  cm  section  in  a 
known  volume  of  mercury,  and  measuring  the  volume  of 
mercury  the  specimen  displaced.  Moisture  content  was 
determined  according  to  the  standard  method  ( ASTM 
D-2216-80).  The  results  of  the  study  are  presented  in 
figures  P.l  through  F.6. 


I 

T 


Fiqure  F.2.  Density  variation  for  specimen 


Moisture  variation  for  specimen 


Figure  F.6.  Moisture  variation  for  specimen 
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REDUCE  CODE 
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e 


0 


average  specimen! 

strain-rate 

•train 

■trass 

interface  particls  velocity 
intarfscs  strass  difference 


input  parameters 
time 
strass 
strain 
strain-rate 
interface  stresses 
interface  particle  velocities 
interface  stress  difference 


r, 
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G.2  Program  Listing 


c*  * 

C*  REDUCE  CODE  LISTING  * 

C*  * 

C 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


This  program  is  used  to  process  the  data  obtained 
from  a  split-Hopkinson  pressure  bar  experiment.  Two 
sources  of  input  data  are  required. 


1)  The  incident  and  transmitter  bar  strain  gauge 
records  hbxxxi  and  hbxxxt  and  the  calibration 
records  calxxi  and  calxxt.  These  files  will  be 
read  in  from  storage. 

2)  The  experiment  number,  gain  settings,  incident 
wave  start  time,  initial  specimen  length,  and  the 
number  of  Fourier  terms  to  be  used  in  the 
computations.  These  data  are  supplied  manually  by 
the  operator. 


The  output  of  the  program  is  the  average  specimen 
strain-rate-time,  strain-time,  and  stress-time 
histories.  In  addition,  the  particle  velocities  at 
specimen-bar  interfaces  and  stress  difference  between 
the  interfaces  are  computed.  These  data  are  formated 
for  output  to  a  printer.  The  results  are  also  plotted 
using  the  DISSPLA  9.0  graphics  package.  The 
graphic  subroutines  and  their  calls  have  been  removed 
from  this  listing. 


c* 

c* 


MAIN  PROGRAM  PARAMETER  IDENTIFICATION 


* 

* 

* 


c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


***********  ARRAYS  *********** 


call  -  incident  bar  strain  gauge  calibration  storage 
file 

calt  -  transmitter  bar  strain  gauge  calibration 
storage  file 

hbi  -  incident  bar  strain  gauge  data  storage  file 
hbt  -  transmitter  bar  strain  gauge  data  storage  file 
barinc  -  incident  and  reflected  wave  stress-time 
history 


v' 

I 


S 

I 


I 


i 


s» 

I 
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c  trsbar  -  transmitted  wave  stress-time  history 
c  diff  -  incident  stress  -  ( -  reflected  stress) 

c  sum  -  incident  stress  +  (-  reflected  stress) 

c  xint  -  average  specimen  strain-time  history 

c  xrate  -  average  specimen  strain-rate-time  history 

c  xstres  -  average  specimen  stress-time  history 

c  barref  -  reflected  wave  stress-time  history 

c  velin  -  particle  velocity  at  specimen-incident  bar 
c  interface 

c  velout  -  particle  velocity  at  specimen-transmitter  bar 
c  interface 

c  stsdel  -  interface  stress  difference 
c 

C  ***********  OUTPUT  FILES  *********** 

C 

c  hopout  -  output  data 

c  hopxxi  -  uncorrected  incident  wave  stress-time  history 
c  hopxxf  -  uncorrected  reflected  wave  stress-time  history 
c  hopxxt  -  uncorrected  transmitted  wave  stress-time 
c  history 

c  hopxxif  -  corrected  incident  wave  stress-time  history 
c  hopxxrf  -  corrected  reflected  wave  stress-time  history 
c  hopxxtf  -  corrected  transmitted  wave  stress-time  history 
c  straxx  -  strain-time  history 

c  strtxx  -  strain-rate-time  history 

c  straxx  -  stress-time  response 

c  sssrxx  -  stress-strain  response 

c 

c  ***********  SYMBOLS  *********** 

c 

c  gaini  -  gain  for  incident  bar  strain  gauge 
c  gaint  -  gain  for  transmitter  bar  strain  gauge 

c  sttime  -  incident  wave  start  time  in  microseconds 

c  srtime  -  reflected  wave  start  time  in  microseconds 
c  tstime  -  transmitted  wave  start  time  in  microseconds 
c  test  -  experiment  identification 

c  expnam  -  output  experiment  identification 

c  scfaci  -  incident  bar  strain  gauge  scale  factor 
c  scfact  -  transmitter  bar  strain  gauge  scale  factor 
c  slen  -  Initial  specimen  length 

c  mm  -  number  of  Fourier  terms 

c  velocb  -  pressure  bar  rod  velocity 

c  densb  -  density  of  pressure  bar 

c  dia  -  diameter  of  pressure  bar 

c  areab  -  area  of  pressure  bar 

c  dzi  -  distance  from  strain  gauge  to 

c  specimen-incident  bar  interface 

c  dzr  -  (-)  dzi 

c  dzt  -  distance  from  strain  gauge  to 

c  specimen-transmitter  bar  interface 

c  tstep  -  data  collection  interval 

c  ndp  -  number  of  data  points  needed  for  computations 

c  tc  -  1/tstep 


c  deltat  -  2  times  the  wave  travel  time  between  the 
c  strain  gauge  and  specimen-bar  interface 

c  ibeg  -  incident  wave  start  time  in  data 

c  ibegr  -  reflected  wave  start  time  in  data 

c  ibegt  -  transmitted  wave  start  time  in  data 

c  nlbi  -  incident  wave  baseline 

c  nlbr  -  reflected  wave  baseline 

c  nlbt  -  transmitted  wave  baseline 

c  ibegi  -  ibeg  -  nlbi 

c  ibegrr  -  ibegr  -  nlbr 

c  ibegtt  -  ibegt  -  nlbt 

c  lend  -  incident  wave  stop  time  in  data 

c  iendr  -  reflected  wave  stop  time  in  data 

c  iendt  -  transmitted  wave  stop  time  in  data 

c  tip  -  duration  of  incident  wave 

c  trp  -  duration  of  reflected  wave 

c  ttp  -  duration  of  transmitted  wave 

c  ni  -  period  of  incident  wave  for  dispersion 

c  correction 

c  nr  -  period  of  reflected  wave  for  dispersion 

c  correction 

c  nt  -  period  of  transmitted  wave  for  dispersion 

c  correction 

c  const  -  integration  constant 

c 

C  ***********  SUBROUTINES  *********** 

C 

c  cal  -  computes  scale  factors 

c  basadj  -  zeros  wave  baseline 

c  four  -  performs  dispersion  correction 

c 

£********************************************************** 

C*  MAIN  PROGRAM  * 

CAM******************************************************** 

C 

program  xx(  tty ,  input=*tty  ,output»tty  ,cali  ,calt, 
lhbi ,  hbt , tapell«cali , tapel2»calt, tapel3»hbi , 
ltapel4»hbt, tape99 ) 
c 

dimension  calt(4000),cali(4000) , barinc( 25000) , 
ltrsbar ( 25000 ) /diff(2500) , sum( 2500) ,xint(2500) , 
lxrateC  2500 ) , xst res  <  2500 ) , bar  ref ( 25  000 ) , velin( 2500 ) , 
lvelout( 2500) ,stsdel( 2500 ) 
c 

1003  £ormat(a) 

1006  format ( f 12 . 4 / lh , £12 .6 ) 

1007  formate* 1  ****  hopkinson  bar  data  reduction  code 
loutput  ****') 

1009  formate 12{ lx, f 9 . 4 ) ) 

1010  formate'l  ****  processed  hopkinson  bar  data  ****•) 

1011  format( f 10 . 1, 8x, f 5.1, 3x, f 4. 1, 3x,f 6 . 1, 4x, f 5 . 1, 
13x,f4.1,3x,f5.1/2x,f4.1/2x,f5.1) 

character* 20  hopxxi ,  hopxxr,  hopxxt,  hopxxif, 


lhopxxrf,  hopxxt f , straxx,  strtxx,  strsxx,  sssrxx, 
lexpnam,  hopout 
character* 4  test 


***********  operator  supplied  input  *********** 

print* , ' input  experiment  identification  number’ 

read  1003 /test 

print* 

print* , ' gain  for  the  incident  and  transmitter  bar 
lstrain  gauges' 
read* ,gaini,gaint 
print* 

print* , ' initial  start  time  for  incident  pulse 
1' microseconds) ' 
read*, sttime 
print* 

print* ,' initial  specimen  length  (m)' 

read* ,slen 

print* 

print* , ' input  the  number  of  terms  to  be  used  in' 
print*,' the  Fourier  analysis' 
read*, mm 

********************************************************** 
*******  aet  file  names  names  for  storage  ******* 

hopxxi- ' hp ' //tea t// ' i ' 
hopxxr- ' hp ' //test// ' r ' 
hopxxt- ' hp ' //test// ' t ' 
hopxxif-'hp'//te8t//' if' 
hopxxrf- ; hp ' //test// ' rf ' 
hopxxt f-'hp' //test//' tf ' 
straxx- ' str ' //test// ' time ' 
strtxx- ' rat ' //test// ' time ' 
a trsxx- ' sts ' //test// ' time ' 
sssrxx- ' prs ' //test// ' strain ' 
expnam- ' hb ' //test// ' x ' 
hopout- ' hp ' //test// ' ou t ' 

********************************************************** 
**********  aet  parameter  values  ********** 

velocb-4886. 
densb-8090 . 
dia-0. 06033 

areab-( 3 . 141?9* (dia**2 ) )/4.0 

dzi-0.6033 

dzr— dz  i 

dzt-dzr 

tstep-0 . 5 

ndp-15000 


o  n  n  o  noon  on 
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tip<*240 

trp-tip 

ttp*tip 

deltat»2 .*dzi/velocb 
srtime»attirae+deltat*le+6 
tstirae-srtime 
tc«l/tstep 
c 

c** ********************************************** ********** 
c  *******  reaa  in  the  calibration  and  data  files  ******* 
c 

readdl,*)  (cali(  i) ,  i-1,4000) 
read(12/*) (calt( i) ,i-l/4000> 
read(13, *) (barinc( i) , i-l,ndp) 
read (14,*) ( trsbar ( i) , i-l,ndp) 
c 

c********************************************************** 

c  ***  calculate  the  scale  factor  from  calibration  data  *** 
c 

call  caKcali, gaini,acfaci , 4000) 
call  cal(calt,gaint,scfact,4000) 
c 

c* ********************************************************* 
c  ****  compute  starting  time  of  waves  in  the  data  **** 
c 

ibeg«tc*sttime+l 

ibegr*tc*srtime+l 

ibegt»tc*tstirae+l 

c 

c** ******** **************************************** ******** 
*****  compute  stop  times  of  waves  in  the  data***** 

iend-tc* ( sttime+tip) +1 
iendr»tc* ( srtime+trp) +1 
iendt»tc*( tstirae+ttp)+l 


*******  length  of  baseline  for  each  wave  ******* 

nlbi»50 
nlbr-20 
nlbt«50 

******************************************************** 
********  compute  baseline  start-times  ******** 

ibegi«ibeg-nlbi 
ibegrr«ibegr-nlbr 
ibegtt»ibegt-nlbt 
c 

c********************************************************** 
c  *******  zero  wave  baseline  ******** 

,  c 


call  baead j ( bar inc , scf aci , ibegi , iendr , nlbi , base) 
call  basad j ( trsbar , scf act, ibegtt, iendt ,nlbt, base) 


c  **  creating  a  separate  array  to  store  reflected  pulse  ** 
c 

do  25  intend, iendr 
bar ref ( i)-barinc( 1) 

25  continue 
c 

c**«* ********************************************  ********* 
c  *******  store  raw  data  ******* 

xi«( ibeg-1) *tstep 
xr-' ibegrr-1) *tstep 
xt«v  _oegtt-l) * tstep 

open( unit-15 , file-hopxxi) 
opent unit-16 , f ile-hopxxr) 
open(unit-17 ,  f ile-hopxxt) 

do  30  i-ibegi,iend 
writet 15 , 1006 ) xi /bar inc ( i)/le+6 
xi-xi+tstep 
30  continue 

do  35  i-ibegrr, iendr 
write(16,1006)xr,barref ( i)/le+6 
xr-xr+tstep 
35  continue 

do  40  i-ibegtt, iendt 
write(17,1006)xt,trsbar( i)/le+6 
xt-xt+tstep 
40  continue 

closet  15 ) 
closet  16 ) 
closet  17 ) 

********************************************************* 
c  *  perform  dispersion  correction  on  each  wave  separatly 
c 

nr«( iendr-iend+50 ) *2 

call  four ( bar ref , tstep, dia , velocb, nr ,dzr , iend , mm) 
c 

ni-t iend- ibegi ) *2 

call  four (bar inc /tstep, dia, velocb /ni,dzi , ibegi, mm) 


c 


nt-t iendt-ibegtt) *2 

call  fourt  trsbar , tstep, dia, velocb , nt,dzt, ibegtt, mm) 
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c  *******  store  dispersion  correction  data  ******* 

c 

42  xi«( ibegi-1) *tstep 
xr-( ibegrr~l)*tstep 
xt*( ibegtt-1) *tstep 
c 

open ( unit-18 , f ile-hopxxif ) 
open (unit-19 , f ile-hopxxrf ) 

.  open ( unit-20, f ila-hopxxtf ) 
c 

do  45  i-ibegi,iend 
write(18,1006)xi,barinct i)/le+6 
xi-xi+tstep 
45  continue 
c 

do  47  i«ibegrr , iendr 
w.itetl9,1006)xr,barref ( i)/le+6 
xr-xr+tstep 

47  continue 
c 

do  48  i-ibegtt, iendt 
writet20,1006)xt,trsbar( i)/le+6 
xt-xt+tstep 

48  continue 
c 

closet  18) 
closet  19) 
closet  20) 
c 

£************************************ ****** **************** 
c  *****  compute  the  following i  ***** 
c  1)  strain, 

c  2)  strain-rate, 

c  3)  stress, 

c  4)  particle  velocity, 

c  5)  interface  stress  difference 

c 

j-ibegr 

ib-1 

do  50  i-ibeg,iend 
diff ( ib)-barinc( i) -bar reft j) 
sumt ib)-barinc( i)+barref ( j) 

5-j+l 
Ib-ib+1 
50  continue 
c 

const-l/t  slen*densb*velocb) 
c 

ipo int- lend- i beg 
xintt 1 )-0 .0 
xrate( l)-0.0 
tstepi-tstep*l . 0e-6 
ic-2 
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j-ibegt+1 

do  100  i-ic,ipoint 
n-i-1 
k-j-1 

xrate( i)-(dif£ (n)-trsbar(k) )*const 

xint( i)-xint( n)  + ( tstepi*0 . 5*(xrate( i)+xrate( n ) ) ) 

j-j+1 

100  continue 
c 

j-ibegt 

do  120  i-l,ipoint 
xstreat i)»( sumt i)+trsbar( j) )/2.0 
j-j+1 

120  continue 
c 

do  127  i-l,ipoint 
n-i-1 

velin( i)-( (diff ( i) )/(densb*velocb) ) 
velout( i)»( ( trsbar ( ibegt+n) )/(densb*velocb) ) 

127  continue 
c 

do  128  i-l,ipoint 
n-i-1 

stadeK  i)«(sum(  i)  -trabar(  ibegt+n)  )/le+6 

128  continue 
c 

C*******************************4************************* 

c  *******  write  output  files  ******* 

c 

open( unit-21, file-straxx) 
open ( uni t-22 , f ile-atrtxx) 
open (unit-23, f ile-strsxx) 
open (unit-24, file-saarxx) 
t-0.0 

do  130  i-l,ipoint 
xstrea( i)-xatres( i)/le+6 
xint( i)-xint( i) *100 . 
write( 21,1006) t,xint< i) 
write( 22,1006) t,xrate( i) 
write (23, 1006) t,xatres( i) 
write(24,1006)xint( 1) ,xstres( i) 
t-t+tatep 
130  continue 
closet  21) 
closet  22) 
closet  23 ) 
closet  24 ) 
c 

C********************************************************** 

C  *******  output  ******* 

c 

xi«( ibeg-1) *tstep 
xr-( ibegr-1) *tstep 


i 


i 


i 


xt«( ibegt-1 ) * tstep 

open ( uni t»27,f ile-hopout) 
write(27,1007) 
write( 27,* ) '  ' 
write(27,*) '  • 
write( 27,* ) ’  • 

write(27,*)'  ***  input  parameters  ***' 

write(27,*)'  • 
write(27,*)'  ' 

writet 27 ,*)' experimental  identification  number 

1 - -  ' , expnam 

write( 27 , *  ) '  ' 

write( 27 , * ) 'gain  for  incident  bar  strain  gauge 

1 -  ' , gaini 

write( 27 , * ) 1  ' 

writet 27 , * ) ' gaint  for  transmitter  bar  strain  gauge 

1 -  ' , gaint 

write(27/*)'  ' 

write( 27 ,*)' incident  wave  starttime  (microseconds) 

1 -  ' ,  xi 

write(27#*)'  • 

write(27/*) 'reflected  wave  starttime  (microseconds) 

1  t 

write(27,*)'  ' 

write( 27,*) 'transmitted  wave  starttime  (microseconds) 
1-  * ,xt 

write( 27  ,* ) '  ' 

write<27,*) 'initial  specimen  length  (m) 

1 - - - -  '  ,  alen 

write(27,*)'  ' 

write( 27 ,*)' number  of  terms  used  in  fourier  analysis 
1 —  ' ,mm 

write(27,*) '  ' 

write( 27 ,*) 'dispersion  bar  correction  length(m) 

1 -  '  ,dzi 

write(27,*)'  ' 

write(27,* )' incident  bar  calibration  factor 

1 -  ',8cfaci 

write(27,*)'  ' 

write( 27 ,*)' transmitter  bar  calibration  factor 
1 -  ' ,scfact 

write( 27 ,1010) 
write(27,*)'  ' 
write(27,*)'  ' 
write( 27 , * ) '  ' 

write(27,*)'  time  stress  strain  strainrate 

linfll  veil  infl2  vel2  stsdif* 
write( 27 ,*)' (microseconds )  (mpa)  (%)  (1/sec) 

l(mpa)  (m/s)  (mpa)  (m/s)  (mpa)' 
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ti-0 

do  137  i-l,ipoint 
n»i-l 

write(27,1011) ti,xstres< i) , xintt 1) , xrate( i) , 
lsum( i)/le+6,velin( i) , trsbar ( ibegt+n ) /le+6 , 
lveloutt i) , atadelt i) 
ti-ti+tstep 
137  continue 
c 

closet  27) 
c 

14  end 

c 

c* ************** ******************************************* 
q********************************************************** 

c*  END  MAIN  PROGRAM  * 

C********************************************************** 

C 

c 

c ********************************************************** 
C*  SUBROUTINES  * 

c* ********************************************************* 

c 

c 

c  This  subroutine  is  used  to  compute  scale  factors 

c  that  are  used  to  transform  the  digitized  numbers 

c  to  stress-time  data 
c 
c 

subroutine  calt blk, gain, scl£ac,ndp) 

dimension  blk(400C) 

calmax»0.0 

calzer«0.0 

ncalmx«0 

ncalze*0 

ibeg-1 

if (blk( 1) .gt . 164 )  goto  130 
100  do  105  i»ibeg+l , ndp 

if ( blk( i) .gt . 164 )  then 
iend-i 
goto  110 
end  if 

105  continue 
iend*ndp 

110  if ( ( iend-ibeg) .gt.8)  then 
do  115  i»ibeg+4 , iend-4 
calzer*calzer+blk( i } 

115  continue 

ncalze«ncalze+ iend-ibeg -7 
endif 
i beg* i end 

if ( ibeg .ge.ndp)  goto  190 
130  do  135  i»ibeg+ 1 , ndp 


o  o 
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i£(blk( i) .It. 164)  then 
iend»i 
goto  140 
endif 

135  continue 
iend^ndp 

140  if ( ( iend-ibeg) .gt .8 )  then 
do  145  i«ibeg+4, iend-4 
calraax«calmax+blk( 1) 

145  continue 

ncalmx»ncalmx+iend-ibeg-7 

endif 

ibeg^iend 

if( ibeg .ge. ndp)  goto  190 
goto  100 
c 

190  calmax»calroax/ncalmx 
calzer»calzer/ncalze 
self ac«l . 2611e+9/gain/( calmax-calzer ) 
c 

return 

end 


c* ********************************************************* 

c 

c 

c  This  subroutine  perform  the  dispersion  correction 

c 

c 

subroutine  £our( fct,tstep,dia,velocb,nn,dz, 

1 ibeg in, mm) 
c 

c  *********  SYMBOLS  ********** 

C 

c  fct(i)  -  array  holding  wave  data 

c  nn  ■  period  of  the  wave  being  analyzed 

c  a  ■  array  to  store  An  Fourier  coefficient  values 

c  b  ■  array  to  store  Bn  Fourier  coefficient  values 

c  t  »  array  to  store  phase  angle  values 

c  c  ■  array  that  holds  the  phase  velocity  at  each 

c  frequency 

c  equa(i)  »  dummy  array  for  computational  purposes 

c  v  ■  bar  radius/wavelength 

c  wo  •  wave  frequency 

c  ao  ■  value  of  Ao  Fourier  coefficient 

c  d  •  array  to  store  the  value  of  Cn  after 

c  transformation  to  a  cosine  series 

c  co  ■  longitudinal  wave  velocity  of  bar 

c 

c  *******  SUBROUTINES  ******* 

c 

c  coeff  -  computes  the  Fourier  coefficients  Ao,  An, 


'  .*  *  * Ajl* i'iA.'i’. 


o  o  o  •-* 
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c  and  Bn 

c  simp  -  performs  Simpson'  rule  integration 

c 

c 

dimension  fct( 25000) , a(100) ,b(100) , c(100) , 
ld(100) , t( 100) 
c 

c  function  statement;  phase  velocity  verses  r/lambda 
c 

Cp(v>-.5764+< .4236/(22. *v**4+12 .8*v**3-2 . 77*v**2+ 
1.92*v**l. 5+1.0) ) 
c 

freq«1.0/( tstep*l . Oe-6 ) 
c 

do  20  j«l,nn 
n-j-1 

fct( j)«fct( ibegin+n) 

20  continue 
pi-3.14159 
wo«2*pi/nn 
c 

call  coeff (mm,nn,wo, fct ,ao,a,b) 
c 

do  60  j-l,mm 

d( j)-sqrt(a( j) **2+b( j)**2) 
t( j)-aco8( a( j )/d( j ) ) 
if (b( j)-0.0)  55,60,60 
55  t(  j)—  1.0*t(  j) 

60  continue 
c 

69  w-(dia/2.0) *wo*f req/velocb 

v-0.0 
i-0 

70  i-i+1 
x»i*w 

80  v-v+0.0001 

xp»2.0*pi*cp(v)*v 
83  if(xp-x)  80,85,85 
85  c(i)-cp(v) 

88  if ( i-mm)70,90,90 
c 

90  do  100  i-l,mm 

c( i)-c( i)*velocb 
100  continue 
c 

do  120  i-1 , mm 

ak«i*wo/velocb* ( velocb/c( i ) -1 . 0 ) *f req 
t( i)-t( i)+ak*dz 
120  continue 

perform  wave  reconstruction 
do  150  i-l,nn 


; 
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fct< i)-ao/2.0 
do  150  j-1, mm 

£ct( i)-fct( i)+d( j)*cos( ( j*wo*i)-t( j) ) 

150  continue 
c 

do  160  j-l,nn 
n«  j-1 

fct( ibegin+n)-fct< j) 

160  continue 
c 

return 

end 

c 

c 

c* *************************  **************************  ****** 

c 

c 

subroutine  coeff  (mm,nn,wo,  fct,ao,a,b) 
c 

c  this  routine  will  compute  the  Fourier  coefficients 
c  for  each  term  in  the  series 
c 

dimension  equa(6000) ,a<100) ,b(100) ,fct(25000) 
c 

c  computing  the  integration  limit;  t/2 

c 

n-nn/2 

do  10  i-l,nn 
equal i)»fct( i) 

10  continue 
c 

c  Simpson's  rule  integration  will  be  performed  to 
c  evaluate  each  of  the  coefficients;  Ao,An,Bn 

c 

c  computing  Ao 

c 

call  simp{equa,nn,result) 
c 

ao-result/n 

c 

c  compute  An 

c 

do  30  j-1, mm 
do  20  i-l,nn 

equa(i)-fct( i)*cos( i*j*wo) 

20  continue 

call  simp(equa,nn,result) 
a( j )-result/n 
30  continue 
c 

c  compute  Bn 

c 

do  50  j-1, mm 


aa^aaa^s tfs 
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do  40  1*1, nn 

equa< i)»fct( i) *sin( i* j*wo) 
40  continue 

call  a imp ( aqua ,nn, result) 
b( j)-result/n 
SO  continue 
c 

return 

end 

c 

c 


c 

c 

subroutine  simp<equa,nn,result) 
dimension  f ( 6000 ), equa ( 6000 ) 
c 

do  2  i-l,nn 
f ( i)»equa( i) 

2  continue 
c 

n-nn/2 

c 

c  this  subroutine  performs  simpson's  rule  integration 
c 

c  parameters 

c 

c  f:  integration  points 

c  n:  number  of  data  points 

c  h:  spacing  between  data  points 

cresult:  estimate  of  integral  returned  to  caller 

c 

c  check  to  see  it  the  number  of  panels  is  even 
c 

c  n-1:  number  of  panels 

c 

npanel-n-1 
nhalf»npanel/2 
nbegin-1 
result-0. 0 
h-1 

if ( (npanel-2*nhalf ) .eq.0)  go  to  5 
c 

c  if  the  number  of  panels  is  odd  use  3/8  rule  for  the 
c  first  four  points  then  continue  using  the  1/3  rule 

c  if  the  numbers  of  panels  is  even  then  the  1/3  rule 

c  will  be  used  exclusively 

c 

result-3. 0*h/8.0*(f (l)+3.0*f (2)  +  3.0*f <  3)+f (4) ) 
nbegin-4 
c 

c  apply  1/3  rule,  add  in  the  first , second,  and  last 
c  values 


Vl  I  i  t  '■  I 
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c 

5  result«result+h/3 . 0* ( f ( nbeg in ) + 

14 .  0*f  (nb«»gin+l)+f  (n) ) 
nbegin-nbegin+i 
if (nbegin.eq.n)  go  to  15 
c 

c  the  pattern  after  nbegin+2  is  repetitive 
c 

nend«n-2 

do  10  i«nbegin,nend,2 
result»reault+h/3 .0*1 2 . 0*f ( i)+4.0*f ( i+1) ) 

10  continue 
c 

15  return 
end 
c 
c 

£********************************************************* 

c 

c 

subroutine  basad j { blk , sc 1 fac , nbeg , nend , n , base ) 
dimension  blk( 25000) 
c 

50  k«nbeg+0 .8*f loat(n) 
base»0.0 
count»0 . 0 
do  100  i*nbeg,k 
count»count+l 
base=*base+blk(  i) 

100  continue 
c 

base»base/count 

c 

c  adjust  data  for  baseline  and  apply  scalefactor 
c 

150  do  200  i»nbeg,nend 

blk( i)»(blk( i ) -base ) *sclf ac 
200  continue 
c 

return 

end 

c 


c 


C*  END  * 


xxxx  hopkinson  bar  data  raduction  coda  output  xxxx 


XXX  input  parametars  xxx 

oxparimantal  idantifieation  numbar  hbl34ax 

gain  for  incidant  bar  strain  gauga - -  2.5 

gaint  for  tranamittar  bar  strain  gauga  -  5. 

incidant  wava  atarttima  (microseconds)  2519. 

raflactad  wava  starttima  (microsaconds)  —  2765.5 

transmitted  wava  starttima  (microsaconds)  ~  2765.5 

initial  spaciman  length  (m) - — -  1 . 26 e-2 

number  of  terms  used  in  fourier  analysis  —  20 

disparsion  bar  correction  langth(m)  - -  0.6033 

incident  bar  calibration  factor  -  7630275.099924 

transmitter  bar  calibration  factor  - -  3816424.365323 
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XXXX  processed  hopkinson  bar  data  xxxx 


tin# 

•trass 

strain 

•trainrata 

inf  1 1 

vail 

inf  12 

val2 

stsdif 

( mi crosaconds) 

(npa ) 

(54) 

( 1/sac) 

(mpa ) 

(m/s) 

(mpa) 

(m/a) 

(mpa) 

0.0 

10.2 

0.0 

0.0 

21.7 

-0.3 

-1.4 

0.0 

23.1 

0.5 

10.2 

0.0 

-21.1 

21.6 

-0.3 

-1.3 

0.0 

22.9 

1.0 

10.1 

0.0 

-18.2 

21.5 

-0.2 

-1.2 

0.0 

22.7 

1.5 

10.1 

0.0 

-14.9 

21.3 

-0.2 

-1.1 

0.0 

22.4 

2.0 

10.0 

0.0 

-11.1 

21.0 

-0.1 

-1.0 

0.0 

22.1 

2.5 

9.9 

0.0 

-6.8 

20.7 

0.0 

-0.9 

0.0 

21.7 

3.0 

9.7 

0.0 

-2.1 

20.4 

0.0 

-0.9 

0.0 

21.2 

3.5 

9.6 

0.0 

3.2 

19.9 

0.1 

-0.8 

0.0 

20.7 

4 . 0 

9.4 

0.0 

9.1 

19.5 

0.2 

-0.7 

0.0 

20.2 

4.5 

9.1 

0.0 

15.5 

18.9 

0.3 

-0.7 

0.0 

19.6 

5.0 

8.9 

0.0 

22.4 

18.4 

0.4 

-0.6 

0.0 

19.0 

5.5 

3.6 

0.0 

30.0 

17.8 

0.5 

-0.6 

0.0 

18.4 

6 . 0 

8.3 

0.0 

38.2 

17.2 

0.6 

-0.5 

0.0 

17.7 

6.5 

8.0 

0.0 

47.0 

16.5 

0.7 

-0.5 

0.0 

17.0 

7.0 

7.7 

0.0 

56.5 

15.9 

0.8 

-0.5 

0.0 

16.4 

7.5 

7.4 

0.0 

66  .6 

15.2 

1.0 

-0.5 

0.0 

15.7 

0.0 

7.0 

0.0 

77.5 

14.5 

1.1 

-0.5 

0.0 

15.0 

8.5 

6.7 

0.0 

89.0 

13.3 

1.3 

-0.5 

0.0 

14.4 

9.0 

6.3 

0.0 

101.1 

13.2 

l  4 

-0.5 

0.0 

13.7 

9.5 

6.0 

0.0 

114.1 

•  3.5 

1  .6 

-0.6 

0.0 

13.0 

10.0 

5.6 

0.0 

127.7 

11.8 

1  2 

-O.u 

0.0 

12.4 

10.5 

5.3 

o.c 

142.0 

11.2 

2.  ) 

-0.6 

0.0 

11.8 

11.0 

4.9 

0.0 

157.1 

10.5 

2 .  _ 

-0.7 

0.0 

11.2 

11.5 

4.6 

0.1 

173.0 

9.9 

2.4 

-0.7 

0.0 

10.7 

12.0 

4.3 

0.1 

189.6 

9.3 

2.6 

-0.8 

0.0 

10.1 

12.5 

4.0 

0.1 

206.9 

8.8 

2.8 

-0.8 

0.0 

9.6 

13.0 

3.7 

0.1 

225.0 

8.3 

3.0 

-0.9 

0.0 

9.2 

13.5 

3.4 

0.1 

243.8 

7.8 

3.3 

-0.9 

0.0 

8.7 

14.0 

3.2 

0.1 

263.4 

7.3 

3.6 

-1.0 

0.0 

8.3 

14.5 

2.9 

0.1 

283.8 

6.9 

3.8 

-1.0 

0.0 

8.0 

15.0 

2.7 

0.1 
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APPENDIX  I 

PLOTTED  EXPERIMENTAL  RESULTS 
WITH  SUMMARY  DATA 
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The  following  nomenclature  has  been  used  in  the 
presents ti in  of  the  summary  data: 

Oa  -  magnitude  of  appli  ed  stress  wave, 

C max  ■  maximum  average  strain-rate, 

Cavg  ■  average  strain-rate  over  the  duration  of  the 
experiment, 

1  a  specimen  length  after  seating  has  been  applied, 
£$*  seating  strain, 

w  a  initial  specimen  moisture  content, 

Wj  a  percentage  if  initial  moisture  lost  during  the 
experiment, 

e  »  void  ratio  based  on  specimen  length  1, 
ng  *  initial  gas  porosity, 

Vt  «  specimen  total  volume, 

7t  a  wet  density, 

S  *  degree  of  saturation. 
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Figure  1.1.  Stress-strain  response  for  experiment  31 
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Figure  1.2.  Stress-strain  response  for  experiment  32 


Figure  1.3.  Stress-strain  response  for  experiment  38 
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Figure  1.5.  Stress-strain  response  for  experiment  40 
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Stress-strain  response  for  experiment  42 
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Figure  1.9.  Stress-strain  response  for  experiment  56 
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Figure  I. 10.  Stress-strain  response  for  experiment  57 
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Figure  I. 11.  Stress-strain  response  for  experiment  58 
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Figure  1.14.  Stress- stra in  response  for  experiment  61 
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Figure  1.15 
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Figure  1.16.  Stress-strain  response  for  experiment  71 
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Figure  1.18.  Stress-strain  response  for  experiment  73 
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Figure  1.20.  Stress-strain  response  for  experiment  75 
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Figure  1.21.  Stress-strain  response  for  experiment  76 
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243  MPa 


response  for  experiment  119 
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Stress-strain  response  for  experiment  132 
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Stress-strain  response  for  experiment  135 


Stress-strain  response  for  experiment 
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Figure  1.43.  Stress-strain  response  for  experiment  162 
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Figure  1.44.  Stress-strain  response  for  experiment  163 
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Figure  1.45.  Stress-strain  response  for  experiment  164 


Figure  1.46.  Stress-strain  response  for  experiment  165 
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Stress-strain  response  for  experiment  166 
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Due  to  the  uncertainty  in  the  evaluation  of  the 
strain  contributions  from  factors  discussed  in  section 
6.4,  a  gross  strain  adjustment,  to  account  for  the  excess 
between  strain  at  peak  stress,  and  the  initial  gas 
porosity  has  been  computed.  This  was  accomplished  by 
plotting  the  difference  between  the  strain  at  peak  stress 
and  the  initial  gas  porosity  against  the  average  force 
felt  by  the  specimen.  The  average  force  is  computed  as  the 
peak  stress  multiplied  by  the  area  of  the  specimen.  A 
linear  regression  line  was  then  fit  to  the  data  for 
specimens  prepared  near  the  optimum  compaction  conditions. 
The  linear  regression  line  was  computed  with  the  dependent 
variable  (Y)  taken  as  the  difference  between  the  initial 
gas  porosity  and  the  strain  at  peak  stress  and  the 
independent  variable  (X)  taken  as  the  maximum  force  felt 
by  the  specimen.  The  equation  of  the  linear  regression 
line  is; 

Y  *  0.98X  +  0.56  .  (J.l) 

The  resulting  plot  is  shown  in  figure  J.l.  The  data  used 
to  prepare  figure  J.l  are  presented  in  table  J.l. 

Using  the  regression  line  and  the  average  force 
sustained  by  the  specimen,  a  strain  correction  was 
computed  for  each  specimen.  The  strain  correction  is  added 
to  the  initial  gas  porosity  and  the  sum  compared  to  the 
strain  at  peak  stress.  The  balance  can  then  be  used  to 


Table  J.l 


Strain  Correction  Data 


Experiment  No. 

Strain  -  Gas  Porosity 

(%) 

Force 

(dynes  10E10) 

112 

11.9 

7.5 

113 

9.7 

7.2 

114 

10.4 

7.2 

115 

7.5 

7.7 

116 

6.2 

4.5 

117 

6.0 

6.1 

118 

3.8 

4.7 

119 

4.2 

4.5 

131 

3.1 

5.6 

132 

3.8 

6.0 

133 

4.5 

5.1 

134 

6.4 

6.1 

138 

6.9 

5.3 

145 

3.9 

3.0 

146 

7.0 

9.0 

147 

6.3 

4.3 

162 

16.2 

9.4 

163 

10.4 

12.2 

164 

4.9 

8.9 

165 

2.2 

2.2 

a)  Strain  is  taken  as  the  strain  at  peak  stress 


gauge  the  success  of  the  correction.  The  results  of  the 
correction  procedure  are  tabulated  in  table  J.2. 

This  procedure  requires  none  of  the  assumptions  made 
in  the  computations  of  strain  contribution  from  moisture 
loss*  pore  water  compression,  or  radial  expansion.  It  is  - 
based  solely  on  the  peak  stress  sustained  by  the  specimen. 
In  addition,  the  strain  contribution  of  soil  loss  is  also 
included. 
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Table  J.2 


Strain  Correction  Results 


H 


i 

I 


Kxparlaent 

NO. 

Oas  Porosity 

<%) 

Strain 

Correction 

(«) 

<%> 

Strain  at 

Peak  S trass 

(4) 

Till  anna 

(«) 

112 

7.91 

7.87 

15.78 

19.79 

4.01 

113 

7.73 

7.67 

15.40 

17.40 

2.00 

114 

9.21 

7.61 

16.82 

19.60 

2.78 

115 

10.55 

8.07 

18.62 

18.03 

-0.59 

116 

8.07 

4.98 

13.05 

14.24 

1.19 

117 

8.50 

6.54 

15.04 

14.50 

-0.54 

118 

10.25 

5.18 

15.43 

14.06 

-1.37 

119 

9.91 

4.93 

14.84 

14.07 

-0.77 

131 

9.57 

6.10 

15.67 

12.68 

-2.99 

132 

6.88 

6.43 

13.31 

10.67 

-2.64 

133 

8.02 

5.52 

13.54 

12.52 

-1.02 

134 

5.80 

6.57 

12.37 

12.22 

-0.15 

138 

4.86 

S.75 

10.61 

11.80 

1.19 

1*5 

6.23 

3.51 

9.74 

10.12 

0.38 

146 

6.26 

9.41 

15.67 

13.30 

-2.37 

147 

6.64 

4.78 

11.42 

12.99 

1.57 

162 

4.44 

9.74 

14.18 

20.60 

6.42 

163 

4.84 

12.52 

17.36 

15.27 

-2.09 

164 

7.49 

9.28 

16.77 

12.43 

-4.34 

165 

7.05 

2.70 

9.75 

9.27 

-0.48 

f* 

/• 

r 

t 


« 


/'/■  av.v.v  ;  ,‘.v. 
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